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NOTICE TO AUTHORS 


1. Communications.—Papers must be communicated to the Society by a Fellow. They 
should be accompanied by 2 summary at the beginning of the paper conveying briefly the 
content of the paper, and drawing attention to important new information and to the main 
conclusions. The summary should be intelligible in itself, without reference to the paper, 
to a reader with some knowledge of the subject; it should not normally exceed 200 words 
in length. Authors are requested to submit MSS. in duplicate. These should be 
as double spacing and leaving a margin of not less than one inch on the 
left-hand side. Corrections to the MSS. should be made in the text and not in the 

in. Unless a paper reaches the Secretaries more than seven days before a Council 
meeting it will not normally be considered at that meeting. By Council decision, MSS. of 
accepted papers are retained by the Society for one year after publication; unless their 
return is then requested by the author, they are destroyed. 


2. Presentation.—Authors are allowed considerable latitude, but they are requested to 
follow the general style and arrangement of Monthly Notices. References to literature 
should be given in the standard form, including a date, for printing either as footnotes or in 
a numbered list at the end of the paper. Each reference should give the name and 
initials of the author cited, irrespective of the occurrence of the name in the text (some 
latitude being permissible, however, in the case of an author referring to his own work). 
The following examples indicate the style of reference appropriate for a paper and a book, 
respectively :— 

A. Corlin, Zs. f. Astrophys., 1g, 239, 1938. 

H. Jeffreys, Theory of Probability, 2nd edn., section 5.45, p. 258, Oxford, 1948. 


3. Notation.—For technical astronomical terms, authors should conform closely to the 
recommendations of Commission 3 of the International Astronomical Union (Trans. 
1.A.U., Vol. V1, p. 345, 1938). Council has decided to adopt the I.A.U. 4-letter abbrevi- 
ations for constellations where contraction is desirable (Vol. IV, p. 221, 1932). In general 
matters, authors should follow the recommendations in Symbols, Signs and Abbreviations 
(London: Royal Society, 1951) except where these conflict with I.A.U. practice. 


4. Diagrams.—These should be designed to appear upright on the page, drawn 
about twice the size required in print and prepared for direct photographic 
reproduction except for the lettering, which should be inserted in pencil. 
Legends should be given in the manuscript indicating where in the text the 
figure should appear. Blocks are retained by the Society for 10 years; unless the author 
requires them before the end of this period they are then destroyed. 


ss 5. Tables.—These should be arranged so that they can be printed upright on 
e page. 


6. Proofs.—Costs of alterations exceeding 5 per cent of composition must be borne by 
the author. Fellows are warned that such costs have risen sharply in recent years, and it 
is in their own and the Society’s interests to seek the maximurn conciseness and simplifi- 
cation of syrabols and equations consistent with clarity. 


7. Revised Manuscripts.—When papers are submitted in revised form it is especially 
requested that they be accompanied by the original MS. 


Reading of Papers at Meetings 


_8. When submitting papers authors are requested to indicate whether they will be 
willing and able to read the paper at the next or some subsequent meeting, and approxi- 
mately how long they would like to be allotted for speaking. 


g. Postcards giving the programme of each meeting are issued some days before the 
meeting concerned. Fellows wishing to receive such cards whether for Ordinary 
Meetings or for the Geophysical Discussions or both should notify the Assistant Secretary. 
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MEETING OF 1952 MAY 9 
Professor H. Dingle, President, in the Chair 


‘The election by the Council of the following Fellows was duly confirmed :— 

William Hallett Day, 22 Nugent Road, Southbourne, Bournemouth (proposed 
by F. O’B. Ellison); 

Anton Linder Hales, B.A., Ph.D., Department of Mathematics, University 
of Cape Town, Rondebosch, South Africa (proposed by W. S. Finsen); 

George Anthony Hole, 44 Sanyhils Avenue, Patcham, Brighton, Sussex 
(proposed by F. J. Hargreaves); 

Reginald Walter Howers, 7 Aboyne Drive, Raynes Park, London, S.W. 20 
(proposed by G. C. Saul); 

Jack Thurston Kent, Department of Mathematics, Agricultural and Mechan- 
ical College of ‘Texas, College Station, ‘Texas, U.S.A. (proposed by R. H, 
Garstang); and 

Albert E. Whitford, Director, Washburn Observatory, University of 
Wisconsin, Madison 6, Wisconsin, U.S.A. (proposed by S. Chandrasekhar). 

One hundred and ten presents were announced as having been received 

since the last meeting. 


ADDITIONAL MEETING OF 1952 JULY 25 


in the Physics Lecture Theatre of the 
University of Leeds 


Professor H. Dingle, President, in the Chair 


The President announced that this Additional Meeting was the fifth to take 
place outside the rooms of the Society in London. He expressed the appre- 
ciation of the Society to the Vice-Chancellor and Council of the University of 
Leeds, through whose courtesy, and that of Professor E. C. Stoner, F.R.S., 
the meeting was being held in the Physics Lecture Theatre. He also thanked 
the Warden of Devonshire Hall, to whom the Society was indebted for the 


34 
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excellent accommodation provided at the Hall; and the President and Council 
of the Leeds Astronomical Society, who as leaders of the astronomical life of 
the City were also the hosts of the Society. ‘The thanks of the Society were 
also accorded to Professor Cowling, who had been responsible for the organi- 
zation of the Additional Meeting and conjoint activities. 


Ninety-six presents were announced as having been received since the last 
meeting, including :-— 
G. D. Rochester and J. G. Wilson, Cloud chamber photographs of cosmic 
radiation (presented by the authors). 


A number of Fellows signed the Admission Book. Papers were read by 
Mr J. Hargreaves, Mr R. Hanbury Brown, Mr L. Mestel and Professor T. G. 
Cowling. 

Other meetings were held in connection with the Society’s visit to Leeds. 
On the evening of July 23, the Vice-Chancellor received the Fellows in the 
Central Court of the Parkinson Building. An Astronomical Colloquium on 
the total solar eclipse of 1952 February 25 was held on July 24, and, in the 
evening, a Public Lecture, entitled Radio Astronomy, was given by Professor 
A. C. B. Lovell. A Geophysical Discussion was held on July 26. Visits to 
places of local interest included an excursion to Fountains Abbey; and another 
to York, where the Fellows had the opportunity of visiting either the Minster 
(through the courtesy of the Dean) and the Yorkshire Museum (under the 
guidance of Mr Sidney Melmore); or the optical works of Messrs Cooke, 
Troughton and Simms, where they were welcomed by Mr E. Wilfred Tayler. 





A STATISTICAL BASIS FOR THE THEORY OI 
STELLAR SCINTILLATION 


S. Chandrasekhar 
(Received 1952 February 26) 


Summary 


Current theories of stellar scintillation and astronomical seeing suppose 
that there is a disturbed region in the atmosphere at a height of about 4 km 
which corrugates a plane wave-front passing through it; and that the observed 
phenomena are to be accounted for in terms of such a corrugated wave-front. 
On the assumption that the disturbed region is a turbulent layer in which the 
refractive index, 4, is subject to irregular fluctuations, the auto-correlation, 
du(r,)Ou(r2)/d24, of the instantaneous fluctuations in the refractive index at 
iwo different points r, and r, is introduced; on the further assumption that 
the turbulence which prevails is homogeneous and isotropic, the auto- 
correlation so defined can be a function, M(r) (say), only of the distance r 
between the two points considered. It is then shown how the statistical 
properties of the corrugated emergent wave-front, such as the angular 
dispersion in the wave normals, can be expressed in terms of M(r). From 
the known facts concerning astronomical seeing it is concluded that we can 
satisfactorily account for the observed phenomena by postulating a turbulent 
layer of a thickness of the order of a hundred metres, a micro-scale of tur- 
buience of the order of ten centimetres and a root mean square fluctuation in 
refractive index of the order of 4 x 107°. 





1. Introduction.—On examining the recent literature on the subject of stellar 
scintillation and astronomical seeing, it appeared to the writer that a proper 
statistical basis for discussing certain aspects of the phenomenon was lacking. 
In this paper an attempt will be made to provide such a basis. But it may be 
useful, first, to clarify the nature of the problem and to indicate the need for a 
theory along the lines outlined in this paper. 

It would appear that in discussing the general problem of stellar scintillation 
we should distinguish between the phenomena which are observed near the 
horizon and at low altitudes and the phenomena which are observed near the 
zenith and at high altitudes. ‘That such a distinction should be made is clear 
already from Lord Rayleigh’s (x) discussion of the problem in 1893. For, as 
Rayleigh showed, at low altitudes the rays of ditferent colours coming from the 
same star and reaching the same point on the Earth are separated by quite 
appreciable distances in the upper atmosphere. ‘Thus, the separation, Ay, 
between two rays of different colours at a height z above the Earth is given by 
(Rayleigh (1), eq. (10)) 
sin@ — 


\, Ap(2’)d2’ = Apz, (1) 


cos" @ 


sin 6 
Nn = 
1 cos? 
where @ is the angle which the direction of observation makes with the vertical 
and Ap is a mean of the difference in refractive index, Ap(z’) (at height z’) in 
the two colours. Over the visual spectral range Au ~7 x 10°° and 


sin @ ‘ 
An =0'7 35? z in km). 
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From this formula Rayleigh computed that at a height of 8km, Ay = 183 cm for 
4=80°, while it is 45cm for @=70°; and separations of this amount are sufficient 
to account for the movements of the bands in the spectra of stars observed near 
the horizon as described by Montigny (2) and especially by Respighi (3). Indeed, 
Rayleigh was able to give satisfactory explanations for all the observations known 
at that time in terms of the separation of the rays due to chromatic dispersion in 
the atmosphere. On the other hand, for 6<40° the largest separations to be 
expected over heights of the order 10km is less than 5cm; it is evident that we 
must trace the phenomena observed at these higher altitudes to a different cause. 
And, in fact, it has been generally agreed (cf. Anderson (4), Strémgren (§) and 
Hansson, Kristenson, Nettelblad and Reiz (6)) that at these higher altitudes 
we must attribute seeing and scintillation to the corrugation of an incident plane 
wave-front by its passage through a “disturbed” region in the atmosphere. 
Expecting that the disturbed region (which is probably to be identified with the 
inversion layer) will be in a state of turbulence, one supposes that the wave- 
front in its passage through it will be subject to an irregular and variable refractive 
index and become, in consequence, corrugated. It is precisely in terms of such 
a corrugated wave-front that one has sought to interpret the phenomenon of 
astronomical seeing. It is to be noted that the variable refractive index has been 
introduced only to account for the corrugated wave-front. While the manner 
of deducing the properties of seeing and scintillation from the assumed existence 
of a corrugated wave-front has not always been rigorous or quantitative, it, 
nevertheless, seems to the writer that Little’s (7) recent description of the current 
ideas under the heading of a “‘ refraction theory”’ which postulates a linear array 
of cylindrical lenses at a definite height is an unfair representation. In saying 
this, it is not the intention of the writer to criticize Little’s application of Booker, 
Ratcliffe and Shinn’s (8) notions of angular power spectrum and Fresnel diffraction 
by an irregular screen. The application is obviously sound and Little’s use of 
the terminology of Booker, Ratcliffe and Shinn represents a distinct advance in 
the subject. But where Little’s discussion is lacking is also the place where the 
discussion of the other writers has been lacking, namely, the manner in which 
the incident plane wave-front gets corrugated in the first place. It is to account 
for the corrugated wave-front (even if only qualitatively) that the earlier writers 
on the subject had invoked the variable refractive index; and Little’s discussion 
of the problem hardly discredits this. Indeed, in the part of the paper where 
Little tries to deduce the variation in the electron concentration in the F, layers 
from the variable intensity of the radio stars, he implicitly recognizes the need 
for a “refraction basis” for the whole theory. 

From the discussion in the preceding paragraph it would appear that we must 
trace to a disturbed layer the initial cause for the corrugation of the incident 
wave-front; and as has been further suggested, the corrugation may be considered 
as resulting from the fluctuating refractive index in the disturbed layer. On the 
assumption that the conditions in the disturbed layer can be approximated by 
those in homogeneous, isotropic turbulence, we shall outline a theory in terms 
of which the statistical features of the emergent corrugated wave-front can be 
described. Once this has been done, we can apply the notions of angular 
power spectrum and Fresnel diffraction as Little has done. 

In large measure, the theory presented in this paper derives from a method 
which Miinch and the writer (9) have recently developed in another connection. 
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2. The correlation functions to describe the fluctuating refractive index in the 
turbulent layer.—Let 5u(r) denote the instantaneous departure of the refractive 
index from the mean at the point r. On account of turbulence du(r) will be 
subject to fluctuations, and to describe this we shall introduce the correlation 
between the fluctuations 5u(r,) and d:(r,) at two points r, and r, in the medium. 
Under conditions of homogeneity and isotropy (which we shall assume) the 
correlation will depend only on the distance, r= |r, —r,], between the two points 
considered. ‘Thus, we shall let 
5u(ty) 54(t2) = 324 M(r), (3) 

where 62 denotes the mean square fluctuation in the refractive index and M(r) 
is some (even) function of r. 

We shall presently see that for following a ray, statistically, through a medium 
of variable refractive index we need the tensor 


0. G 
R= (5) (are) “s 


giving the correlation of a component of grad dy at the point x, with a component 
at x,’. The definition of this tensor does not require the introduction of any 
additional scalar function; for, clearly, 

— ou 


—- uss 





R 


a 


We therefore have 


— {/M’ M' - M’ 
Ry = aT (= — =) EE 5+ ne (6) 
is it r 


where primes denote differentiation with respect to r. From (6) it follows in 


particular that 
Ns Fo 
aoe ae rf L -— On . 
Ox, . Ox, . (0, 0,0; 0, 0, a5’=r) , r 


We shall denote this transverse correlation by 





—— dM 
d2uR(r), where R(r)= — - ——< (8) 


‘The correlation function M(r) defines a ‘‘ micro-scale”’, 75, such that M(r) 
becomes negligible for r>r,. When we wish to illustrate the general formulae 
by a particular case, we shall assume that 

M(r) =exp (—r*/7r,°). (9) 

3. The equations governing a ray traversing a medium of variable refractive 
index.—The general equations governing a ray traversing a medium of variable 
refractive index can be written in the form (cf. Synge (10), p. 101) 

d px; Tere. . 
: = (7 + XQ? + x,7)!? — 10 
fe \ 1 2 3 ) Ox , ( ) 


du \(x,? + X97 + x5”)"” a, 





where 
- Cl 
A=] 


and u is a single-valued parameter along the ray. 


(11) 
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For the problem on hand, we may clearly regard the variations in » from its 
mean value as small, so that the departures from linearity of the rays traversing 
the medium may also be considered small. Under these conditions the general 
equations governing the rays can be greatly simplified. ‘Thus, letting u =s where 
s is the distance measured along the undeflected linear trajectory and denoting 
the displacements (considered small) along two fixed directions, x, and x, normal 
to s, by é(s) and »(s), and ignoring all quantities of order higher than the first in 
equation (10), we find that 


(25) ma M7-(25 ) 
ds? ax, a on ae Ox, ! ay (12) 


It may be readily verified that equations (12) represent the proper generalizations 
of the equation which Rayleigh has made the basis of his discussion ((1), eq. (7)). 
One of the principal quantities in which we shall be interested is the angular 
deflection of the ray from the linear trajectory. Let ys, and ys, denote the angles 
which the ray makes with the direction s in planes containing x, and x, 
respectively. ‘Then in the approximation in which equations (12) are valid, 


i a _% , 
_— , 4 and tbo = (13) 


ds” 
The equations governing ys, and #, are therefore 


dp, _ oO : df, (0. 
y's =n (= by) and ry = (5 by) . (14) 


4. The mean square deflection and the mean square displacement of a ray after 
traversing a known distance in the medium.— Using equations (12) and (14) we can 
readily express the mean square deflection and the mean square displacement of 
a ray, after it has traversed a known distance in the medium, in terms of M(r). 
Thus considering the deflection, y/,(s) after the ray has traversed a distance s, we 
have (cf. eq. (14)) 


yp 1 


"3s ra} : 2 
b,(s) = ie by) ds). (15) 


From this equation it follows that 


(Ss) =O. 


On the other hand, 


— ‘3 8 0 0 
ws 2 = — 5 L =——— ( So. 


‘The averaged quantity under the integral sign ts clearly the transverse correlation 
(7). We can therefore write 





h?(s) =32u | | R(|s;—52|) ds, dss. (18) 
“0-0 


By a known technique (cf. Chandrasekhar and Miinch (g), eq. (18)) the double 
integral in (18) can be reduced to the single integral 


s,*(s) = 28%u i R(r)(s—r) dr. (19) 


It is, of course, evident that f(s) will be given by the same formula. 
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In order to obtain the mean square displacement, £*(s), we must integrate 
equation (15) once more. Thus 


‘eo fh a 

£(s) - 5 

&(s) Jods |, dt, ( au), : (20) 
From this equation it follows that 


&(s) =o. (21) 
But 


= —ae £0 78 rai (8s 
&*(s) Fn | | ds ds, | J, dt, dt, R(|t, —ts|). (22) 


The quadruple integral representing £(s) can also be reduced to a single integral 
by reductions analogous to those by which a similar quadruple integral was 
reduced by Chandrasekhar and Miinch ((g), eq. (60)). We find: 


Bs) = 25% |” Rls) — br(st—)] dr. (23) 


Since there will be no correlation between the deflections y,(s) and _y,(s) in 
the two directions x, and x,, the mean square of the total deflection, /*(s), will 
be given by 

f(s) = yh,7(8) + Yoo"(8) = 2yh,"(5). (24) 


Similarly, the mean square of the radial displacement, p, will be given by 


p(s) = 2s) +7%(s) = 225). (25) 
Finally, substituting for R(r) in terms of M(r) (eq. (8)) in equations (19) and (23), 
we obtain 
* M'(r) 


f(s) - 4% | (s—r)dr (26) 


0 
and 


p(s) = 45% | MO (93-18) — Jnl? 1?) |dr, (27) 

As we have already pointed out, the function M(r) introduces a micro-scale 
ry such that M(r) becomes negligible for r>r,. Now in the applications, we 
shall be interested, principally, only in values of s which are large compared to 
ry: thus, while s will be measured in units of 100 metres, 7g is expected to be of 
the order of t1ocm. Under these circumstances we may approximate equations 
(26) and (27) to a high degree of accuracy by 


— (* M'(r) 


L?(s) =— 40%us | dr (28) 


0 


% ~o dr. (29) 


p(s) = —$ 5%" | 
0 


M(r) u(*) (30) 
Yo 


we can express equations (28) and (29) alternatively in the forms 


(s) = 40.8%. (=) (31) 


0 


Writing 
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and 


P= 5% (2), (32) 


Basi 3 aS : (»- "), (33) 


M(x)=exp(—x?) and a=y7. (34) 
Using the value «= 4/7, measuring s in 10cm and 7, in 10cm, we find that 
equations (31) and (32) take the numerical forms 


where 


For the case (9), 


V/ p(s) = 1-74 x 107 (=)" 4/5*u seconds of arc 
0 
and 
3\12  —— 
p*(s) = 4°86 x 10° (=) Vo*ecm. (36) 

We shall return to the application of these formulae in Section 7. 

5. The correlation in the deflections experienced by two parallel rays.—In 
Section 4 we obtained an expression (eq. (19)) for the mean square deflection, 
y,°(s), in a direction normal to the ray. We shall now consider the correlation 
in the deflections experienced by two parallel rays incident on the atmosphere 
and separated by a distance X. ‘There is no loss of generality in supposing 
that this separation, X, is in the direction of x,. 

Let s, and s, measure the distances along the two (undeflected) linear 
trajectories. ‘Then, clearly, 








Halos ans = J" J" (3) (8H), ded 


| It may be noted here that the validity of this formula requires that +/£*(s) <7; 
we shall see that this is the case in the applications we shall consider.] Since 
s, and s, are measured along two lines parallel to the x,-axis and separated by a 
distance X, the value of the correlation which occurs under the integral sign in 
(37) is not given by (7); instead, it is given by (cf. eq. (6)) 


M” M’ M) = 
— Fy (> =" or) xt+ 5 =8*Q(r) (say), (38) 


where, now, 


r=[X2+(s,—s,)?]"2. (39) 
Hence, 





y(0; SX; 5) = 25% | ds, i. ds,Q(r). (40) 


Using r as the variable of integration in place of s,, we can rewrite equation (40) 
in the form 


-4/[X*+(8—8,)*) 


W,(0; s\y(X; s)= 28% ‘ ds, \. dr") (aH = 2° (41) 


Letting y=s—s, and inverting the order of the we we find that the 
integration over y can be performed and we are left with 


Vv (X* +8") rdr 


Jil; SW(X; 8) = 25% I. Q(Is— (2 -X*)""] a (42) 
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For s>79, equation (42) reduces to 





Yx(0; s)by(X; s)=25%ps [ Panes dr 


The corresponding correlation coefficient is given by (cf. eq. (28)) 


oe ee B cher rata oss 


p(s) ts a ar 
“0 
For the case (9), formula (44) reduces to 
Ya(0; Sal; in (: ~2 =) exp (— X?/r,?) ; 
f,7(s) ‘0 


this is the same as the longitudinal correlation, 


7) r) M’ d {M'’ 
jai ‘aipake ae — — 92 — 3 6 
(= “(seo , ee r, 0, 0) i [- r vai (+ r )] (4) 


for M given by equation (9). 
6. The retardation in phase on traversing the turbulent layer.—The change 
in phase as the ray traverses the medium will be given by 27L(s)/A, where 

















L(s) = | 5,(51) as (47) 


and A is the wave-length of light considered. ‘lhe mean square change in phase 
will, therefore, be determined by 


L*(s) = i, |" 5u(s,) du(s.) ds, ds, Po [" ; M(|s,—5S2|) ds, dsy. (48) 
0/0 Jo/0 


On reducing this integral over M in the usual fashion, we have (cf. eq. (19)) 
L*(s) 287 u {" M(r)(s —r) dr; (49) 
“0 


and for s>7», we have 
«o 


L*(s) = 25% | M(r) dr. (50) 


0 
For the case (9), the numerical form of equation (50) is 
/ L*(s) = 4°21 x 10%Xsry)! 24/52 cm, (51) 


if sis measured in 10cm and7,in tocm. Also, in this case (cf. eqs. (31) and (34)) 
a :. = 
VPs) = — V Ls); (52) 
ti) 


i.e. the root mean square change in path length is therefore one-half of the root 
mean square angular deflection times the micro-scale. ‘This is a physically 
understandable result. 

The correlation in the changes in phase experienced by two parallel rays 
incident on the turbulent layer and separated by a distance X can also be found. 
Thus, it will be determined by 


L(o . s)L(X ; s)= | "| Bale ) dpu( Sy) ds, ds., (53) 
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where s, and s, are measured along the two (undeflected) rays. Therefore 





L(o; s)L(X; s) BY M(r) ds, dso, (54) 
where pide 

r =[X? + (s.—s,)?}'2. (55) 
Reducing the integral on the right-hand side of equation (53) in the manner 
the integral in equation (40) was reduced in Section 5 (eqs. (40) to (43)), we find 
—— pV (+8) 
L(o; s)L(X; s) = 25% | M(r)[s — (v2 — X2)12] 





r 
(2 — xia dr. (50) 


For s>r,, this tends to 





=- (°% rM(r) 


Lio; s)\L(X; s) 287us . (2 — X22 dr; 


and tor the case (9), this gives (cf. eq. (50)) 
L(o; s)\L(X; s) 
L*(s) 

‘Thus, in this case, the random changes in phase experienced by the different 

rays in passing through the medium directly reflect the variable refractive index 
in the medium. 

7. Numerical illustrations.—For the sake of definiteness we shall suppose 

that the disturbed region in the atmosphere responsible for astronomical seeing 

is a horizontal turbulent layer of thickness d. For a star observed at an angle 0 


from the zenith, the length of the path traversed by a ray in its passage through 
the layer is 








= exp (—X?/r,?) = M(X). 3 (5) 


s=dsec8. (59) 
Inserting this value of s in equations (35), (36) and (51), we have 
d 


Yo 


1/2 
) (sec @)'? seconds of arc, (60) 


_ — 3\ 1/2 
Vp? = 4°86 x 10° (5) (sec 4)? cm, on 
- 


0 
and V L? = 4°21 x 107 1/82u (dry)"? (sec @)'? cm, (62) 
where it may be recalled that s is measured in 10cm and 7p in lo cm. 
Now according to the measures of Strémgren (§) and Hansson and others (6), 
the angular radius of the “ tremor disk’’ js about 0-75 seconds of are under conditions 
of average seeing and near the zenith. Inserting this value for +/¥ in (60), 


we find: 
i ae 
(-) \/ 04 & 4°3 x IOS. (63) 


"oy 
On the other hand, a number of independent lines of evidence lead one to suppose 
that 


ry~ 10cm, 1.e. I in the units adopted. (64) 
Hence, combining (62), (63) and (64), we infer that 
V/ L? = 4:2 x 10? x 4:3 x 10° 8 =1°8 x 10° = 1 800A. (65) 


A change in path length of this order is entirely sufficient to account for scintillation 
in colour. 
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Again, for rg =I (i.e. 10cm) equation (63) gives 
D24/52n ~ 4 x 10-8, (66) 


If we assume that the thickness of the turbulent layer is 100 metres (i.e. d=1 in 
the units adopted) then we should conclude that 


Ven 4x 10 e, (67) 


‘This corresponds to a root mean square fluctuation in density of amount 10 


times the mean density; since a variation in density of this amount may be 
expected to occur over a distance of the order of r, (10cm), we are postulating 
gradients of density of the order of 10->cm™ in units of the mean density; the 
assumption that density gradients of this order of magnitude occur in the 
turbulent layer does not seem unreasonable. Finally, we may note that for 
d=rg=1 and 1/4 > 4x 10-8, equation (61) leads to 4 pi~e2 x1o-%cm; and 
this is indeed small compared to ry (~ 10 cm). 

From the preceding discussion it would appear that we can satisfactorily 
account for normal astronomical seeing by postulating a turbulent layer of a 
thickness of the order of 100 metres, a micro-scale of turbulence of the order of 
10cm, and a root mean square fluctuation in refractive index of_the order of 
410%. And it would follow that to explain bad seeing when 1/7 can become 
as large as 10 seconds of arc (or even larger, occasionally), we should suppose that 
ry is smaller and/or d is larger. And finally, it would also appear that in astro- 
nomical seeing we have a means of exploring the spectrum of turbulence in the 
atmosphere at heights of the order of 3 or 4km. 


Yerkes Observatory, 
Williams Bay, 
Wisconsin, U.S.A. : 
1952 February 19. 


References 
(x) Lord Rayleigh, Phil. Mag. [5] 36, 129, 1893; also Scientific Papers, Vol. 4, p. 60. 
(2) C. M. V. Montigny, Mem. de l’ Acad. de Bruxelles, 28, 1856. 
(3) L. Respighi, Bull. Assoc. Francaise pour l Avancement des Sciences, 1, 148, 187 
(4) J. A. Anderson, Your. Optical Soc. America, 25, 152, 1935. 
(5) B. Strémgren, Matematisk Tidsskrift, B, Festskrift Norlund, Part I, p. 15, 1945. 
(6) N. Hansson, H. Kristenson, F. Nettelblad and A. Reiz, Ann. d'Astrophysique, 
13, 275, 1950. 
(7) C. G. Little, M.N., 111, 289, 1951. 
(8) H. G. Booker, J. A. Ratcliffe and D. H. Shinn, Phil. Trans. Roy. Soc. A, 242, 579, 
1950. 
(g) S. Chandrasekhar and G. Minch, Ap. F., 115, 103, 1952. 
(10) J. L. Synge, Geometrical Optics, Cambridge Mathematical 'T'ract No. 37, Cambridge, 
1937: 





MODELS FOR UNMIXED STARS 
A. E. Roy 


(Communicated by the Director, University Observatory, Glasgow) 


(Received 1952 April 4) 


Summary 
A set of point-convective models of a star with a range of discontinuities 
of molecular weight between core and envelope has been derived using a 
modification of the Bondis’ method. ‘These models are used to examine the 
evolutionary spread of unmixed stars to the right of the H—R diagram. It is 
found that the radius of an unmixed star could increase by a factor of ten 
owing to the exhaustion of hydrogen during evolution. 





1. /ntroduction.—I\t has been pointed out by Hoyle and Lyttleton (1, 2) 
that a discontinuity in molecular weight in a star will have an effect on the radius 
of the star. Models have been calculated by them in which the discontinuity 
is formed by the accretion of hydrogen by the star and it has been shown that 
giant stars may be formed in this way. In this case the discontinuity lies in 
the outer part of the star. Again, in the models calculated by Li Hen and 
Schwarzschild (3) a difference in chemical composition also occurs in the envelope, 
while in the set of papers by C. M. and H. Bondi (4, 5, 6) all possible models for 
red giant stars are discussed systematically. While we are not concerned in 
this paper with the problem of the red giants, it may be noted that the set of 
models calculated in the present paper are limiting cases of the above where the 
discontinuity occurs at the convective core interface, due to the conversion of 
hydrogen into helium in the core. 

‘he existence of such a discontinuity will depend upon the star’s rotational 
velocity. In recent years the importance of rotation in stellar structure has been 
pointed out by various authors. Stars can be mixed or unmixed, depending 
on their rotational velocities. In a mixed star, no discontinuity of molecular 
weight forms at the core interface, since the stirring allows hydrogen from all 
parts of such a star to be consumed and spreads the resultant helium throughout 
the star. But in general, stars of spectral type later than F5 are likely to be 
unmixed (7). It therefore seemed of interest to calculate a set of models of 
unmixed stars and contrast their evolutionary tracks on the H—R diagram with 
those of mixed stars. It is to be noted that stellar models with discontinuities 
in molecular weight due to accretion are not examined in this paper. 

In this paper the homology-variables introduced by C. M. and H. Bondi (8) 
are used, viz., 

5= Frye a \ dinp 
ONC On ¥"om’ “~*~ tar 
We introduce as an independent variable u, where 
u= —In(M/ M,). (1.02) 


(1.01) 








No. 5, 1952 Models for unmixed stars 485 


One advantage of this is that the physically important quantity M/M, is 
obtained automatically during the integration of the model; we also obtain 
differential equations better suited for treatment by numerical integration than 
the equations (4) and (6) of the Bondis’ paper (hereinafter referred to as Paper 1). 

The equations of structure in radiative regions are 

dP/dr = —GMp/r’, (1.03) 

dM/dr=4nr*p, (1.04) 

dT/dr = — Bp®/r*T”. (1.05) 
Here o, v take constant values in two important cases when radiation pressure 
can be neglected: 

(a) radiative envelopes of uniform composition, with Kramers’ opacity law 
and guillotine factor approximated to by a power law; 

(6) radiative envelopes of uniform composition, with electron scattering 
opacity. 

For (a), o=2, v=6°5, neglecting the guillotine factor; for (b), o=1, v=3. 
In the absence of radiation pressure 
P=pAT|p. (1.06) 
Also 


dinT 
= ye SD (1.07) 


N 


We have 
(—rQ/N)(dN/dr)=S+oa-—(o+v)N. (1.08) 
Now 
— S/Or= —4nr*p/M= —d|n M/dr. 
But 
du=—d\|n M. 
Thus — Qr/S =dr/du and equation (1.08) becomes 
S aN 
N du 
; _SdN S+a—(v+o)N " 
Now we have Nas” 1+2S-40’ from Paper I, and so 
dS/du=1+2S—40Q. 


=S+a—(a+v)N. 





Also, by Paper I, 
SdQ_ S+N-Q 
OdS  1+2S—4Q’ 
so we obtain 
$49 
O du 
Rewriting (1.09), (1.10), (1.11), 
dN/du=(N/S)\(S+B-«N), 
dS/du=1+2S—4Q, 


dQ/du=(Q/S)\(S’+N-Q), 


=S+N-Q. 


where 8 =o and x =1' +0. 
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With w as independent variable the series expansions of the solutions in the 
neighbourhood of the origin are 


Pee. 


4fh?aué 





V=- 4+ 


u— 
% oa(1+f) 


42°BlaB(38 —2) 4 


(1+B\a+Pya+ B+ ap) 


+ 6% 


11f? 


- «7 |u6? 





(1+ P)(a4 


4 
. i 
“a+ fp 


403(38 — 


Q 


pies + 2PB)\(4+P4 
2(38 —«)ué? 





x)[88? + 3a8 —o 


Bla B+ 2B) 


27 |ub 


aB)(a+2P +a) 





’ 2 
Te + 





t+ 3P) 
(38 —a)[8A? + 348 — a*|0% 





— IQa* 38)64 


B?(a + B)*(a + 28)(x 
B = % pe “’ 

-0+(5) 5 a+B : +(j (a + B)*(a + 2B) 

I (~\8(38 - mat? 

? 5 (5) 


«)[ 32684 + 304 — 418?x? + 211485 
Ss [: - 


(~ +P) (%+28)(a + 38) 
Pe Ae ué 4 
a(1 +f) 
where @ = Au’ and A is a constant. 
The terms 6 in the expression for QO and ué* in S, which have been omitted, 
are small compared with the other terms in the range of wu used. On substituting 
the appropriate numerical values for «, 8, we obtain 








(a) Kramers’ opacity model: 
N =0°235 294 + u(0-078 431 —Or1 
S : — 3°238 100 + 2°752 4067 + 0-825 9447 +u+...), 
O=6 I°O11 9050 — 0-352 2467 — 0-011 720° + 0-843 14u+...). 


(6) Electron-scattering model: 
N =0:25 + u(0-125 — 0°177 780 — 0-177 7867 +...), 
S =u(1 — 3°20 + 2°133 3307+ 0-975 25f +u+...), 
O = 0(1 — 0-80 — 0-426 676? — 0-268 196% +... +0°875 00u+...). 


2. The boundary conditions at the core surface.—Consider a discontinuity in 
molecular weight at tue interface between convective core and radiative envelope. 
The following boundary conditions hold at the interface. Suffixes 7 and o 
refer to values of a quantity on the core side and envelope side of the interface 
Since P, T are continuous, 


57 000 — 0°:056 5507 + 


respectively. 

Poli=HelHe, since P=p&T |p, 

pw, and yw, being the values of the molecular weights in envelope and core 
respectively. 

Also, by definition, S and »Q are continuous. 

transfer then determines the discontinuity of N at the interface, giving 


N, e (e2)’ - (H:) PK oe 
Pi Ke Koc ‘ 


The equation of radiative 


Pe Koe 


He oc 


(2.01) 


N, 


where Kk ,, Ko. are the opacity constants in envelope and core respectively. 
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In this paper the Kramers opacity model is discussed. We then have « =8°5, 

8 =2 and (2.01) becomes 

Koe He 


Koc Me 


N,=0°4 


0 
Introducing X, Y, U, where U=1—X-— Y and X, Y are the masses of hydrogen 
and helium respectively in the star’s composition per gram (9), 
p=(2X+ 3V+$U)4. (2.02) 
In the model dealt with here, U is taken to be the same in core and envelope and 
kg=h(1+X)(1—X— Y), 
where & is a constant taken to be the same in core and envelope. ‘Thus 
I+X, 2X,+2Y,+4U 
I+X, 2X,+2Y,.+4U 





N,=0°4 (2.03) 
In order to examine the widest range of discontinuities in molecular weight 
possible under the above circumstances, a set of models was chosen consisting 
of 99 per cent hydrogen, 0 per cent helium and 1 per cent other elements in the 
envelope with cores ranging from 99 per cent hydrogen, 0 per cent helium and 
I per cent other elements to 0 per cent hydrogen, 99 per cent helium and 1 per 
cent other elements. ‘Then 


b, =(2X,+ }U)* =0°503 78, 
f= (2X, +0°75Y,+0°005) 1? =(0°747 5 +1°25X,) 3, 


and (2.03) becomes 





N, =0"401 0 (CAS), (2.04) 


The resulting changes in N,, », and p,/u,, as X, changes from 0-99 to 0, are given 
in Table 1. A graph of y,/u, against V,, taken from ‘Table I, helped to determine 
where to cut off the set of solutions required for various values of constant A. 


TABLe I 
dhe No fre 


0°400 00 0°50378 
"70 0°382 71 0°616 33 
*50 0°36691 0°72860 
*30 0°34625 0°890 87 
‘10 0°31807 1°14613 
‘00 0°29975 1°33779 


A graph was also drawn of QO against S for N =o0-4 by using a table constructed 
from the relations between the Bondi variables and the Emden variables x, y, 
tabulated (10) for m=1-5. Hereinafter these graphs are referred to as graphs | 
and II respectively. 

Using the modified form of numerical integration outlined in the Appendix, 
eight integrations were performed of equations (1.12), (1.13) and (1.14) for 
eight values of A, each value corresponding to a value of p,/p,. 





488 A. E. Roy Vol. 112 


3. The p-variables in terms of the h-variables.—The p-variables are given in 
terms of the h-variables using equations (12) of Paper I, with u replacing S as the 
independent variable. We then have 


(3-01) 
(3-02) 
(3-03) 


(3-04) 


The luminosity and radius of the resulting star can now be obtained in terms of 
its mass by means of the following relations. 


sri _ ("30-8 o5)* 
(a) In 3M, Je 5 du. (3.05) 





(6) Adopting the law of energy generation 
e=EpT", (3.06) 


where E is a constant throughout one star and from star to star we have 


L= | “4nr°p edr 
0 


am, (ATA(n+1))3? (*. 
= 4nE pT | , [ ‘strata, 


47.Gp, 


_, { 474 Gee | 
~ |AT(n+1) 
is the radial coordinate in Emden units. Now 
d*x 2dx 


dz* _ zdz 


hence x=I—}27+.... 


0 


where x= 7/T, and 


+x"=0, 


On account of the strong temperature dependence of the rate of energy 
production we can take 
x=1—}2? 
in the region contributing most of this energy. Hence, on using x as independent 
variable, we can write 


2 1 
| "x2" 4" 3%dz = 34 o| x” (1 — x)" dx = 31/6 B(m+ 1, 3), 
0 0 


where m=2n+7. 
With » =17 and n = 3 in the core, m= 20. 
Hence equation (3.07) becomes 
, 3/2 
L=0:2661 7Ep? T," eon Cy 
47.Gpe 
which is the required relation. 


* After equation (13) of Paper I. 
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(c) We have 
Thus 
47 Par 


In M2 


(3.09) 
Also, by equation (14) of Paper I, 


’ 


P+. hae 
In —— InO-—u-+ | ~=du+ — du. 
cVil4s / 0+ /u S 


Equation (3.09), multiplied by N,, plus In A minus equation (3.10), gives 


B, “© N —N AS’ |. “O 
In(47) =} oe + In Oo H(4N,—1) | du (2N,—1)u, 


2 “ 


4n P.r,* Ns ; x 
b, (e a) and B, 


\s u>o, 4S*8/O-+1, and on noting that V,= f/x, we obtain the required relation 


PJ-Play 1-40)a r°N—N 
: >=IndA+ | bt, (3.11) 
Jo S 


where 





J 
In 4 4r)Piz Gib 9, M2 28 af 


(d) ‘The equation of radiative transfer 
-_ > - B 
dl 3 ob _?p 
dr 167 ac r* Te? 
can be treated as in Paper I to give, in the notation of the present paper, 


> " J p \ 
yx 281 > Ko Lf 2&\* (47) B SS pt 4B 
167 ac \yw,G NO* 


S8 t, 
NO® > B 


/ 


lim 
So 
hence at the surface, (So), (3.12) becomes 


yxei% 3 Koel ( # 
B° 428° plBac\uG 
which is the fourth relation. 
Let 


(7° 30-S } 
eXP 4 | <i dus=e (A). 


‘The method of evaluating ey and e, is given in the Appendix. Then, using 


\- 


relations (3.05), (3.08), (3.11), (3.13), we obtain 


D Vi (x 
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F Gu, : . 
ppiycs Air) 7390 yee 


+B 148 
(-) A* yp4b- a nisi Ais... HG" 
| i 4 


> 
3) 


where 
and 


where 


a 


a= —-3, 
—(e~ Pp) 
2(a—pf 
—a(2n + 3) 
2(a— PB) . 
I 
a—4B—n-3° 
Denoting the continuous solution, (j,/, = 1), which is the normal Cowling model, 
by suffix a, the luminosity L, and the radius R, of a star without discontinuity at 
the core surface are given in terms of its mass M by } 





a; = 
a= 


as = 


(a—2B—n—1)(a—48—n—3) 
R,=®2, M“ 1 vies 


L _M- (28 +1)—a—2B—3)/(a—48 n—3) 


a 


Then for a general wena star of the same mass with a discontinuity at the core 
surface, the luminosity L and radius R are given by 


© _ 0 ™ 
R=—R,=9R,; L re OL, (3.19) 


®, 
where for stars with a given law of opacity and energy generation, ® and © 
are known functions of «,, uw, only, given by 


a if be n A ay ey a fe.\% ts 
wr Ue Wet et 
(ey eel ay 


where 





In practice the luminosity L can be derived, once R has been found, from 
following relations. 
(1) For Kramers’ opacity 
LR"? =(A/A,)** LR,"". 
(2) For electron opacity 
L=(A/A) L.. (3.21) 
Calculating values of the radius, luminosity and central temperature for each value 
of A in units of the corresponding v ralues for the continuous C owling model by 
using the expressions (3.15), (3.16), (3-17), (3-18), (3.19) and the expression 


Z. LM, R e ox /(a—f) e. Bi(a—B) / 4 \a(a-B) 
Eee ow) 


the changes that take place in an unmixed star were plotted in Fig. 1 
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It is of interest to calculate the rate at which stars of different luminosities use 
up the hydrogen in their cores. In ‘Table I] below is shown the approximate 
time, in units of 10’ years, taken by stars on the present model, when, = 0-503 78, 
to reach the value of j,/, at the core interface tabulated in the top line of the table. 

R T T T T Y T me 
Rc owling Lcowling 

2 
TeCowling 





Luminosity 











i L 1 i 1 
1-25 1-50 1-75 2:00 2-25 2:50 2:75 
he 
Fic. 1.—Values of the radius, luminosity and central temperature of a star with Kramers’ opacity, 
in terms of the Cowling continuous model, at various values of the discontinuity in molecular weight. 





TABLE II 





j 


i 
| Mbol 





| 


"92 II 000 16 000 
"04 600 930 
‘O02 250 430 
“41 35 131 
‘O52 13 21 : 
“712 1°9 2°9 Oo ' 4°4 

34 0'042 0°065 . 0'098 


























4. The evolutionary tracks in the H-R diagram.—-\n order to obtain the 
evolutionary tracks of the unmixed stars in the H—R diagram, values of m),,, and 
In 7, were calculated, where m,,, and 7’, are the absolute bolometric magnitude 
and effective temperature defined in the usual manner by the equations 


M,,., = constant — 2°5 logy, L, (4.01) 


¥ L \wa 
(= =) (4-02) 
35* 
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1 ga, (CHe\"* EE iy, ~ 
AN R (477)!!23"s¢ yee, ad 6 (3-17) 


J 


where 


and : 


42+8 71 +Pac (u,G\* +! ; 
Si gn RP ’ (3.18) 


“) 


where 





as = 


Denoting the continuous solution, (j,/u, = 1), which is the normal Cowling model, 
by suffix a, the luminosity L, and the radius R, of a star without discontinuity at 
the core surface are given in terms of its mass M by 
R 
L 


Then for a general unmixed star of the same mass with a discontinuity at the core 
surface, the luminosity L and radius R are given by 
@ — 1) _ 
R D. R, @R.; &L 57 ©L,; (3-19) 
a 


a 


(a—2B—n—1)/(a—4B—n—3) 
®,M n—3) 


a 


(-) M! (28 +1)—a—2B—3]/(a—48——3) 


a 


where for stars with a given law of opacity and energy generation, ® and © 
are known functions of p,, pw, only, given by 


Snffey (2) (2"(s)\"\* 
1563 | be, A, CNa esa 

SICKO KORA 
\ oe CENa Ca 4 Se 


) j 
3a] o + 26(1 4 n)| 


where 





ae 
* 2(a— fp) 


In practice the luminosity / can be derived, once R has been found, from the 
following relations. 
(1) For Kramers’ opacity 
LR‘? =(4/A.)** LR.™. (3.20) 
(2) For electron opacity 
L=(A/A,) L,. (3-21) 
Calculating values of the radius, luminosity and central temperature for each value 
of A in units of the corresponding values for the continuous Cowling model by 
using the expressions (3.15), (3.16), (3.17), (3-18), (3.19) and the expression 


: ue\ Rf ex \2 ( e, \P' )/ A \=Ka-B) 
See 


the changes that take place in an unmixed star were plotted in Fig. I. 
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It is of interest to calculate the rate at which stars of different luminosities use 
up the hydrogen in their cores. In ‘lable IL below is shown the approximate 
time, in units of 10’ years, taken by stars on the present model, when, =0°503 78, 
to reach the value of y,, 1, at the core interface tabulated in the top line of the table. 

L 
Lcowling 
T. 
TeCowling 





R T “T T T T T 


Re owling 


Luminosity 





wy 


perature 
Cent 


ral Tem 





i 1 ) 1 rt j 
1-25 1-50 1:75 2:00 2-25 2:50 2:75 4 
Me 








Fic. 1.—Values of the radius, luminosity and central temperature of a star with Kramers’ opacity, 
in terms of the Cowling continuous model, at various values of the discontinuity in molecular weight. 


TABLE I] 








II 000 16 000 

‘ 600 930 
"O02 280 430 
“41 85 131 
‘oS2 13 21 
“712 19 2°9 3° : 4°4 
$4 0°042 0'065 >*o8 . 0098 























4. The evolutionary tracks in the H-R diagram.—\n order to obtain the 


evolutionary tracks of the unmixed stars in the H—R diagram, values of m,,, and 
In 7, were calculated, where m,,, and 7, are the absolute bolometric magnitude 
and effective temperature defined in the usual manner by the equations 


m,,., = constant — 2°5 log, L, (4.01) 
(4.02) 
5* 
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Kuiper’s bolometric corrections (11, Tables 9, 13) can be used to reduce absolute 
bolometric magnitudes to absolute photovisual magnitudes. Fig. 2 gives the 
tracks. 

4:3 , 3-9 3:7 3°5 log Tp 





' Lf T 




















l L. i | a | i j 
BO B5 AO AS FO GO KO KS5*+ M2 








Fic. 2.—Evolutionary tracks of mixed and unmixed stars. 


In order to contrast these tracks of unmixed stars with the tracks of mixed 
stars, a set of mixed stars was placed in the H—R diagram. Once again Kramers’ 
opacity was used. Also »,=s,= where p varies from 0-50 to 1-34 as the star’s 
composition changes from X =0-99, Y=o to X=o0, Y=o-9g. At the interface 
between radiative envelope and convective core P, T are continuous and so O, S, V 
are continuous. ‘Thus the required model is the continuous Cowling model. 
Using equations (3.14), (3-15), (3.16), (3-17), (3.18), we obtain after some 
reduction 

R =(14/p.)"” R,, (4.03) 
ne 
L = (1/14) - L,, (4.04) 
12 — (a, fy, \7°5 1/2 Pi 
LR* =(u Ha) * Egle ’ (4.05) 


where R,, L, are the radius and luminosity when the mean molecular weight is 
p, and R, L are the appropriate values when p is some other value. Using the set 
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of values of ./u, that we had for y,/u, in the unmixed case, the set of mixed stars 
obtained was plotted in Fig. 2. In addition the main sequence and giant sequence, 
adapted from Otto Struve (12) is shown. It is to be noted that Struve’s curve for 
the main sequence is highly selective and very narrow. 

Line AB was fixed by adjusting the value of the opacity constant so that the 
line of solar mass in the unmixed case passed through the Sun’s position in the 
H-R diagram. It is of interest to note that, on the present model, if the Sun was 
originally a star of mean molecular weight one-half, its present age would be of 
order 5 x 10% years. Discrepancies between the line AB and the main sequence 
must be expected, since it is possible that changes take place in the laws of opacity 
and energy generation as one moves up the main sequence. For example, 
electron opacity becomes more important in the stars of higher mass and a change 
to proton-proton energy generation has been shown to be probable by H. Bondi 
(13) for stars of one-sixth the luminosity of the Sun. But in general there ts 
satisfactory agreement, since stars of high luminosity will have had time to spread 
to the right of the line AB, while low-luminosity stars have not had time to deviate 
far from line AB. 

It is seen that the two populations of mixed and unmixed stars do not mingle in 
the H—R diagram. 

APPENDIX 

5. The integration of the equations.—The integration of the set of equations 
(1.12), (1.13), (1.14) is commenced at u=0, that is, at the surface of the star, by 
using the series (1.15), and then continued either numerically or graphically. 
The equations are easily and accurately treated with a hand-operated calculator 
using the following method of numerical integration. 

In the expression 

Vy =Y_y + 2e(1 + 45? — 54+ ...) yp’, (5.01) 
zw is the interval of tabulation and 6? is the double central difference operator, etc. 
(5.01) was used at such an interval that the third term in the bracket was always 
negligible. ‘Io avoid errors due to the necessity of estimating the second term, the 
whole integration was repeated. ‘The intervals in u ranged from 0-000 25 to 
0-128 by doubling or halving the interval at appropriate stages. 

6. The boundary values at the interface.—The eight integrations of equations 
(1.12), (1.13), (1.14) were cut off at the correct values of N, (see Section 2) in the 
following way. 

It would be possible to choose a set of values of the discontinuity in molecular 
weight (for example, the set given in Table 1) and then find those values of the 
parameter 4 that would give solutions satisfying the boundary conditions at the 
interface. However, a less laborious method was adopted, consisting essentially 
of finding the value of j,/, which corresponds to a given A. 

‘The following two equations of condition hold at the interface 


1-985 ) I 

=( ~-~ —0-7475)—., (6.01) 

: G Q, 1°25 

B 0°796 

No= —- 6.02 
eo (1+X,)0,/Q, ey 
‘hese equations are equivalent to the single equation (2.04). ‘lo derive a solution, 
a value of A is chosen and the integration of equations (1.12), (1.13), (1.14) for 
this value of A is continued until the values of N approach 0-4. Forming the 
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quantity Q,,/Q, for various values of N from the integration and graph II, Q,/Q, is 
plotted against N in graph I. It is seen that the correct boundary values will 
occur where the graphs intersect. When this happens the integration is cut off 
and more accurate values of N,, u,/u,=Q,/O,, X, are obtained by choosing that 
value of Q,/Q; which will satisfy equations (6.01) and (6.02) and also give an \, 
equal to the one corresponding to that value of O,/Q, in the integration. 

‘Table III below summarizes the results: the meaning of the last two columns 


is defined in Section 7 following. 


TABLE III 


M, 


[p/ ft, é No u —In M, é os 


m 
I°OI9 O 
11541 
1°427 5 
1°624 1 
1915 6 
2°131 6 
2°363 2 
2°637 9 


‘235 1 0°398 5 "925 8 0°265 18 3°655 0 
"241 5 0°387 9 ‘952 0°283 43 3°987 0 
‘251 3 0°368 2 ‘985 0°320 00 4°578 7 
256 8 0°355 2 998 0°346 97 4°966 8 
263.5  0°337 6 "999 5 0°387 04 5°4929 
268 o 0°325 6 ‘995 6 O°419 24 5°884 7 
‘2725 03137 "983 0 0°454 83 6295 8 
‘278 0 0°*300 6 954 4 0°499 36 6-800 o 


ooo0o00000 


As a test of the accuracy of the method, values for the Cowling model, (,/, = 1), 
were extrapolated from the set of solutions and compared with those obtained by 
J. G, Gardiner (14) and others. 


a - 
R (Emden units) ’ R core Pa Pmean 


M sar Peentral 
Present paper 

J. G. Gardiner 
Cowling 

Hoyle and Lyttleton 
C. M. and H. Bondi 
C. M. Bondi 


"O14 5 ‘169 938 
‘OII 7 "170 04 “146 38 0°027 206* 
"027 *169 ‘144 0°027 


0'146 38 0°027 162 
° 
° 
*967 “E71 0'146 0°025 
° 
° 


anrnxn 


"026 ‘170 *147 0°027 
‘662 "163 ‘180 0'020 


sis 


* Value computed from 2,, X,, %., as given in Table V of Miss Gardiner’s paper. 


7. Evaluation of ey and e,.—Let 


[f° N-N, , | 
exp, =ey=exp;| —>o—du}, 


~ 0 a J 


: ((°30-S , ) 
exp C,. =e, = exp4 | ~~ di 
J0 S J 
\lso, we can write 
in c 
Si TS2> 
; HN ~ J 
where £ | ———— du 
/0O ‘ 
- eg N -N, 
and f, — du, 
Jue (4S 


with similar expressions in ¢. ‘Then, by considering equations (1.12), (1.13), 
(1.14) it can be shown that 


¢ ee Neg ai U, + InN, N, (7.6 I) 
rae }, _——_ au = rs f* 
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“e30-—S 3 \ ; a 2(a - a—B)+3 3. 
du = Ss? 7.02 
4 Jo S = a(n) (a) J 4-48 si 


where suffix c denotes the value of the quantity at the interface on the envelope side. 

The right-hand members of (7.01) and (7.02) being known, these integrals can 
be evaluated. The value of €, found in this way is very accurate, but because of 
the powers appearing on the right-hand side of (7.02) and because the right-hand 





TABLE IV 


Q a be = et 


0°266 61 
0°302 33 
0°337 68 
0°372 67 
©°407 29 
0°441 56 
0°475 86 
O°510 OI 
0°544 00 


(7) 
CY 


—0°642 35 
—0°516 45 
—0°426 22 
—0°358 35 
—0'305 46 
—0'263 15 
-701 I —0°'229 69 
612 7 —0°202 22 
540 2 —9O°179 34 


0°580 71 
0°523 05 
0°476 23 
"437 06 
*404 06 
"375 62 
“351 14 
"329 48 
"310 54 


ieee. te. By, Be B.- 
wnNne OVO CONTI AM 
Or + Se eH He me 
C0CCOnNHHN Z 


o'611 69 
0°679 31 
0°746 68 
0813 81 


‘278 41 
"252 36 
*23° 74 
"212 59 


429 2 —0'143 81 
349 6 —o'118 13 
290 3 —o0-098 70 
2449 —0'083 60 


me Soxiu 


00200 
0000 


o'914 62 
I‘O15 22 
1-115 62 
1‘215 


194 6 —0'066 90 
158 3 —0°054 73 
131 4 —0'045 56 
1107 —0°038 52 


‘190 15 
‘172 00 
“157 02 
"144 45 


4 
7 
me) 
3 


22900 
0000 


"116 
‘097 55 
083 
°073 
"065 55 
‘059 


"516 
"817 33 
‘1177 
“418 
"718 
‘o18 


‘071 46 0°025 06 
"049 QI 0'017 56 
"036 79 —0'013 00 
‘028 29 —O'OIO 02 
"O22 40 —-0°007 95 
o18 17 —0'006 462 


aS) 


‘ 


233 S28 
se 
000000 


3°618 
4°219 
4°819 
5°419 


o12 66 —0'004 508 
009 32 -0°003 228 
007 15 —0°002 557 
‘005 65 0'002 026 


"O49 « 
"O42 
037 


2000 
0000 


6°619 
7°819 
g°o19 
10°219 


003 792 OOO! 359 § 
‘002 719 —0°000 972 7 
‘002 044 0°000 733 § 
‘OOl 593 ©'000 572 3 


20900 
coooo$o 


12°620 

15°020 3 
17°420 4 
19°820 5 


‘OO! O45 —0'000 3764 
‘000 738 —0°000 266 3 
"000 549 —o°000 198 2 
‘000 424 —0°'000 153 3 


ooo0°o 


24°620 ° ‘ —0°000 275 —0*000 099 3 
29°420 : , —0°'000 192 —0°'000 069 3 
oe 
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side is also the difference of two similar numbers, the accuracy of ¢, is not as high 
as that of €,. It is accurate enough, however, for most purposes. 

For example p mean’? cores USing the calculated value of ¢, for the Cowling model, 
is given as 0-027 30. Miss Gardiner’s value is 0-027 2. 

Now in the core, i.e. where S,<S< oo, (N—VN,) is constant; thus, after 
changing the independent variable from u to S and using the values of Q and S for 
N =0-4 obtained from the tables (10), the functions 

és { . dS Se | 2 (30-—S)dS 
N-N, Js S(1+2S—4Q)’ ° s, (1+ 2S —4Q) 


were tabulated for S,>0-5 in Table IV. ‘This table therefore can be used for all 
stars with convective cores in which radiation pressure is neglected. 

Detailed tables of the set of models for unmixed stars have been deposited with 
the Royal Astronomical Society and are available on loan to anyone who desires 
them from the Assistant Secretary. 
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THE DETERMINATION OF THE POSITION OF A RADIO ST? 
F. G. Smith 
(Communicated by M. Ryle, F.R.S.) 
(Received 1952 April 18) 


Summary 


A number of new methods of using spaced aerial interferometers to 
letermine the position of a radio star have been devised, and their accuracies 
are here compared with those of methods which have already been used. 
he new methods are particularly applicable to the accurate determination 

t the positions of a small number of the most intense radio stars. 

The apparent position of a radio star as determined by all such methods 
will differ from the true position because of various refraction effects in the 
Earth’s atmosphere; the magnitude of these effects at various wave-lengths 
is estimated in this paper. 

It is shown that the results may be confused by the simultaneous 
‘eception of radiation from other nearby radio stars, and the resulting erro 
is discussed. 





1. Introduction 

‘lwo main types of interferometer have so far been used for the detection and 
location of radio stars: the ‘ Lloyd’s Mirror’’ interferometer based on the use ot a 
single aerial on a high cliff overlooking the sea (Bolton and Stanley (2)) and the 
‘““Michelson”’ interferometer which uses two aerials spaced laterally in an East- 
West direction (Ryle and Smith (7)). ‘The purpose of this paper is to describe a 
number of new methods of using the spaced aerial interferometer, with particular 
reference to improved methods of determining the declination of a radio star. In 
Section 3 the limiting accuracies of the original and the new methods will be 
derived and compared. 

Consideration is also given to the etfects of refraction in the troposphere and 
ionosphere, and to errors in the observed positions of radio stars due to the 
simultaneous reception of radiation trom other nearby radio stars. 


2. The spaced aerial ‘nterferometer 


(a) The fundamental method of observation._\he power intercepted by a 
spaced aerial interferometer from a radio star varies periodically as the star passes 
through successive maxima and minima of the interference pattern. ‘The spacing 
of the maxima depends on the separation (d) between the aerials, and the total arc 
over which the pattern extends depends on the East-West aperture (5) of the 
individual aerials; the number of cycles of the pattern is thus related tod/b. ‘The 
time at which the central maximum occurs is the time at which the star crosses the 
plane of symmetry perpendicular to the line joining the aerials; this plane is here 
referred to as the axial plane. If the axial plane corresponds exactly with the 
meridian plane, the local sidereal time at which the transit occurs is the apparent 
right ascension (x) of the star. It will later be important to examine the possible 
errors in determining this time. 
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(b) The relation between the observed times of maximum signal and the position of a 
radio star for an East-West interferometer.—Consider first an interferometer aligned 
exactly on a horizontal East-West line. If a signal is received from a distant source 
of hour angle H and declination 6 (Fig. 1) the phase difference 4 in the path to the 
two aerials is given by: 

2nd . . 
db sin H cos 6, 
where d is the separation of the aerials and A is the wave-length. 


P 





Q 


Fic. 1.—The hour angle H and declination 8 of a star. PSQ is the observer's meridian. 


As the Earth rotates the hour angle of the star increases, passing through zero 
at upper culmination, and the received signal will vary periodically as the phase 
difference between the signals received by the two aerials changes. ‘The central 
maximum will be recorded* when H =o0, and other maxima will occur when 


: na 
sin H= — sec 4. 
d 


A circumpolar star will also produce a similar record at lower culmination 12 
hours later, provided that the aerials can be made sufficiently receptive in the 
direction of the star at this time. 

(c) The relationship for a general orientation of the interferometer.— Suppose 
now that the two aerials of the interferometer pair are not exactly on the East-West 
line and are not at the same level. ‘The axial plane will no longer contain the 
meridian; the great circle in which the axial plane cuts the celestial sphere will be 
called the “fictitious meridian’’ (Fig. 2). The relation between the true and 
fictitious meridian is given by the angles here referred to as “ projected azimuth 
error’’ (m) and “ projected level error’’ (m); these are known as Bessel’s angles in 
conventional positional astronomy. ‘They are related to the alignment of the 
interferometer by the equations given in Appendix A. 

* In practice improved accuracy is obtained by the use of the phase-switching system, when 
observations are made of the times when the radiation intercepted by the two aerials is in phase 


quadrature (Ryle (6)). For simplicity, however, the present analysis will assume that the times of 
the maxima are measured. 
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It is seen in Fig. 2 that a star moving round the small circle LM will be observed 
to make its transit at a time later than the time of transit across the meridian, the 
difference being the time taken to traverse the arc I'l”. ‘This time is a function of 
the declination of the star and in some of the methods described in this paper use 
is made of this fact to determine the declination. 





Fic. 2.—Celestial sphere showing fictitious meridian XY with East point at A. 
Projected azimuthal error m= BE. Projected level error n= BA. 


It is shown in Appendix A that a star of right ascension « and declination 6 will 
cross the fictitious meridian at a local sidereal time © given by : 
O=a+H,, 
where H, is given by: 
sin (m— H,)=tan 6 tan n. 


3. Methods of using the spaced aerial interferometer for determining the position of 

a radio star 
Spaced aerial interferometers could be used in several ways to determine the 
position of a radio star. Some of the possible methods are listed below, and an 
indication is given of the way in which they could be applied to determine the 
declination of the star. In all the methods the right ascension is found from the 


local sidereal time of transit ; the transit may be across a true or fictitious meridian. 
In the latter case the value of the declination must be known for computation of Hp. 

(a) The ‘ Periodicity Interferometer”’.—-In the “ Periodicity Interferometer” 
the declination is determined from a measurement of the rate of angular movement 
of the star near transit. If the periodicity of a record due to a star at declination 6 
is denoted by ¢, and that of a record which would be produced by a star in the 
equatorial plane is denoted by t), the approximate relation is given by: 


t =t,secd. 
The quantity f, is calculated from measurements of the aerial spacing (d) and of 
the radio frequency. Since large frequency band-widths are often used for 
improving the signal-to-noise ratio on the record, it may be difficult to determine 
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the effective frequency with sufficient accuracy. ‘This difficulty is overcome by a 
development of the Periodicity Method, which has been called the method of the 
“* Displaced Collimation Plane”. It is described below. 

When the signals received from the two aerials are exactly in phase for all 
frequencies, the star is said to be on the collimation plane of the interferometer. 
‘This plane coincides with the fictitious meridian if the two aerial systems including 
their transmission lines are identical. If the electrical path in the transmission 
line to one aerial is artificially increased, the collimation plane will be displaced to 
one side, and the observed time of transit of the star will change by an amount 
depending on the declination. ‘The angular rate of movement of the star may then 
be observed by measuring the time difference between the transits with the 
increased electrical path first in one aerial system and then in the other. The angle 
traversed is determined by the spacing between the aerials and by the electrical 
length of the inserted transmission path.* The frequency characteristics of the 
receiving system do not affect the measurement of declination by this method and a 
wide band-width may be used. 








Q 
Fic. 3.—The two fictitious meridians of a Crossed Axis Interferometer. The two times of 
transit of a star are separated by an amount T,T, which depends on the declination. (Thu 
projected azimuth errors are taken to be equal in this figure.) 


(b) The Crossed Axis Interferometer.—In the Crossed Axis Interferometer 
two separate interferometers are arranged on the same site with their axes at a 
considerable angle to each other. It is convenient but not necessary to arrange 
the axes symmetrically about a true East-West line. ‘The difference between 
the two times of transit across the two fictitious meridians is now a measure of 
the declination, as may be seen from consideration of the effects of the differing 
projected azimuth and level errors of the two instruments. Differences in the 
projected azimuth errors cause a difference in the times of transit which is 
independent of declination, but the effect of the different projected level errors is 
dependent on declination. This is illustrated in Fig. 3, which shows the 


* This length may be determined by measurement of the frequency at which the extra 
transmission line resonates. 
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fictitious meridians of two interferometers of the same projected azimuth error 
but with differing projected level errors. Such an alignment would be found in 
a crossed axis interferometer on level ground at the equator; more generally 
the projected azimuth errors would also be different. 

(c) The Variable Azimuth Interferometer.—The Crossed Axis Interferometer 
suffers from two practical disadvantages : 


(i) The necessity for the accurate determination of the collimation errors 
of the two separate interferometers. 
(ii) ‘Two complete systems are necessary. 


The use of a single interferometer whose azimuth may be varied by moving 
one of the aerials in a North-South direction has certain advantages over the 
crossed axis interferometer. Since the same aerials are used throughout the 
experiment the determination of declination is not dependent on measurements 
of collimation error, whilst the full gain and resolving power of the whole system 
is used for each observation. ‘The declination is determined by successive 
observations of the star with the interferometer at different azimuths ; the 
change of the time of transit is then related in a known way to the declination. 

In this system the radio star is always observed at the same sidereal time, 
since the azimuth of the individual aerials is unchanged, so that when the azimuth 
of the interferometer axis is changed the maximum of the interference pattern 
corresponding to the transit of the star across the fictitious meridian is no longer 
the maximum of greatest amplitude. It is therefore necessary to make observa- 
tions at a series of different azimuths in order to determine the number of maxima 
through which the interference pattern has moved. Observations over a wide 
range of azimuths may also be of value, since they may sometimes allow the over- 
lapping patterns from neighbouring radio stars to be distinguished. On the 
other hand, considerable mechanical problems are involved in moving large 
aerial structures, and for this reason the method is not applicable to the deter- 
mination of the position of radio stars of very small intensity where aerials of 
large aperture would be necessary. 

(d) The Double Transit Interferometer.—If an interferometer is mounted so 
that its axis makes a considerable angle with an East-West line, it is possible to 
determine the position of a circumpolar star by observing its two times of 
transit across the fictitious meridian (Fig. 4). 

The mean time between the transits is the right ascension +6 hours, and 
the time between them is related to the declination. Furthermore, as will be 
shown in Section 4, the right ascension may be determined without a knowledge 
of the collimation error. 

In this system it is of course necessary for the aerials to have appreciable 
receptivity at both upper and lower culminations; alternatively they might be 
variable in elevation to receive radiation from either direction. 


4. The accuracy of determining the position of a radio star by interferometer 
methods 

(a) General considerations.—In the following sections each of the systems 

described in Section 3 will be examined in turn, and their accuracies in right 

ascension and declination will be deduced in terms of the reading accuracy of a 

single interference maximum, as a function of declination. The sources of error 

in the determination of the time of transit of a radio star will first be examined. 
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(i) Collimation error.—lf the phase delays introduced by the two aerials and 
their associated preamplifiers and transmission lines are unequal, the central 
maximum will occur at a time differing from the time of transit across the axial 
plane; this error corresponds to the ‘collimation error” of visual transit 
instruments. ‘The difference between the phase delays of the two aerials may in 
some cases be determinated by interchanging them, and observing the change 
in recorded time of transit of a star: differences in the electrical length of the 
transmission lines may be determined by direct measurement. It may also be 
possible to measure collimation error by observations of both upper and lower 
culminations of a circumpolar star. 


P 
\/ 


ES 
f Jy 3 


Q 


Fic. 4.—Upper (T) and lower (T’) culminations of a star on a Double Transit} Interferometer. 


(11) Random fluctuations in the record.—The accuracy of a determination of 
the time of transit of a radio star may be seriously limited by the random fluctua- 
tions of the record; the reading accuracy is then determined by the ratio of signal 
intensity to the fluctuation intensity. Increasing the response time of the 
recording device reduces the reading error provided that the response time does 
not become comparable with the periodicity ¢. It is then necessary, however, 
to determine the precise value of the response time, which may introduce other 
difficulties. 

(iti) Confusion between neighbouring radio stars.—In an earlier paper (Ryle, 
Smith and Elsmore (g)) it was shown that the positions of weak radio stars could 
often not be obtained precisely owing to the overlapping of traces from adjacent 
radio stars which might themselves be too weak to be detected; the apparent 
positions would then correspond to the “centre of gravity’’ of two or more 
neighbouring stars. In accurate observations of more intense radio stars a similar 
effect can occur, although the unwanted radio star may have a relatively much 
smaller intensity and may be difficult to detect. ‘Thus a second radio star whose 
intensity is 5 per cent of that of the wanted radio star may alter the phase of the 
interference pattern by 3°, which is equivalent to a displacement in right ascension 
of ;3; of the angular separation of the interference maxima. This source of 
error is examined further in Section 6. 








No. 5, 1952 The determination of the position of a radio star 503 


(iv) Irregular diffraction.—A discussion of the effects of refraction in the 
ionosphere will be given in Section 5. In addition to refraction effects which 
alter the mean apparent position of the source, it has been shown that at times 
irregular diffraction gives rise to rapid fluctuations of the apparent position as 
well as of the intensity. The occurrence of irregular diffraction, which is analogous 
to the scintillation of visible stars, limits the number of nights* when accurate 
observations may be made, and at the longer wave-lengths is likely to be the most 
serious limitation in accurate positional observations. 

(v) Response time of recording apparatus.—The response time-constant (r) of 
the recording apparatus will introduce a delay in the apparent time of transit. 
If r is small compared with the periodicity of the record, this delay is nearly equal 
to tr and may therefore be allowed for. 

In the interferometers discussed here, the reading accuracy is normally 
expressible as a definite fraction of the period of the record produced by the star, 
and thus represents a definite angular accuracy of position measurement 
irrespective of the declination of the star. Collimation error (c) also causes a 
fixed angular error, and an uncertainty in the value of c enters into the positional 
error in the same way as does the reading error. 

The following terminology will be used: 

Azimuth error 

Level error 

East-West aperture of each aerial 

North-South aperture of each aerial 

Separation of aerials 

Reading accuracy of a single maximum expressed as a fraction 
of the period 

Reading accuracy of a single maximum expressed as an angle 
subtended at the interferometer 

Collimation error 

Accuracy of measurement of collimation error 

Angular accuracy of a complete transit observation 

\ccuracy in right ascension 

Accuracy in declination 


In general the angular accuracy Aé is given by 


rA\2 b 
gy2—(—— ) — + (Ac) 
(Aé) (3) 7 + (Ac)?, 


where the reading accuracy r is, of course, itself a function of the aerial dimensions 
as well as of the intensity of the radio star. 
For an interferometer aligned approximately East-West, 
Az =(A@) seco. 


As will be seen from the general considerations outlined above, it is preferable 
to deduce Aé@ rather than Az for position determination with an East-West 
interferometer. 


(b) The Periodicity Interferometer 
(i) Accuracy in right ascension.—The accuracy Az is given by: 


: \\ 2 h 12 
Ad = (Ax) cos 8 ~~ (=) “ n (Ac)? f 





* The effect has been found to be greatest in the night (Ryle and Hewish (8)). 
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‘The accuracy Ac with which the collimation error may be determined depends 
on the type of aerial system employed. If the two separate aerials may be inter- 
changed it is possible to determine c directly; however, this is not normally 
possible as the aerials are either too large to move bodily or else have such a broad 
receptivity pattern that radiation reflected by the ground is added to the direct 
radiation, so that the resultant phase angle depends partly on the contours of the 
ground and is not reversed when the aerials are interchanged. A value of c is 
most conveniently found from observation of an intense source whose position 
has already been found by an independent method such as the Double ‘Transit 
method. If Ac may be made small compared with the angular reading accuracy, 


then 
wn [2 
d d 


(ii) Accuracy in declination.—A number (d/b) of maxima are timed each with 
an angular accuracy rA/d during the movement of the radio star through the 
radiation pattern, which covers an arc A/b radians. ‘The periodicity t of the record 
may thus be determined with an accuracy Af given approximately by: 


At ra 4 a b\ 38 
td (A/b)\/(d/b) \G) ° 


‘ 2 gn 
Aé =~ 4r (3) coto. 
d 


The accuracy Aé is thus related to A@ by: 
A6 


Hence 


b 
AV) = . cotd. 


This relation is shown in Fig. 5 where the aperture of each aerial, 6 = 10A. 


Aé 
Aé 


80 











0 -6° -3° © 3° GF 03 

Fic. 5.—Periodicity Interferometer. Accuracy AS in declination expressed in terms of Aé, 

the angular accuracy of a transit observation, as a function of declination 8 (b= 10d). 

(c) The Crossed Axis Interferometer and the Variable Azimuth Interferometer. 
The analyses of these two methods are similar, and they are considered together. 
For the Variable Azimuth Interferometer the following analysis therefore 
corresponds to the use of only two separate aerial positions, whereas several 
positions would in practice be used ; the resultant error would therefore be reduced. 


* It should be noted that r is itself a function of 6, since it depends on the area ab. 
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(i) Accuracy in right ascension.—A single observation of the transit of a radio 
star locates it on a great circle with an accuracy Aé. ‘This great circle will not 


necessarily coincide with an hour circle, but the angle between them will in 


practice usually be small and the accuracy in right ascension may therefore be 
considered to be 


Ax = A@ sec d, 
as in the Periodicity Interferometer. 
(ii) Accuracy in declination. —The local sidereal time of transit © is given by: 
©=a+Hp, 
where H, is a function of the alignment of the interferometer and is given by: 


sin(m— H,) =tan6tann. 


Vy, 











90° —60° —30° 0 30° 60° 90° 3 


Fic. 6.—Crossed Axis and Variable Azimuth Interferometers. Accuracy AB in declination 
expressed in terms of AO, the angular accuracy of a transit observation, as a function of declination 3. 
{Azimuth difference=0-1 radian.) The relation is shown for latitudes 6=0° and ¢=45°. 


In practice the interferometer axes will both lie in an approximately horizontal 
plane at the latitude 4 of the observing station, so that both m and n will differ 
for the two interferometers. ‘The time interval between the two transits of a 
star will thus be made up of a part due to the change in m, this part being inde- 
pendent of declination, and a part due to the change in , which will be dependent 
on declination. An analysis valid for small differences of azimuth between the 
two interferometer axes shows that 

dH, 
dk 
and it may be shown that for this case the accuracy in declination is given by: 
cosé Ad 
Ad = ——— 


cos k,—k,’ 


- sind —tandcosd, 
, 


where k,—k, is the azimuth difference between the two interferometer axes. 
If the difference k,—k, is produced by a movement x of one aerial then 


k, —k, =x/d so that 
. cosd ra Db 
yaad a 
cosd x d 


The graph in Fig. 6 shows how A8é/A@ varies with declination for interfero- 
meters at latitudes 6 =0° and¢=45°. The value of k, —k, is taken as o-r radian. 


» 
30 
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(d) The Double Transit Interferometer.—Observation of each time of transit 
locates the star on a great circle; the two great circles are equally inclined to an 
hour circle at an angle related to the projected level error of the interferometer 
(Fig. 7). 

The effect of collimation error is also shown in the figure: the fictitious 
meridian becomes a small circle parallel to the axial plane, and the apparent 
declination is changed while the measured right ascension is unaffected. A colli- 
mation error ¢ results in an error in declination equal to ccosecy, where ¢ is the 
angle between the fictitious meridian and the hour circle passing through the star. 
When collimation error is zero, the declination is given by: 


T’ — T=12 hours — 2(m— H,), 
where 
sin(m— H,)=tanodtann. 





Fic. 7.—Intersection of a fictitious meridian with hour circles H and parallels of declination D 
at upper culmination (FM) and lower culmination (F’M’). The broken lines show the effect of 
collimation error. 


From Fig. 7 it is seen that if the reading accuracy of one transit observation 
is Aé, then 
Aa = \/2(A@)cosysecd, Ad = 1/2(A@)cosec ys, 


where is given by: 
sin ys =sinn seco. 


If ys =45°, the area defined by the limits of accuracy is least, and 
(Ax) cos 5 = Ad = 248. 


This degree of accuracy is of course only attainable if the alignment and 
collimation error of the interferometer remains constant between the two transits. 

(e) The Lloyd’s Mirror Interferometer.—The use of this interferometer was 
first described by McCready, Pawsey and Payne-Scott (4), and in greater detail 
for observations of radio stars by Bolton and Stanley (3). ‘The technique allows 
observation of the times of rising and setting of a star at the same or at separate 
observing stations overlooking the sea. ‘These times are deduced from the 
observed times of maximum signal as the star passes through a series of inter- 
ference maxima above the horizon. ‘The mean of the times of rising and setting 
gives the right ascension of the star, and the time interval is the semidiurnal are A, 
related to the declination of the star and the latitude of the station by 


tand=cos A cot¢. 
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Analysis of the accuracy of this method is here made on the assumption that 
the actual elevation of the source at the observed time of rising is known to be 
zero with an accuracy of A@, corresponding to the accuracy of a complete transit 
observation in the previous methods. In practice observations made at small 
angles of elevation are subject to more serious refraction effects in the terrestrial 
atmosphere than is the case for observations near the zenith; the determination 
of the absolute position of a radio star by this method may therefore be limited 
by the accuracy with which the observations may be corrected for refraction. 

The accuracies of determination of right ascension and declination are: 


Ax = 1/2(A@) sec’? (¢ +8) sec!* (6 —5), 
Ad = 1/2(A0) cos 8 cosec d. 

The graphs in Fig. 8 show the variation with declination of the angular 
accuracies (Ax)secd5 and Ad expressed in terms of \@, for latitudes ¢=0", 30° 
and 45°. 

(Aa) sec 64 
0 
A sue 
Nw) “i 
Benet 2 
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Fic. 8.—Lloyd’s Mirror Interferometer. Angular accuracies (Ax) sec @ and Ad tn terms 
of the angular accuracy A®é of the observations of both rising and setting. 


It is possible to determine the position of a radio star by observing the rate 
of change of its elevation, together with its time of either rising or setting. It is 
usually possible, when a single site is in use, to observe at only one of these times, 
and the resulting accuracy of position-finding may be somewhat worse than that 
obtained from observations at both rising and setting. 

‘The method is analogous to that employed in the Periodicity Interfesometer, 
in that an angular rate of movement is derived from a measurement of the 
periodicity of the record obtained when the star is near the horizon. ‘The 
accuracy of measurement of this periodicity is given in terms of the reading 


accuracy Aé by: 
At b 
eee). Ll 
7; = (;). 


where 5 is now the vertical aperture of the aerial. ‘Then it may be shown that 


Aé = : b Adé 
os: tan (6+ 4)—tan (o—¢) r : 





360* 
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This relation is shown graphically in Fig. 9 in which the aperture 6 is taken 
as 5A. 

As the accuracy in declination is usually considerably less than that of the 
time of rising (or setting), the star is located by this method inside an area which 
is considerably elongated. The direction of elongation is not necessarily parallel 
to an hour circle, so that observations made with East-West and cliff interfero- 
meters may be used together for position finding without the measurement of 
periodicity in either. 


A> 49+ 
Aé 











—9° -60° —30 0 90° 8 


Fic. 9.—Lloyd’s Mirror Interferometer. Accuracy Ad in declination expressed 
in terms of @ for observations either of rising or of setting of a radio star (b= 5A). 


5. The effects of atmospheric refraction on the apparent position of a radio star 


Appreciable refraction may occur at radio wave-lengths both in the troposphere 


and in the ionosphere. ‘The deviations in the apparent position of a radio star 
may be considerable; the refraction in the ionosphere may in fact be sufficient 
for measurements to be made of the total ionization of the F region. 

(a) Tropospheric refraction.—Normal tropospheric refraction has the effect 
of making the star appear nearer the zenith than it really is. ‘The refractive 
index of dry air is practically the same for radio waves as for light waves, and the 
refraction effects therefore correspond closely to the standard astronomical 
effects. ‘The refraction angle R may be expressed for dry air at N.'T.P. as 


R=56":3 tan Z —0"-005 tan*® Z, 


where the zenith angle Z is less than about 75°. 

The refractive index of water vapour, however, is greater at radio wave-lengths 
than at visual wave-lengths, and at large zenith angles appreciable deviations 
ray be caused by a rapid gradient of refractive index caused by a rapid lapse 
rate of water vapour. ‘lhe problem of refraction produced by different vertical 
distributions of water vapour has been dealt with in detail in the Physical and 
Meteorological Societies’ Report (1947), and it may be seen that a detailed 
knowledge of the vertical profile of meteorological conditions in the lower tropo- 
sphere is essential for the computation of refraction angles at large zenith angles. 

The magnitude of the correction that must be applied is illustrated by an 
example quoted by McCready, Pawsey and Payne-Scott (4), in which the total 
atmospheric refraction at an elevation of 5° was calculated to be 15’, of which 5’ 
was due to the presence of water vapour. ‘Their measurements on solar radiation 
at 200 Mc/s agreed well with this value. 
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(b) Jrregular ionospheric diffraction.Ryle and Hewish (8) have shown that 
the rapid fluctuations in the signals from radio stars are associated with irregularities 
in the ionosphere, and that they are accompanied by changes in the apparent 
position of the radio star. ‘The position may vary in a period of a few seconds 
by 1-2 minutes of arc on 3°7m wave-length; the magnitude of the effect is pro- 
portional to the square of the wave-length. This limitation to the accuracy of 
measurement of position is of considerable importance, since it limits the number 
of days and the time of day at which accurate observations may be made. 

(c) Regular ionospheric refraction. Steady deviations of the apparent position 
of a star from the true position may also be caused by ionospheric refraction. 
‘l'wo causes may be distinguished : 

(i) A horizontally stratified layer over a curved Earth. ‘The apparent zenith 
angle of a radio star is increased because of the curvature, and the apparent 
declination as measured by a Periodicity Interferometer is decreased even at 
vertical incidence. 

(ii) Departures from horizontal stratification of a layer. A_ horizontal 
gradient in ionization may be caused by diurnal variations in ionization density, 
and the resulting small-angle wedge will alter the observed right ascension of a 
radio star. 

‘The magnitudes of these steady refraction effects are considered in Appendix B, 
where it is shown that they are primarily determined by the total ionization per 
unit column of the ionosphere. Measurement of these effects may prove to be 
of value as a method for exploring the ionosphere (Bailey (1); Smith (11)). 


6. The effect of neighbouring weak radio stars on positional determinations 
using interferometers 


‘Lhe record may be considered as consisting of a sinusoidal trace of amplitude 
R, from the wanted star, together with smaller sinusoidal traces from all other 
stars in the receptivity pattern of the aerials. ‘These smaller traces are here 
considered as having the same frequency as that of the wanted signal, although if 
stars are received simultaneously over a wide range of declinations this will not 
be so; in such a case only those traces with a frequency close to that of the wanted 
trace will be likely to cause confusion. 

If the probable distribution in intensity of these interfering stars is known, 
it is possible to estimate the most probable amplitude R, of the resultant sinusoidal 
trace on the record. 

If R is the record amplitude from the mth star in order of decreasing intensity 
in the solid angle of the aerial receptivity pattern, then it may be shown that the 
most probable value of R is related to n by 


R°?n = K, 


where A is a constant (Ryle (5§)). 

Then An=}KR °*AR, where An is the number of stars giving records with 
amplitudes between R and R+AR. The resultant amplitude of the sinusoidal 
traces, which will add in random phase, is then R(An)"?. The square of the 
resultant R, is the sum of the squares of all these contributions, and R;, is there- 
fore given by: 


R,=[f{R* dn", 





510 F. G. Smith 
where the integration includes all stars of intensity less than Ry. Hence 
By 1/2 
-| | IKRAUR | 
“UV 
=(3KRs!™}, 
therefore, = 4/(3p,), 


where p, is the probability of finding a star giving a record of amplitude Rg or 
greater in an area of sky equal to the solid angle over which the aerial is receptive. 
When the ratio R,/Rg is considerably less than unity the effect of the resultant 
sinusoidal trace R, on the measured right ascension of the star may be simply 
investigated. ‘The phase modulation of the sinusoidal trace Rg due to the addition 
of a sinusoidal trace R; in random phase has a maximum value of R,/Rg radians. 
The r.m.s. value is 1//(2)(R,/Rs) radians, representing an error in right 
ascension of 
I I R; 
on . V2) ‘ Rs as 
where ¢ is the periodicity of the record. 
For example, in a recent determination of the positions of four radio stars 
(Smith (z0)) an aerial system was used which at a wave-length of 1-4m was 
1 


receptive over a solid angle of about ; steradian. ‘The probability of finding a 


radio star of the intensity of the radio star in ‘Taurus (05-01) in this solid angle 
is about ;. The lobe spacing of the interferometer was 17’ arc, so that the 
resultant uncertainty in the right ascension of this star was about two seconds. 
For the more intense radio stars (19°01) and (23-01) the uncertainty was less than 


} second in the right ascensions. 


7. The application of the different types of interferometer to various observations 
in positional astronomy 


It is important to decide which is the best way of using an interferometer for 
a given purpose. In discussing this question it is convenient to divide the 
positional observations of radio stars into two categories, which may be termed 
accurate position determinations and survey observations respectively. The 
former are normally only possible for a few of the most intense stars, whereas the 
object of the latter observations is primarily to detect the maximum possible 
number of radio stars, determinations of the position of most of the stars observed 
being necessarily inaccurate owing to their low intensity. 

(a) Accurate determinations of the positions of intense radio stars.—For 
observations of the apparent positions of radio stars of low declination the 
Lloyd’s Mirror Interferometer has many advantages; the considerations of 
Section 4(e), however, have shown that for accurate observations over a consider- 
able area of sky the latitude of the site must be in the region of 30° to 45° (see 
Fig. 8). There are considerable practical advantages both because only one 
aerial is used and also because there is no need to determine a collimation error. 
This interferometer has been used to observe radio stars between declinations 
+ 40° and — 40° (Bolton and Stanley (3); Stanley and Slee (12)). 

With a spaced aerial interferometer observations can be made when the 
zenith angle is least, which minimizes the corrections for refraction; and a phase- 
switching system (Ryle (6)) can be used so as to provide considerably increased 
accuracy in reading. ‘The choice between the various methods of using spaced 
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aerial interferometers to measure declination largely depends on the range of 
declination in which the star lies. ‘The Displaced Collimation Plane method, 
applied to a Periodicity Interferometer, provides reasonable accuracy in 
declination for radio stars of declination greater than 30°. For circumpolar 
stars the Double Transit method allows an independent determination of 
declination provided that collimation error may be determined by the inter- 
changing of the individual interferometer aerial elements. Both these methods 
have recently been used by Smith (ro). 

For stars of low declination the Periodicity Interferometer is inaccurate. 
The use of the Variable Azimuth and Crossed Axis Interferometers would then 
appear to present the most satisfactory solution. 

(b) Survey observations. Observations designed to detect radio stars of small 
intensity require the use of large aerial systems. It is not normally possible to 
determine the azimuth, level, and collimation errors of such an aerial as accurately 
as is desirable for absolute determinations of position, but it may still be possible 
to determine the relative positions of a large number of stars with considerable 
precision. 

The practical difficulties of constructing a large aerial for use in a Lloyd’s 
Mirror Interferometer are considerable, as it must be rotated about a vertical 
axis for observations on different parts of the horizon. ‘T‘he spaced aerial inter- 
ferometers, on the other hand, can make use of large East-West apertures in the 
aerials without the necessity for rotation of the aerials about a vertical axis. 
Various ranges of declination may be covered either by making the North-South 
aperture small, or by rotation round a horizontal axis. Aerials with appreciable 
aperture in a North-South direction could be used to determine declination by 
the use of their resolving power in this direction to an accuracy which might be 
sufficient for survey purposes, but this will not be considered here. 

For declinations greater than about 30° the Periodicity Interferometer clearly 
has advantages for survey observations. At lower declinations a further 
possibility is presented by the use of the Crossed Axis Interferometer which has 
considerable accuracy over a wide range of declinations including low declinations. 
Positional measurements with this system do not depend on movement of aerials 
as in the Variable Azimuth Interferometer, and large structures may therefore 
conveniently be used. 


APPENDIX A 
1. The effect of azimuth and level errors.—The “ projected azimuth error” 


(m) and “ projected level error’? (m) are related to the azimuth error (&) and 
level error (7) by the equations: 


sin m=sin/sind + cos/sink cos 4d, 
sin m cos n =cos7sink sind — sinicos 4d, 
cos m=cosn=cosicosk, 
where ¢ is the altitude of the interferometer. 
2. The time of transit of a star across the fictitious meridian of a skewed 
interferometer.—The hour angle H of a star when the received signal is a 
maximum (see Section I) is given by 


sin (H —m) = 


sec d secn—tandtann. 


na 
d 
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The central maximum thus occurs at a sidereal time © given by 
O=a+H,, 
where H, is given by 
sin (m— H,)=tandtann. 


APPENDIX B. Jonospheric refraction 


1. Steady refraction in uniform ionospheric layers.—In optical positional 
astronomy there is an important refraction effect due to the departure from unity 
of the refractive index of air. ‘The magnitude is complicated by the curvature 
of the Earth. The refractive index of dry air is the same for radio waves as for 
light, but additional refraction occurs in the ionosphere. As the ionization at 
ground level is negligible, deviations of apparent position are only caused by 
the curvature of the ionospheric layers. ‘The analysis of the effect is therefore 
an extension of the detailed treatment of optical astronomical refraction. 

It may be shown that the refraction angle R suffered by a ray in passing from 
one spherical shell with a refractive index » to one with refractive index + yu 
is given by 


SR= bu ty SIN ZF 





(w= po" sin? 2+ 252)?’ 
where /9 = refractive index at Earth’s surface, 
= apparent zenith angle of star, 
= height above Earth’s surface in terms of the Earth’s radius as unity. 

‘This expression applies only where s is small. For the ionosphere py =1 
and =1—Ap, where Apu = }(f,/f)® and is small (except when the frequency f is 
near the critical frequency f,). For Ayu small, 
OV. sin 2 
apo t__sins__ 
pe (cos® z+ 2s)"? 
The total refraction angle R is then given by 
R= 





sin = 


du, 


jt (cos? z+ 2s)" 
where the integration is taken over the depth of the ionosphere. 
Restricting the analysis to zenith angles where cos? z is large compared to 2s, 


q $ 
R= franz (1- - ) a. 
cos* z 


The tirst part of the integral is the part which applies to plane layers; it is 
zero when , is unity at ground level and outside the ionosphere. (This part 
gives the amount of tropospheric refraction R =) tan z referred to previously.) 

The second part is due to the curvature of the layers: 

R= —tanszsec® zfsdp. 

Since x =1—kN, where k =e? 27mf? and N is the electron density, it follows 

that 


kn 
R=—tanzsec?s—, 
. 


where » is the total electron content of the ionosphere per unit area and r is the 
radius of the Earth. The effect of this type of refraction is to make the star 
appear farther from the zenith than it really is. 
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The order of magnitude of the refraction angle is illustrated in ‘Table I 
below. This refers to an ionospheric layer of parabolic distribution, with halt- 
thickness = 100 km, and penetration frequency 10 Mc/s. In this layer there are 
°6 x 10 electrons per sq.cm column. 


TABLE | 
Refraction angle R (minutes of arc) as a function of frequency and zentth angle for 
an tonospheric layer in which there are 1°6 x 10'* electrons per sq. cm column 
Zenith angle 10 20 30 40 50° 60 
Frequency 
in Mc s Refraction angle in minutes of ar 
40 0°94 ‘74 3°26 
so 0°24 49 ‘$2 
160 0°06 12 ‘20 
320 is oo! 03 05 


‘The accuracy of the measurements of declination which have so far been 
described has not allowed any measurements of this type of refraction to be 
made. It should eventually be possible, however, to measure the total electron 
content of the ionosphere in this way. It has been shown by Bailey (1) that 
detailed measurements of this effect at large zenith angles might lead to knowledge 
ot the electron distribution above the maximum of the F layer. 

2. Wedge refraction in an ionospheric layer..-Suppose that the total ionization 
x changes linearly in East-West direction with a gradient dn dx. ‘Then it may be 
shown that the refraction angle R is given by 


R=k 


where k= ——- as before. 


d. 2nmf* 
Some observations of the effect of this wedge refraction will be published 
elsewhere (Smith (11)). 


dn e 
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ON SOME POSSIBLE MECHANISMS OF RADIO S'TARS 


I. ‘HE CRAB NEBULA 


F. D. Kahn 
(Received 1952 August 7)* 


Summary 


It is assumed, first, that the radio noise emitted by the Crab Nebula 
is due to the deflection of fast electrons (with a velocity of about 10° cm/s) in 
fields of force with atomic dimensions. ‘The total energy output necessary for 
this process is found to be improbably large. 

It is further shown that it is not possible to explain the observed intensity 
of radio noise by free-free transitions due to electrons with thermal velocities 
only. 





1. Introduction.— Most of the radio stars for which observations are at present 
available are thought to be objects of small angular diameter. ‘The intensity of 
their radiation is found to vary inversely as the frequency. An exceptional 
object, however, is the Crab Nebula, in which the intensity of noise emitted is 
apparently constant with varying frequency, and which, in the visual spectrum 
at any rate, subtends at the Earth an angle of some minutes of arc. These results 


are due mainly to Stanley and Slee (1). 

It is tempting to jump to the conclusion that the Crab Nebula is an optically 
thin cloud of hot gas emitting radio waves by the free—free transitions of electrons 
in the field of ionized atoms. In such a process the intensity emitted is very 
nearly constant per unit frequency interval for sufficiently low frequencies. It 
will be shown that this mechanism is hardly likely to be efficient enough. In fact 
any reasonable mechanism which is based on the emission of radiation from 
independent accelerating charged particles requires an improbably large output 
of energy in the nebula as a whole. 

2. The radiation from an accelerated electron.-It is not difficult to derive a 
formula for the intensity of the radiation emitted by an electron which suffers 
a sudden change in velocity. Consider a particle which, referred to suitable 
axes Oxy, has velocity components (u, —v) from time t= — o until t=0, when 
the velocity suddenly changes to have components (u,v). ‘The acceleration of 
the electron at time t becomes 

f(t) = 2v8(t) 
parallel to Oy, where 8(t) is the Dirac delta function. Now f(t) is an even function 
of t, hence we may write 
~ 00 


f(t) = 2v8(t) =(2 mua a,, cos nt dn. 


“0 
oe at 
I'he Fourier components a,, are 


a,=(2 7)! 2 2vd(t) cos nt dt = v(2/7)*?. 
~ 0 


* Received in original form 1952 May 10. 
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The energy radiated in the frequency interval between m 27 and (n+ dn) 27 
is easily shown to be 22 
U(n) dn = —, . — . 2ma,? dn, 

3¢ 7 
oe 166772 
—; a,* dn = —— dn, 
3¢ 3rc8 
where e = electronic charge, 
c= velocity of light. 

‘The total energy radiated is, of course, infinite. However, an instantaneous 
change in velocity, such as is described here, cannot occur in nature. Suppose 
then that an electron enters a deflecting field which is such that the acceleration 
parallel to Oy at time ¢ is given by f=/(t), where f(t) is an even function of f, and 
that ro - oO 

20 =| fit)dt=2| fit)dt. 
-o 0 

We require further that f(t) does not become infinite at any time, that it 

reaches its largest value when t=o and decreases monotonically to zero when 


t>o. It follows that | | f(t)|? dt converges. ‘The total energy radiated by the 


electron is now ee ne 
é = (22/304) | | f(t)? dt =(4e? 3c3)| [| f(t) |? dt. (1) 
tant J0 


Once again let sa 
f(t) =(2 nr | a,, cos nt dn, 
“0 

@ 
then a,, = (2/7)! | f(t) cos nt dt 


0 


and Ja, | <(2/m)'2[ f(t) dt =0(2/m)"*. 


16€2v? 


3nc8 


As before 2s 
U(n) = — 34," < 


(and it may be shown that the equality sign is approached as n-> 0). 
With the aid of (1), the inequality (2) becomes 


qa © 2 - 
fg A t, 


[nora ” 


P v? ; : ' ; bie sre 
where t = —————— , and 1s a measure of the duration of the collision in which 


“x 


| | f(t) |? at 


the electron is deflected. 
For consider the curve f = f(t). Since z f(t) dt, v is the area under the curve 
a 
between t=0 andt= %. On the other hand, if fis the distance of the centroid of 
2S 
=] | f(z) |? dt 
° . ‘ 3 2 0 ‘ ba v~ Uv 
this area from the t-axis, f= ————— : “nce ———— Yi so that 
v an 2 
| | f(t) |? dt 
0 
t is the width of a rectangle which has sides parallel to Ot and Of, has an area 
equal to v and whose centroid is at the same distance from Of as that of the area 
under the curve f =/(?). 
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‘To give two examples: if f=(2v/\/m)e, the value of # defined in this sense 
. es a 2U - i ee ‘ , ; : , : 
is y (7/2). If f =———., t=z. Thus 21 gives an indication of the time which 
: am(1+ ft) 
the electron spends in the deflecting field. 
Expressed in terms of the frequency v, relation (3) becomes 
U(v) =27U(n) <8 61. (4) 


3. Application to the Crab Nebula.—Suppose now that the radiation from 
the nebula were due exclusively to fast electrons. Each time an electron is 
deflected in a field of force, of linear dimensions /, it emits electromagnetic waves. 
The electron has a velocity of the order of 10!° cm/s, so that the period 2¢ in 
which the change of velocity takes place is given by 2#~10°'/. If a total 
energy & is radiated in the encounter then U(v), the energy density per unit 
frequency interval, obeys the inequality 

U(v) <8té =4 x10 C6 
or 
6 >2°5 x10%" U(r). (5) 

Now Stanley and Slee (1) find that at 100 Mc’s the intensity of the radiation 

from the nebula received at the Earth is 

2 x10 73 watt metre? (c/s)! =2 x10°" ergs! cm=*(c,s) 1. 

If the distance to the nebula is 10° parsec = 3 x 107! cm, then the total intensity 
of radiation emitted by it is 2 x 104 erg s™! (c/s) '. By formula (5) the total energy 
radiated at all frequencies is at least 2-5 x 10%! times as large and therefore 
exceeds 5 x 10/1 erg/sec. With /~10~-% cm this becomes 5 x 10*' erg/sec. 
Now the mass of the Crab Nebula is estimated to be 15 Mg =3 x 10**g, so that 
the energy output required must exceed 10’ erg g''s"'. If the Crab Nebula has 
existed 1000 years, i.e. since the supernova of 1054 A.D., each gram of matter 
in the nebula must have liberated at least 3 x 10'’ erg. This appears to be most 
unlikely. 

A comparison may be made with a hypothesis which has recently been put 
forward by W. H. Ramsey (2), and according to which the Crab Nebula derives 
its energy from the radioactive decay of *C atoms. ‘These atoms have a half-life 
of about 10" sec, and the energy of their decay is 1-5 x 10° eV, or 2-4 x 10 7 erg (3), 
most of which is given to the f-electron. This corresponds to a release of 10'* erg 
per gram of 'C, and is clearly not sufficient. If, as seems probable, only about 
OI per cent of the nebula consists of *C atoms, the difficulty becomes much 
yreater. 

We must conclude that radiation in atomic fields by fast electrons is of little 
importance. 

4. The efficiency of free-free transitions.—-It is still possible that free—free 
transitions of thermal electrons in the field of ionized atoms may provide sufficient 
radiation to explain the observations. It will be recalled that the nebula is 
optically thin, and in order that a thermal process may provide sufficient radiation, 
the actual temperature must exceed 2-5 x 10® deg. K, the equivalent black-body 
temperature. 

‘To take the most favourable case, suppose that the gas consists entirely of 
fully ionized atoms of some element with atomic number Z. Let there be » such 
atoms per cm® and let V be the volume of the nebula. The atomic weight of 
this element will be approximately 2Z, and if « is the mass of an atom 

iu 2Zmy, M=np V= 2ZnmyV, 
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where 
my = mass of a hydrogen atom, 
M=mass of the nebula. 
Hence 


VW 
The corresponding electron density will be 
M 
2m_V" 


n = Zn = 


But at frequencies v, small compared with kT/h, the energy spontaneously 
radiated per unit frequency interval per unit time in a unit volume of gas at 
temperature 7’ is, by the usual formula, 


ony Zee Bet (m\' (AT 
a ia V3 ma \kRT "\iy }° 


Since we are attempting to find an upper limit to the intensity of radiation 
emitted, there is no need to consider stimulated emission, which is always 
less important than absorption. Insertion of numerical values gives, for 
v=100 Mc/s = 108 c/s, T > 2-5 x 108 deg. K, 

J <8-2 x10-"n,nZ?. 
Che total energy output, per unit frequency interval, for the whole nebula 
is then 
8-2 x 10° 4° M?Z 
4my"V 
following relations (5) and (6). With M=3 x10* g and my =1-6 x10 * g 
JV <82 x107Z/V ergs (c/s)"!. 


‘The observations require that 





JV <8-2 x 10°"n,nZ?V 


’ 


JV =2:2 x 10%, 
so that V/Z<4 x10, 

Taking Z=10 as a reasonable value for the atomic number of the gas 
comprising the Crab Nebula we find that V <4 x 10°* cm’, so that the radius of 
the nebula must be less than 4-5 x 10'7 cm =} parsec, or, at least, the bulk of the 
matter must be concentrated in a small region such as this. Since the radius 
of the nebula is observed to be six times as large as our upper limit this inference 
contradicts the experimental results. 

It follows that the radio emission remains unexplained whatever the energy 
sources of the nebula may be. It will perhaps be necessary to consider some 
process which transforms mechanical energy of the turbulent gas into radiation. 
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ON ‘THE FORMATION OF INTERSTELLAR DUST 
F. D. Kahn 
(Received 1952 August 27)* 
Summary 
A process is suggested which leads to the formation of elongated and 


electrically polarized ice grains in interstellar space. Some values are found 
for the dipole moment to be expected on grains of different sizes and shapes. 





1. Introduction.—'Vhere is, at present, little agreement among astrophysicists 
about the nature and composition of the interstellar dust grains. ‘The recently 
discovered polarization of light from stars in certain regions of the galaxy is 
almost certainly due to the grains (1, 2) but no mechanism to explain it has as 
yet been generally accepted. Spitzer and ‘Tukey (1) have shown how the 
polarization may be explained if it is assumed that a large number of grains 
consist mainly of ferromagnetic material. ‘They have suggested a series of 
reactions which can produce such grains but, for a satisfactory explanation, they 
require a rather intense magnetic field to exist in interstellar space. 

Davis and Greenstein (2) have developed a theory which requires dust 
grains consisting mainly of dielectric material, with an admixture of iron. ‘The 
grains have an elongated shape and are kept in rotation by collisions with the 
free atoms in space. ‘They show that a frictional couple acts on each grain owing 
to the interaction between the iron atoms and the interstellar magnetic field. 
This couple ultimately tends to set the grain so that its principal axis with the 
greatest moment of inertia is parallel to the field, and so that the grain rotates 
about this axis. In this paper it will be assumed that the theory of Davis and 
Greenstein is correct. 

‘The main object of the present work is to show how a dust grain may be 
formed so that it possesses a large electrostatic moment. Owing to this the 
grains will be very elongated, and their existence can be shown to lead to the 
formation of a source of galactic radio noise (3). 

2. The water molecule and its properties.—According to Harrison Brown's 
table of cosmic abundances (4) oxygen is the third most common element in the 
universe, the first and second being hydrogen and helium. ‘The next most 
abundant elements are nitrogen, carbon and iron. It is, therefore, reasonable 
to assume that the material in a dust grain will consist largely of ice (H,O), with 
smaller admixtures of ammonia (NH;), methane (CH,), iron and its oxides. 
This was pointed out by van de Hulst (5). Any peculiar property of the H,O 
molecule will have an important effect on the nature of the grain. 

The water molecule has an electric dipole moment of 1-8 x 10-8 e.s.u. 
A grain of ice in which the dipoles of the constituent molecules are parallel will 
therefore have a large electrostatic moment (equal to 6 x 10‘ e.s.u. per cm® in a 
pure specimen). This is not found in ordinary ice, but when a crystal grows 
in a sufficiently strong electric field it is possible for it to be fully polarized. 


* Received in original form 1952 May 17. 
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‘The peculiar nature of ice formed in an intense field was first noticed by 
Bleekrode (6). He found that when a crystal of tourmaline was cooled in liquid 
air a strong electric field was set up at each end. When the specimen was 
removed from the liquid air, long crystals of ice condensed at each end and gave 
it the appearance of a magnet dipped into iron filings. In the absence of a field, 
the ice freezing on to a cold surface was found to form a smooth coat. 

The German meteorologist Rossmann (7) has postulated that polarized ice 
crystals exist also in thunder-clouds and has succeeded in growing a polarized 
crystal in a refrigerator. ‘The conditions necessary for the formation of such 
crystals will be examined later and will be shown to exist in interstellar space. 

According to Lennard-Jones and Pople (8, g) a water molecule is approximately 
spherical with a radius of 1-38 A. ‘The oxygen atom (Fig. 1) is at O, the centre of 
the sphere, and the hydrogen atoms are at A and B on the surface such that 
ZAOB~109° =cos-!(—4). Near A and B are regions of high potential which 
will attract negatively charged particles. ‘The orbits of the two pairs of electrons 
which do not take part in the formation of the OH bonds are repelled to the far 
side of the molecule, forming a region of high electron density near the points 
C and D of the regular tetrahedron ABCD, whose centre is at O. The direction 
of the dipole is along the perpendicular from O on to AB. 





/ 
/O, 


Fic. 1.—The H,O molecule and its four nearest neighbours. 


In an association of water molecules the end A or B of one molecule will 
tend to be in contact with an end C or D of a neighbouring molecule. If O, is 
the centre of one of the neighbouring molecules the position of smallest potential 
energy is that in which OAQ, is a straight line. 

If the molecule at O, is joined only to that at O it ts free to turn about OO, 
and to oscillate so that 7 OAQO, varies (between 160° and 180° at 273 deg. K). 
In the H I regions of space the temperature is only about 60 deg. K and we shall 
ignore the possible bending of the bonds. 

In any crystal of ice the molecule with its centre at O is surrounded by four 
other molecules with centres at O,, O,, O, and O,, each at a vertex of a regular 
tetrahedron and each 2-76 A distant from its centre O. Nothing has yet been 
said about the dipoles belonging to the molecules with centres at O,, O,, Os 
and O,. However, it is possible to construct a model of a specimen of ice in 
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which all the bonds are formed and are straight, and in which all the dipoles 
point in the same direction. A diagram to illustrate this is given in Fig. 2. 
in this model each positive pole is in contact with a negative pole of a neighbouring 
molecule. All the dipole moments are seen to be parallel. 

Once formed, such a structure will be stable at 60 deg. K. In order to 
destroy it a large number of molecules must break at least three bonds, since 
two bonds with neighbours are enough to hold it fixed. According to Pople (9) 
the energy required to break a bond at 273 deg. K is 11-5 kcal/mole, or 
Sx 10° erg per pair. The structure will therefore persist at 60 deg. K, which 


corresponds to a thermal energy of 10°! erg. 


ae 
“i 











Fic. 2.—Diagram to illustrate an ice crystal with a regular structure. The symbol enclosed 
hy the circle represents one molecule. The positive ends of each molecule are marked with small 


circles. 


3. The growth and strength of a grain.—Unless there is a nucleus on which 
it can begin to condense, an ice crystal cannot form at all.* Inthe H I regions ot 
space the initial condensation may be expected to take place on small grains 
made up of heavier elements which would condense earlier than the lighter 
substances. Such small grains will adsorb a certain number of atoms (mainly 
H atoms, since these are the most common). Any O atom which strikes the 
grain will combine with the hydrogen to form an H,O molecule. In Section 5 
it will be shown that each grain may be expected to have an electric charge, and 
the H,O molecules, by Bleekrode’s experimental result, will tend to collect 
where the field is most intense. ‘The H,O molecules will be held parallel to the 
tield if mE >kT where 

m = dipole moment of the H,O molecule (=1-8 x 10°'* e.s.u.), 

E =intensity of the electric field, 
k = Boltzmann’s constant, 
T =temperature (=60 deg. K). 

* A process has been described by Kramers and ter Haar (B.A.N., 10, Nos. 361 and 371, 1946) 
by which grains can grow without metallic condensation nuclei, but are formed by agglomerations 
of CH and CH? molecules. It is unlikely that our theory will apply to grains which begin their 
growth in this way. 
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The necessary minimum field intensity is then 6 x 10% e.s.u., or 1°8 x I0® volt 
per cm. 

Once the crystal has begun to grow with a regular structure it will continue 
to do so. ‘The positively charged nucleus will be in contact with the negative 
ends of the H,O molecules so that the H atoms become the outer end of the 
growing crystal. ‘Ihe electric field at the outer end will now be approximately 
that due to a surface layer of charge of 6 x 104 e.s.u. per cm? (equal to the dipole 
moment , per unit volume). Any molecule present at the positive end will be 
added to the crystal with its dipole parallel to the existing polarization. 

The field at the positive end is sufficiently strong to hold some free H atoms. 
With a polarizability « =3 x 10°*° cm* the energy of an H atom in a field E is 

-4aE*. In the present case E~47p=7:2 x 10° e.s.u., so that the potential 
energy becomes —7:5 x 10 erg. Since kJ =10 ' erg, some H atoms may be 
expected to remain near this end of the crystal. Any O atom arriving there will 
combine with them and the resulting H,O molecule will become part of the 
grain. Any colliding N or C atom may similarly form NH, or CH, molecules. 
However, these have no intrinsic dipole moment and are therefore less firmly 
held than an H,O molecule. 

It is important to know the exact composition of the grains; and to find it 
one has to consider the effect of this weaker binding on the eventual abundance 
of methane and ammonia. It is hoped to attempt this in a later paper. Other 
important impurities may be expected to be silicon, iron, and possibly molecular 
hydrogen. 

It is hard to say what proportion of foreign molecules can be present before 
the ice-like structure of a grain is destroyed. Since each molecule is held fixed 
by two (out of four) nearest neighbours it will presumably be possible for the 
regular model to exist when up to about one-half of its constituents are impurities. 

A regular ice crystal free from impurities will have a very great tensile strength. 
In the electrostatic model of the O—HO bond given by Lennard-Jones 
and Pople (8), the force of attraction between neighbouring molecules is 
37x 10° dyn. ‘To fracture a grain all bonds in one plane of the crystal must 
be broken. Now the distance between two adjacent molecules whose centres 
lie in a plane perpendicular to the length of the crystal is d= R(2 +4)! where 
R (=2-76 A) is the distance between nearest neighbours in the crystal. Hence 
d= 4°5 A and the number of molecules per unit area in the plane is approximately 
5x10! cm’. In order to fracture the crystal across such a plane ro” bonds 
must be broken per unit area. Since each bond is inclined to the crystal axis 
at an angle cos '(1/\/ 3) the breaking stress becomes 


(3:7 x 10° x 10!)/4/3 ~2 x 10!” dyn per cm’. 


The strength sets an upper limit to the possible speed of rotation of a grain. 

If a typical grain is a prolate ellipsoid of revolution, with axes of length 2a and 
2b (a>6), and rotates with angular velocity w about a minor axis, the internal 
stress at its centre due to centrifugal force will be P/zb? where 

ra 

P= | pA(x)w*x dx 

~0 
where p = density of the grain, A(x) =cross-sectional area at distance x from the 
centre, taken to be small compared with a’. 


37 
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Now A(x) =76*(1 — x*/a*), so that 

P/rb* = tpa*w?. 

‘The grain will not fracture if }pa*w?<2 x 10'° dyn percm®. If mis the mass 
of the grain and J the moment of inertia about a minor axis, J = 4m(a* + b?)~1ma* 
for a highly prolate grain. If V is the volume and @ the kinetic energy of 
rotation we now have 

spa*w* = Blu, V=5 é V 
and the grain will not break as long as 5 &/V <4 x 10" dyn per cm?. 

4. The electrostatic potential on a grain.—Owing to its polarization there will 
be a considerable electrostatic potential (denoted by y,) near the surface of an 
ice grain. From the microscopic point of view this potential is due to the charges 
present at those ends of H,O molecules in the surface which have not been able 
to form any bonds. Strictly speaking one should investigate the electrostatic 
conditions near the H,O molecules themselves. However, this would be difficult 
since the neighbouring molecules will always have a considerable influence at 
any one place. We shall, therefore, take the macroscopic point of view. 

To facilitate the evaluation of the potential we assume that a typical grain 
has the shape of a prolate ellipsoid of revolution with its long axis (of length 2a) 
parallel to the direction of polarization, with minor axes of length 26 and a dipole 
moment per unit volume. 

Let O be the centre of the grain and Ox lie along the major axis. ‘The 
potential due to the grain at any point is equal to E,, which we define as the 
x-component of the intensity due to a body with the same shape as the grain, 
but having a volume density of charge ». Following Kellogg (10) the value 
of E, at interior points is given by 

E,=2Ax 
poe ds 
o (a? +s)3(b? +5)” 





where A =zab* | 


if e is the eccentricity of the ellipsoid 


3 I—e? I+e , 
E,=27px - | in (=<) - 2e| =2mpxF(e), Say. 


‘The potential is greatest where x=a. Some values of F(e) are given in 
Table I, for large e, or small b/a | =(1 —e?)"?}, 
TABLE | 
b a F 

o'l 0040 
0°032 0°006 3 
oO'o! 0000 86 
0°003 2 0°000 II 


Some typical values of yp at the ends of grains with different volumes and 
eccentricities are given in ‘Table II. 

In the case of all ellipsoidal grains (such as we are considering) the potential 
will be positive on half the surface and negative on the other half, the dividing 
plane being x=0. ‘To save troublesome calculations later we now make the 
following approximation. If, for a given grain, the potential at the end x =a has 
been found to be xo, then that grain will be treated as though the potential had 
the constant value x) on the whole of the positive half, and the value — yx, on 





>? 
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the whole of the negative half. ‘This approximation is not too violent, for it is 
easily shown that half the positive surface charge of a grain resides on the region 
for which a/\/2<x<a, and on this portion the potential varies by a factor of 
\ 2 only. 

5. The ionization equilibrium in the H I regions. Its effect on the grains and on 
the condensation nuclet.—The growth of the grains is assumed to take place in 
the HI regions of space. ‘The free hydrogen atoms here are almost entirely 
non-ionized and they will absorb all radiation with a wave-length shorter than 
the Lyman limit (A=A*=911 A). Any free electrons will be supplied largely 
by carbon, magnesium, silicon and iron. Since these elements are not very 
abundant, the electron density 7, will be relatively low. Seaton (11) has recently 
found that 7,~0-003 cm. 

TABLE II 
he potential in volts at the positive end of a grain 
Volume 


Ratio of axes h/a 
fem") 


>°O032 ° 





60 

121 

»$2°¢ 

ov / 

601 16°S 
I 210°0 121° 36°0 
2 530 253°0 Fy ae 


id 





6 O10'O 601°0 1608-0 


Now it has been shown in Section 4 that there will be regions of high electro- 
static potential on the ice grains. ‘lhe free electrons will tend to settle there, 
thus reducing the potential yp) to x, say. Van de Hulst (5) has discussed the 
capture of free particles by crystalline grains at low temperatures. He has pointed 
out that the accommodation coefficient will be small on a very pure grain. We 
therefore assume that electrons may be captured by an ice grain only by radiative 
transitions; conversely we also assume that electrons may only be removed by 
photo-ionization. It is then permissible to use Saha’s formula, modified to 
allow for the dilution of radiation, to estimate the extent to which electrons have 
settled on a grain. (Van de Hulst, on the other hand, has considered grains 
which are very impure and which capture practically all atoms, ions and electrons 
on collision.) 

The nature of the radiation field to which the grains are exposed depends, of 
course, on the distribution of stars nearby. For our present purposes we shall 
assume that its intensity is that found in a black body at a temperature 
T =10' deg. K, diluted by a factor 6=10~-', and cut off at the Lyman limit. 

The number of places available for electrons on a given area of a grain will 
be defined as that number of electrons which must settle there to reduce the 
potential from xy» to zero. Now let x be the fraction of these places which are 
occupied by electrons and let y be the corresponding potential which results on 
the grain’s surface. Then the modified Saha formula states that, in the absence 
of the cut-off, 


I—x (27m)*'? 


— een ome (eT (5 
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Here P,=electron pressure, 7’ =temperature of the free electrons, 
m = electronic mass, n=€X; 
« = electronic charge, k = Boltzmann’s constant. 
h= Planck’s constant, 
With the values here assumed 
(1 —x)/x =10°%e-™*?, 


This formula must now be modified to allow for the cut-off. Following 
Pannekoek (12) we assume that the absorption cross-section for the detachment 
by radiation of an electron from a grain varies as 1/v*, when v, the frequency, 
exceeds /h. ‘The absence of radiation with a frequency beyond the Lyman 
limit (v = v*) reduces the rate at which electrons are removed by a factor 


| © WET dy 
-- ={—e— kT 
| ewkT dy 
7] 
where we have put hv* =»*. 

Thus (1—x)/x=109(e"*? —e-”*?), when 7 <7*, and so x ~o. 

When 7>7*, no radiation reaching a grain has photons of sufficient energy 
to detach an electron from it. On the other hand the Saha formula shows that 
when equilibrium has been reached practically no electrons can remain at a place 
where 7 <7*. 

We let x*=7*/e. ‘The simple situation now arises: (i) If yy<x* on a given 
grain, no negative charges will settle. (ii) If yg>x* on a given grain, then 
electrons can settle until the potential is reduced to x*. Any further electrons 
which settle now are removed again by radiation. In this case the fraction of 
places occupied is given, on the average, by x =(x»— x*)/xo- 

Similar reasoning may be used to find the electric field on the condensation 
nuclei. Suppose that these are spheres, of radius a, composed of some 
conducting material. Electrons will be removed by radiation until equilibrium 
is reached when the potential at the surface is x* volt. ‘The electrostatic field E 
is then x*/a volt per cm = x*/300a e€.s.u. 

If a=107 cm, E=4:6~x 10° e.s.u. 

If a=10°* cm, E=4-6 x 10° e.s.u. 

In either case there will be an electric field which is sufficiently strong to 
initiate the growth of a regular crystal. 

Dr H. C. van de Hulst has kindly pointed out to the author that Spitzer has 
made a different estimate of the probable charge on a metallic grain (13). In this 
it is postulated that electrons and ions are captured largely by collision and that 
their removal by radiation is a relatively tnefficient process. With an assumed 
temperature of 10* deg. K for the matter in interstellar space Spitzer finds that 
there should be a potential of about — 2 volt at the surface of a grain. The author 
has used the same theory but assumed the temperature to be 60 deg. K. The 
surface potential then becomes +0-15 volt. This has the same sign as the 
potential worked out by the method of the present paper, but is considerably 
smaller. For a grain with a radius of 10~* cm the electrostatic field at the surface 
becomes E=5 x 10° e.s.u. This is just sufficient to initiate the formation of a 
regular ice-crystal. If the condensation nucleus is any larger there would 
probably not be any regular growth. 
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6. The electrical state of grains of various sizes.—In the same way that electrons 
collect at the positive end of a grain, positive ions may be expected to collect at 
the negative end. Here one can no longer assume that radiative capture will be 
important and Saha’s formula cannot be applied. If, however, the total charge 
of the positive ions settled on a grain were less than the total charge of the electrons 
the grain would carry a net negative charge and the capture cross-section for 
positive ions would become larger. It will therefore be assumed that the total 
charge of the ions and electrons on the surface is zero, and that, if x volt is the 
potential on the positive end, the potential on the negative end will be — x volt. 

Some values of (1 — x), the fraction of places not occupied by electrons, and /, 
the dipole moment, are given for some typical grain sizes in Table III. 


Tas.e Ill 


Ratio of axes b/a 


Volume 
(cm*) . 0'032 O'O! 








1°O 1'O 1-O 
6xro"* ro" 10 
0°67 1'o ire) 
4xto™ so7** 10 
03 1°o 1'0 
réxs0°™" 6x te" <10 
O'l4 0°48 1°0 
$-4x 107" 29 x 107! <10 
0°07 0°27 o'8 
i4xX16°° 4°8>» 10 
ovll 0°38 
6-6 10° 2°3 x10 
0°05 018 
3x10~° I*l X10 
0°023 0-08 

4x10°° 4°8x 1077 


‘Table 111 can be used to predict the size which a grain will reach betore the 
yrowth with a regular structure ceases. At 60 deg. K the field intensity has to 
exceed 6 x 10% e.s.u. for regular growth, while the intensity near the surface ot 
a grain will be approximately 47p(1 —x) =7°2 « 10°(1 — x) e.s.u. when a proportion 
(1 —x) of the available space is unoccupied. ‘The condition is (I—x)>0o-or and 
according to the values in the table the regular growth will cease when V = 107-1? cm’, 
tor ba=orl, 

The ice grains here considered have the long ellipsoidal shape envisaged by 
van de Hulst (14), who considers grains for which b=3 x 10°° cm, of the order 
of the wave-length of blue light. If b/a =o-1 this value of b gives V =1-1 x 10-" cm’, 
and the grains whose formation has been described in this paper can grow to 
the size required by van de Hulst. 

It may also be shown that the growth of a regular crystal ceases when 
6 =3x10-5cm even for smaller values of the eccentricity, e.g. when b/a=0-32 
and ¢=0°9. 
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Summary 

The paper examines the theory which ascribes the periodic reversal of 
magnetic field of some stars to internal oscillations. It is shown that, to get 
sufficiently rapid oscillations, a large internal magnetic field need not be 
postulated ; mechanical forces unaided can always produce a sufficiently 
rapid oscillation. Indeed, only for a particular class of oscillations do the 
mechanical forces produce a sufficiently long period : in oscillations of this 
class the motion is shown to be mainly confined to the surface layers. Such 
oscillations are examined to see if they can explain the periodic reversal of the 
observed magnetic field. The results so far are negative, but do not cover all 
the possible oscillations. The problem of the maintenance of the oscillations, 
and other related problems, are briefly discussed. 





1. In 1947 H. W. Babcock found that the star HD 125 248, already known to 
be a spectrum variable, is also a magnetic variable. Babcock has recently given a 
detailed account of the properties of the star(x). It possesses a magnetic field which 
reverses its sign periodically ; the period of magnetic changes, which is also the 
period of spectrum changes, is about 9:3 days. The polar field H, inferred from 
the observations, assuming a dipole field, ranges from +7000 to —6000 gauss ; 
the variations in H, are not exactly sinusoidal, the negative maximum being much 
sharper than the positive. The radial velocity of the star shows periodic variations, 
which appear to be correlated with the magnetic changes. ‘The magnitude of the 
radial velocity fluctuations is difficult to assess (different elements give discordant 
values), but it is of order one or two km/s. 

A magnetic field in a star takes a time of the order of the age of the galaxy to 
decay, and can hardly be created and reversed in a few days. ‘he apparent 
reversal of the field must accordingly be due to motions which alter the field’s 
relation to the line of sight. One obvious suggestion is that rotation is carrying 
the field round so that first the positive pole, then the negative, is being presented 
to the observer. Babcock rejected this suggestion because the required rotational 
velocities, of order 11 km/s, would imply a far greater degree of line broadening 
than is observed, and because the corresponding negative maximum would be 
much less sharp than is actually observed. Instead, he favoured a suggestion 
that the fluctuations in the field are due to free oscillations of the star in the presence 
of its field. 

The suggestion that magnetic variation is due to oscillatory motions was first 
studied by M. Schwarzschild (2). Schwarzschild made a number of simplifying 
assumptions in order to render the problem tractable. He assumed a uniform 
liquid star, with a uniform internal magnetic field in the absence of the oscillations, 
and supposed that the restoring forces making the star oscillate are mainly electro- 
magnetic, the mechanical forces due to gravity being negligible in comparison. 


* Received in original form 1952 April 29. 
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He found that to get a g-day period an internal field of order 10° gauss was 
necessary. In interpreting his results he assumed the 10® gauss to represent 
the field near the centre of the actual magnetic variable, the field falling off to a 
value comparable with 10000 gauss near the surface. 

Schwarzschild’s analysis is open to several cbjections. Ferraro and Memory 
(3), making assumptions similar to Schwarzschild’s save that they took explicit 
account of the decrease of the field near the surface, showed that a field far in excess 
of 10° gauss at the centre would be required to get an oscillation with a 9-day period ; 
a field consistent with surface observations would give a period comparable with a 
year. They suggested that a 9-day period must be due to excitation by a 
mechanical oscillation of the star. A similar suggestion was followed up earlier 
in an incomplete investigation by Miss Gijellestad (4). Another objection to 
Schwarzschild’s theory, in the form in which he presented it, was that he did not 
indicate how his oscillations produced a periodic reversal in the field; indeed, 
his analysis appears instead to indicate a field which increases and diminishes 
periodically, but without reversing sign. 

In the present paper the magnetic oscillations of a star are considered further. 
It is shown that, far from a large internal magnetic field being necessary in order to 
produce a 9-day oscillation, the mechanical forces evoked by the oscillation are 
quite strong enough on their own to produce a sufficiently short period, without 
any appreciable assistance from forces of electromagnetic origin. Indeed, the 
mechanical forces normally produce too rapid an oscillation; only by considering 
a special class of oscillation is it possible to increase the period up to that observed. 
Oscillations of this class are examined to see if they can actually produce a reversal 
in the observed magnetic field. ‘The conclusions drawn from the examination 
are not wholly favourable to the oscillatory theory. 

2. Consider the oscillations of a gaseous star under the joint action of electro- 
magnetic and mechanical forces. ‘The star is supposed to be non-rotating; its 
magnetic field in the absence of oscillations is supposed to possess general 
similarities with a dipole field or the field of a uniformly magnetized sphere, in that 
the field has a diameter of the star as axis of symmetry, and the lines of force are 
in planes through this axis. No assumption is, however, made about the precise 
variation of the field inside the star, save that it is presumably several times as 
great near the centre as atthesurface. As usual in considering problems of motion 
in a large part of a star, the electrical conductivity of the material is supposed to be 
infinite ; also pressure changes during the oscillation are supposed to be adiabatic, 
so that dissipative effects are neglected. 

In the equilibrium state about which the oscillations are taking place, the 
mechanicai equation of equilibrium is 


Vp =p8+ jaH, (1) 


where p is the pressure, p the mass-density, g the vector acceleration of gravity, 
j the current-density and H the magnetic field. ‘The electromagnetic equation of 
equilibrium is 


VaH = 4nj. 
The gravitational equation is 
V.g= —4nGp. 
In small oscillations the small changes in pressure, density, magnetic field, etc., 
at a particular point due to the oscillation are denoted by 5p, dp, 5H, etc. ‘The 
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velocity of a given element of mass, and its displacement from its equilibrium 
position, are denoted by v and & respectively. As usual in problems of small 
oscillations, squares and products of the small quantities 5p, 5p,..., § v are 
neglected. 

If E is the electric intensity during the oscillation (measured in electromagnetic 
units), the total electric field on the moving material is E+vaH. Since the 
conductivity is assumed infinite 

E+vaH=o. 
But, from the equation of induction, 


05H 
— —— =VaE. 
ot 
Hence 
03H 
ot 
or, integrating with respect to the time, 
8H =Va(EaH). (4) 
This equation represents the fact that the changes in magnetic field are due to the 
displacement of the lines of force, which are ‘“‘ frozen”’ into the material. 
The remaining equations of oscillation are the following : 


=Va(vaH), 


(i) The equation of motion. If the period of oscillation is 27/o, this equation 


po*& = Vip — pdg — dpg — jadH — djaH. (5) 
(ii) ‘The perturbed equations corresponding to the equilibrium equations (2) 
and (3); these are 
VadH = 479j, (0) 
V .8g= —4nGbp. (7) 
(ii) ‘The equation of continuity, 
= + V.(pv) =o. 
This can be integrated to give 
bp + V.(pE) =o. (8) 
(iv) ‘he adiabatic equation. In this it is to be noted that magnetic 
‘‘ pressures”’ are not to be taken into consideration ; for any work done by these 
is done at the expense of the magnetic energy, not of the heat. ‘Thus the adiabatic 
relation takes the usual form p cc p’, or in small oscillations 
Ap/p = vApi/p, (9) 
where Ap, Ap are the changes in pressure and density following the motion of a 
mass-element. ‘I'hese are connected with dp and dp (the changes at a particular 
point) by 
Ap=85p+&.Vp, Ap=5p+&.Vp. (10) 
From (8) and (10), another expression for Ap is 
Ap = —pV.6&. (11) 
3. The orders of magnitude of the different terms on the right of equation (5) 
can readily be estimated. We can speak of a wave-length L characteristic of the 
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oscillation; roughly this means that, if ¢ is one of the variables characteristic of 
the oscillation (like &, 5H, dp, 5p, etc.) then V¢ is comparable with ¢/L. If the 
oscillation is one in which there are no internal nodal surfaces, L is comparable 
with R, the star’s radius ; if it is one with many nodal surfaces, L may be small 
compared with R. 

In equation (4) the gradient operator operates on both § and H. Operating 
on &, it gives a term of order &AH/L ; operating on H, it gives one of order 
EaH/R, which may be smaller, but is not larger. Hence if the symbol ~ is used 
to denote “‘ is of the same order as’’, in general 

dbH~E,H L~EH'L. 

Similarly from equations (8) and (11), in general 

bp~pl/L~Ap. 
From equations (9) and (10) it next follows that 

dp~pé L ~Ap. 
Finally, from equations (3) and (7), 

g~4r7GRp, bg~4rGL dp, 
and from equations (2) and (6) 
4rRji~H, 4nLdj~sH. 
Using these results, we see that, of the terms pdg, gdp on the right-hand side ot 


equation (5), the first is at most of the same order of magnitude as the second. 
Also 


Vip~op L ~pé i. 
and, using the equilibrium equation (1) and noting that the electromagnetic term 


jsH is unlikely to be as large as the gravity term, 


gx Ppa tie PE 
Opg Spe oR RL’ 


Hence the terms pdg, gdp are at most of the same order of magnitude as Vop. 
Again, of the magnetic terms 5jaH, jadH in the same equation, the second ts at 
most comparable with the first. 

A pure magnetic oscillation, in the sense studied by Schwarzschild, is one in 
which the magnetic terms on the right of (5) are more important than the mechanical 
terms. As has been seen, Vdp is normally the most important of the mechanical 
terms; the most important of the magnetic terms is 5jaH, which is such that 


. Hb5H~ €&H? 
sjaAH~ — ~—F- 
47r7L 4nL? 
Thus the ratio of the two terms is of order 
H2 
P qa’ 
i.e. is comparable with the ratio of gas pressure to the magnetic “ pressure”’ 
H?* 87 in the equilibrium state. It can readily be deduced trom the equilibrium 
equation (1) that the magnetic pressure is not large compared with the gas pressure ; 
indeed, no one expects that the magnetic pressure can be more than a small fraction 
of the gas pressure. It follows that pure magnetic oscillations, in Schwarzschild’s 
sense, do not in general take place. 
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The qualification ‘‘in general”’ is necessary because of a possible breakdown 
in the order-of-magnitude argument. If the motion is nearly non-divergent 
(V.& small) and horizontal (&. Vp, §. Vp small), the quantities dp, 5p are also small, 
and the mechanical terms on the right of (5) are much smaller than the argument 
given above would suggest. In particular, for torsional oscillations, in which 
twisting takes place about the magnetic axis, dp, dp and 5g are rigorously zero; 
thus torsional oscillations are pure magnetic oscillations. But for practically 
every other class of oscillation the mechanical terms on the right of (5) are at least 
comparable with the magnetic terms. 

It follows that the periods of normal oscillations are in general not very different 
from those in the absence of the magnetic field. ‘The order of magnitude of these 
periods is given, according to (5), by 


pore ~Vip~pé/L2~REV p/ L?-~pgRE/L? 


o2~gR/L?. 


‘That is, o? is in general at least of order g/R, or GM/R®, where M is the total mass 
of the star. ‘This accords with known results; for example, for second-order 
spherical harmonic oscillations of a liquid sphere o? =4GM/5R*, and tor two 
similar oscillations of a polytropic gas sphere of index 3 the author (§) has tound 
the values o? =8-7G MR and o? = 4-7GM/R®. The values of o? found in radial 
oscillations are similar, but perhaps a little smaller. 

For the Sun, GM _R* is about 4 x 107-7 sec *; taking o? as 2G MR’, this gives an 
oscillation of period about 2 hours. For Babcock’s star, the corresponding period 
is probably longer, since the radius may be more than twice that of the Sun; but 
the period is unlikely to be more than about half a day. A similar result is got by 
extrapolating the observed periods for Cepheids (assuming these to be pulsating 
stars). ‘The Cepheids are all more massive than Babcock’s star, and their periods 
decrease with decreasing mass. [Extrapolating to the mass of Babcock’s star 
(likely to be, say, three times the mass of the Sun) one gets a period a good deal less 
thanaday. ‘Thus the more fundamental free periods of Babcock’s star should be 
less than one-tenth of the period of the observed magnetic changes. 

4. ‘Two conclusions at once follow. ‘The first refers to the size of the internal 
field. Schwarzschild demanded a central field of order 10° gauss, since otherwise 
he could not get a sufficiently rapid oscillation; Ferraro and Memory showed that 
an even larger central field was necessary. However, the above analysis indicates 
that the mechanical forces normally need no assistance from the electromagnetic 
forces in producing a rapid oscillation. Hence there is no reason to postulate a 
central field as large as 10° gauss. ‘The central field which one would infer trom 
the surface observations is not likely to be much above 10° gauss, and may be smaller. 
Thus the periods of oscillation are likely to be determined by the mechanical 
forces; in comparison with these, the electromagnetic forces should be negligible 
throughout most of the star’s volume. 

This conclusion would, of course, be unfounded if the oscillations were torsional 
oscillations round the magnetic axis, since in such oscillations the restoring 
forces are wholly electromagnetic. However, it can readily be seen that torsional 
oscillations could not give observable magnetic variations. ‘The disturbance 
magnetic field 5H to which such oscillations lead is directed along circles round the 
magnetic axis; its line-integral round any such circle is 47 times the current 
through the circle. Since electric currents flow in closed circuits inside the star, 
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the total current through a circle at the surface of the star vanishes. Hence 5H 
vanishes at the star’s surface; that is, a torsional oscillation produces no 
observable effects on the surface field. 

The second conclusion refers to the possible types of oscillation. Not only 
do the mechanical forces need no assistance from electromagnetic forces in 
producing an oscillation as rapid as that observed; they are normally too strong. 
It was found above that the periods of the major free oscillations of Babcock’s 
star should be less than one-tenth of the observed period of the magnetic changes. 
‘lo get agreement with the observations, an oscillation must be found in which the 
etfect of the mechanical forces is much diluted. As noted earlier, this demands 
that the motion during the oscillation should be nearly horizontal and non- 
divergent. Oscillations corresponding to such motion were considered earlier 
by the author (§), ignoring electromagnetic effects; they were then called gravity 
oscillations. ‘The oscillations possess a number of nodal surfaces (spheres con- 
centric with the star) on which there is no radial motion. ‘The nodal surfaces 
being only a small fraction of the radius apart, the motion between two of them is 
mainly horizontal; the mechanical effects of gravity in the small vertical parts of 
the motion are much diluted by having to drive the horizontal motions as well as 
the vertical ones. ‘To get a free oscillation as long as that actually observed for 
Babcock’s star, the oscillation must approximate to one of these gravity oscillations. 

An appropriately long period is easily provided by this means; indeed, since 
the observed period is in the range-covered by high-order gravity oscillations, there 
are many free periods not very different from the observed period. ‘This in itself 
introduces a difficulty; why should a particular free period be excited when 
neighbouring free periods are not excited ? This question must be given an 
answer before the oscillation theory of magnetic variation can be regarded as satis- 
factory. For the present, however, the possibilities of the theory are simply 
being explored in a general manner, and this question will therefore not be studied 
closely. 

5. Consider now the general properties of the oscillation. For definiteness 
the star 1s considered (as by Schwarzschild) to be viewed along its magnetic axis, 
and to be oscillating in meridian planes through that axis. Forces of electro- 
magnetic origin are neglected compared with the mechanical forces; their neglect 
is not valid in the extreme outer layers of the star—a magnetic field of 5 000 gauss 
leads to a magnetic pressure H?/8z of about one atmosphere—but it is likely to be 
valid in layers not far below the surface. Thus magnetic effects during the 
oscillation are mainly kinematical, due to the lines of force being dragged about 
by the moving material. ‘The velocities observed by Babcock (of order 2 km/s) 
would, during an oscillation, transport matter through a distance comparable with 
the stellar radius (say through 5 x 10!°cm) to and fro along the star’s surface. 

It is sufficient to use (slightly modified) the methods used by the author in 
studying oscillations of a polytropic star (§). Assume the star to be a polytropic 
sphere of index n; we put 


P=Poy'*", P=Pey"s T=To%, (12) 

where P,, p, are central values, and 
ro” = (1 +2) po/47Gp,’. (13) 
‘Lhe oscillation is supposed to be one in which 5p, dp and the radial component 
£, of the displacement & are each the product of a surface harmonic of order s 
and a function of r. ‘The cases s=1 and s=2 are those which give the greatest 
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magnetic variations when integrated over the whole of the visible hemisphere ; in 
applying the results we shall take s = 2, corresponding to a surface motion which is 
everywhere poleward, or everywhere equatorward, according to the phase of the 
oscillation. 

The perturbation of gravity, 5g, is neglected; this is amply justified here, since 
{as will be shown) the oscillation has little effect save in the surface layers, whose 
mass is small. ‘Thus, approximating by neglecting the electromagnetic forces, and 
using (1), the radial component of equation (5) becomes 
Odp _ 5p dp 
po" &,= 3 r r+ (14) 


Also, substituting from (5) into (8), 


> 28 dp dp 
ao dp = —V*5p+V. (* 4) 
0 (,00p _ 20P dp s(s+1) . 
~ Por ( Or p 4) * eS i 


s(S+1). 


(po Ww" €.) + 2 Op. 


~ Pr 


Again, combining equations (9) and (10), 


d 
+ (6p + é.. $)- : (2 +€ 2). 


Equations (14), (15) and (16) are the equations of the oscillation. 
Non-dimensional quantities Z, ys, x, « are now introduced, defined by 


E,=16, Sp=Pe, Sp=Pex, a= 47 Gpgx/(I +n). (17) 
Using these definitions and (13), equations (14), (15) and (16) can be transformed 
to ub 
yrat= (140) xy’, (18) 


a oy, (19) 


vyx=e+U(1+2—- we (20) 
where y’ =dy/dx. ‘lhese equations are equivalent to equations (11), (12) and (13) 
of the paper already quoted (5). 

Using the value of o corresponding to a period of 9-3 days, and an estimated 
value of p,, the parameter « is found to be not much greater than 10°°. Because 
of its smallness, further approximations can be made. ‘The terms in y and % in 
equation (18) are at least comparable with the corresponding terms in equation (20) 
(Oxb/Ax is at least comparable with , because the oscillation has a large number of 
nodes, and the distance between successive nodal spheres is much less than that 
corresponding to one unit of x). On the other hand, the term in ¢ in (18) is far 
smaller than the corresponding term in (20); the term in (18) may accordingly be 
neglected. ‘Thus (18) becomes 

oy 


(I+n)y'x= = (21) 


0 
_? - ee (ay"x*l) + 


and, substituting for x from this equation into (20), 


> y 0 , 
(1+ n) y"*y" {= - nyt —(I—pn)y i, (22) 
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where 


vy 
/ 


—_———_———-. (23) 
ny—I-n 


fe 

Here j is a parameter which is positive in a stable region, where the temperature 
gradient is less than the adiabatic value, and negative in an unstable region. 

Substitute for y and € from (21) and (22) into (19). Then, after some 


reduction, 
0 (x*uf dap hs {0 (x ' a+ 1) (2 
—— 6a | Pa ny w t= i = rast i} ——?. a 
ais" \ a * ox \y f » 


, 0 (x* oe as : 
In this equation the term — { -; } on the right is clearly negligible compared with 
Ox \y , 

s(s+1)(a+1)/x. Again, in considering the magnetic effects of an oscillation we 
are concerned only with the behaviour near the surface. Because the oscillation 
is one with many nodes, x? and y’ do not vary greatly in a “ wave-length” of the 
oscillation near the surface, and so may be treated as constant in carrying out the 
differentiation 0/dx on the left of (24). Accordingly, we put 


x? = X%, y’ =(y)eud) 


where x = X corresponds to the surface r = R of the star; and we also approximate 
by putting 


<n 2(V')rax 


where > =X—wx. ‘Then equation (24) approximates to 


0 0 
as {us + 5, (8 ~w) = Aus, 


where 


Ss +Il( n+), , 
A ieee . sia (y’)e=x: 


6. ‘To illustrate the general properties of the oscillations, suppose first that 
j. is constant and positive. ‘Then the general solution of (26) can readily be 
expressed in terms of Bessel functions; it is 


Az\"2 
ys F((*) ) P,(c0s ) 00s (ot te), (28) 
fL 


where F(u) is a solution of the Bessel equation of order n, and P,(cos@) is the 
Legendre polynomial of order s in terms of the polar angle 6. If m is not an 
integer, F(u) is given by 

F(u) = AJ,(u) + BJ_,(u), 


where A and B are constants; if m is an integer, J_, must be replaced by J,. 
‘The ratio A/B of the constants in this expression is determined by the boundary 
conditions. If the analysis of the last section were valid right up to the surface, 
the solution would be such that B =o, this corresponding to ys =0, i.e. dp =0, at 
the surface. In actual fact, because electromagnetic forces must be taken into 
account in the surface layers, the boundary condition must be rather different. 
Suppose, say, that y= 3, n=3. ‘Then the parameter yp is 3, and A is 
about 2000. ‘The large number of nodes and loops of the oscillation follows 
from the appearance of the large factor A/u in (28); save in the immediate vicinity 
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of z =0, successive zeros of F((Az/)"*) correspond to values of (Az/)"” differing 
roughly by z. Thus nodal surfaces are close together, and become increasingly 
closer as the surface is approached. At z=1, corresponding to a depth ot 
2 x 10cm, the distance between nodal surfaces is only 3-5 x 10°cm—only 24 per 
cent of the star’s radius; the smallest zero of F((Az/)"*) corresponds to a depth 
of only a few thousand kilometres, but here equation (28) is not strictly applicable, 
since the neglect of the electromagnetic forces is not permissible. 

It can readily be verified that (28) corresponds to motion which is mainly 
horizontal. It can also be shown that the horizontal motion is largely confined 
to the surface layers. If & is the 6-component of the displacement, with our 
approximations equation (5) gives 

2. = dsp 
po :e™ roo 
Os 


y" x00 


nae ((Az\t?\ dP,(cosé) 
const. x 3~"/?} ((~) see cos (at +e). (29) 


€g=const. x 


For large values of its argument, F((Az/z)!*) is an oscillatory function whose 
amplitude varies roughly as z=". Hence the horizontal displacement has an 
amplitude varying roughly as 7 “+24; the energy per unit volume of the 
horizontal motion is proportional to p£,*, and so roughly to =". 

If uw is large (corresponding to a temperature gradient nearly equal to the 
adiabatic) the oscillation possesses fewer nodal surfaces, and may be rather less 
of a skin effect. Regular oscillations are not possible if is everywhere negative 
(temperature gradient exceeding the adiabatic), but they are possible if y is 
negative in the surface layers and positive lower down. In the latter case it can 
readily be deduced from (26) that the oscillation has at most one nodal surtace 
in the region in which y is negative; in fact, often it has no nodal surface in this 
region, and z "#, which is proportional to the horizontal displacement, increases 
with increasing depth in the region. Even in this case, however, the oscillation 
affects dominantly the surface layers unless the unstable layer in which p <0 
extends to great depths. 

The fact that the oscillation is normally confined mainly to surface layers 
introduces a further difficulty of the oscillation theory. In the surface layers 
the assumption of adiabatic changes in temperature is not so well justified as at 
greater depths; this is especially true for an oscillation with a large number of 
nodes. One effect of departures from adiabasy is to damp the oscillations. 
To maintain the oscillations, as they are seen to be maintained, some source of 
energy must be available. ‘This source of energy, like the dissipation, must be 
situated in the surface layers; if it were in the interior it would tend instead to 
stimulate faster oscillations, with fewer nodal surfaces, less completely confined 
to the surface layers. 

One is driven to seek the source of energy in some instability of the outer 
layers. Whether an unstable layer, such as a hydrogen convection zone, could 
be present in the surface layers of Babcock’s star (an A-type star) is a question 
requiring further investigation. But an oscillatory theory of magnetic variation 
certainly demands a source of energy like that supplied by an unstable layer. 

7. ‘The fundamental question is, however, whether an oscillation like that 
studied above can explain the observed magnetic reversal. Up to the time 
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when this paper was first submitted, I was firmly of the opinion that it could. 
My belief was based on the expectation that electromagnetic forces, though not 
negligible near the surface, were not likely to upset the picture given above, of an 
oscillation mainly confined to surface layers, and possessing several nodal surfaces 
near to the star’s surface. ‘The argument then ran as follows. 

The lines of force crossing the surface of the star must be substantially the 
same at all phases of its oscillation. ‘To explain the reversal of polarity these 
lines of force must be, where they cross the surface, directed predominantly away 
from the observer at one phase of the oscillation, and towards him a half-period 
later. ‘This is secured by the rapid change in the horizontal motion with 
increasing depth; the total horizontal displacement during an oscillation may be 
5 x 10°km at the surface, and it drops to zero at a depth amounting to, at most, a 
few tens of thousands of kilometres. ‘The material drags the lines of force with 
it; thus at the instants of extreme displacement the lines of force are inclined to 
the surface at a small angle. ‘The oscillation which has been considered (corre- 
sponding to s = 2 in equation (29)) is such that at one extreme of the displacement 
the lines of force near the surface are bent towards the poles; at the next, they are 
bent towards the equator. ‘This creates the appearance of a reversal of polarity 
to an observer viewing the field from a point on the magnetic axis. ‘lhe pattern 
of the lines of force at the extremes of displacement would, on this theory, be 
roughly as shown in Fig. I. 








An analytical justification of these statements can be advanced. ‘The change 
5H in the field, due to the oscillation, is given by equation (4), i.e. 
5H =Va(EaH). 
In this, the spatial variation of & is far faster than that of H; hence the equation 
approximates to 


SH =(H. V)E—H(V.£). 


Again, the motion is nearly non-divergent (YV.§ small) and & varies far more 
rapidly in the vertical direction than horizontally. Hence if H, is the vertical 
component of H, approximately 

5H =H,(0€/dr). (30) 
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That is, 5H is nearly horizontal, like §; and it is large compared with H,, since 
the surface value of & is large compared with the depth at which & first changes 
sign. ‘This means that the sign and magnitude of the field seen by an observer 
looking along the magnetic axis is mainly determined by 5H, not by the undis- 
turbed field H. ‘The observed field reverses sign nearly with 5H. 

8. However, although the above argument is plausible, I have come to 
realize that it is vitiated by an oversight. With an arbitrarily assigned distribution 
of values of &, the field 5H given by equation (4) cannot in general be fitted 
smoothly, at the surface of the star, to the potential field which must exist outside. 
The condition of continuity of normal magnetic flux implies that the radial 
component of 5H is continuous at the surface; but if this radial component is 
once assigned over the surface, the potential field outside the star is completely 
determinate, and the horizontal component of 5H just outside the surface cannot 
be adjusted to give continuity with the horizontal component just inside. ‘That 
is, if equation (4) is used, with an assigned distribution of §, to determine 5H 
inside the star, the field so found cannot be fitted on to the external field without 
a discontinuity at the surface. Such a discontinuity implies a surface current- 
sheet (a finite current flowing in a layer of infinitesimal thickness). It can 
readily be shown that the distribution of & indicated in Sections 5 and 6 does, in 
fact, demand such a current-sheet. 

A current-sheet cannot in practice exist; it would imply finite for at the 
surface, due to the flow of the current across the magnetic field. We must 
therefore conclude that near the surface the distribution of § departs from that 
considered in Sections 5 and 6 in such a direction as to ensure continuity of 5H 
at the surface. Instead of a discontinuous change in field at the surface, a 
continuous transition takes place near the surface. ‘The electromagnetic forces 
become important near the surface, and ensure that, instead of being suddenly 
bent round at the surface, the lines of force are progressively turned round, in a 
certain transition zone, into the direction desired for continuity of 6H. 

The transition zone is fairly thick. It begins at the depth at which the 
electromagnetic force 5jaH first becomes comparable with the inertia force 
po*—. ‘The current corresponding to the field 5H given by (30) has approximate 
magnitude H, a&, 

= Or 
and so 0jaH has a horizontal component 
H,? 0&5 
7 or?” 
By (29), this is comparable in order of magnitude with 
H? Xo 
47 4p37q" 
The ratio of this to the horizontal inertia force po*&, is 


H*X 


lomp Z1y"po" 


or, using equations (12), (17), (25) and (27), 
H? »# PS 
16mp" mw” s(s-+1)(n +1)" 
As before, we take A as comparable with 2000, and p as $; X?/s(s+1)(m+1) is 
about 2. Hence the expression (31) becornes comparable with unity when the 
38 


J 


31) 
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gas pressure p is 2-4 x 10° times the magnetic pressure H?/8z, i.e. at a pressure 
of a few million atmospheres. ‘The depth at which this occurs, though small, 
corresponds to a layer far below the visible layers. 

It follows that by the time the visible layers are reached the adjustment of 
the lines of force to the direction of the external field is nearly complete. ‘Thus 
the surface field cannot be nearly horizontal, whatever may be the state of affairs 
in deeper layers. ‘The main effect of an oscillation is to displace the lines of 
force; the modification of their direction is a secondary phenomenon, and cannot 
by itself produce a reversal of the observed field. 

Suppose, for example, that the undisturbed magnetic field is a dipole field 
near the surface, the radial component of field on the surface being given in terms 
of the polar angle 6 by H,=fcos@. Suppose that at the extremes of the 
oscillation this is replaced by 

H,=2fBcos*@ or H,=2fsin*6@cos8, 


these being fields derived from the first by displacing the lines of force bodily 
towards or away from the poles. ‘The displacement implied is about as extreme 
as can possibly occur; in the first case no lines of force are left at the equator, 
and in the second none are left at the poles. ‘The horizontal field at the surface 
is then 

H,=%4fsin@ 


in the undisturbed state, and 


H,=%,Bsin0(5cos*é+1) or H,=+fBsin0(7 — 15 cos*@) 


10 

at the extremes of oscillation. ‘hus, in the last case, when the lines of force 
have been swept bodily towards the equator, they are not also tilted down towards 
the equator; in fact, they are inclined away from it if cos?@>. Similarly 
in the phase when the lines of force are swept towards the pole they tend to tilt 
equatorward. ‘lhe mean field along the line of sight @ = 0, averaged over the visible 
hemisphere, is }f for the undisturbed field, and 48 and 38 at the extremes 
of oscillation, if no regard is paid to darkening at the limb; the corresponding 
expressions are 0°303/, 0°342 and 0-264), if a darkening law 0:55 + 0°45 cos @ is 
used. ‘Thus the oscillation makes the average field fluctuate about its undisturbed 
value, but not reverse its sign, the amplitude of the fluctuations being at most 
13 per cent of the undisturbed value. ‘This being the result for as extreme a case 
of oscillation as that considered above, it is clearly idle to expect an oscillation to 
produce a reversal in the observed field. 

g. The above discussion has been given at some length because I have 
previously given publicity to the mechanism suggested in Section 7, and now 
desire to indicate clearly the reasons which have led me to reject it. ‘The 
mechanism may appear to have been given undue importance; it refers only to 
one particular type of field (a dipole-like field with magnetic axis pointing towards 
the observer) and to one particular kind of oscillation (one in which motion takes 
place in planes through the magnetic axis). Nevertheless one important fact 
emerges, which is not limited to the particular model assumed. ‘This is that the 
reversal of the observed field can hardly be brought about by a modification of 
the orientation of the lines of force crossing the visible hemisphere, without a 
change in the actual lines of force crossing that hemisphere. Now an oscillation 
cannot, in the short times available, alter the field appreciably by lifting through 
the surface lines of force which previously ran completely below the surface ; 
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the lines of force are too firmly imbedded in the material (6). Hence the oscillation 
can produce the observed reversal of the field only by transporting patches of 
different polarity from the invisible hemisphere to the visible at different phases 
of the oscillation. 

Whether the observed velocity-changes and the reversal of field can be 
explained in this manner requires further investigation. However, they can be 
explained only if the orientation of the field, and the mode of oscillation, are 
both different from those envisaged earlier. ‘The axis of the undisturbed field 
cannot point towards the observer; it must be roughly transverse to the line of 
sight. ‘The oscillation must be one which displaces the “ poles” of the field in 
opposite senses; in the surface layers it may perhaps be rather like a torsional 
oscillation round a diameter perpendicular to both the line of sight and the axis 
of the field. Once again it is necessary to ask why such oscillations should be 
excited when the other oscillations of which the star is capable are not excited. 

To summarize, the difficulties of the oscillation theory of magnetic variation 
are as follows. First, a mechanism for maintenance of the oscillations remains 
to be found; second, one must explain why one oscillation should be preferred 
above neighbouring free oscillations; third, a satisfactory proof that oscillations 
are able to produce the observed reversal of field has still to be given. ‘The theory, 
moreover, gives no immediate explanation of the differences between lines of 
Fe, Cr, and the rare earths, differences which, Babcock suggests, imply that the 
lines of Cr and of the rare earths are emitted by different regions of the star’s 
surface. However, in spite of these difficulties the theory cannot lightly be set 
aside; the observed variations of radial velocity strongly suggest an oscillation. 
Clearly a far more detailed investigation is required than has yet been attempted. 

One last comment may be made. Some workers have suggested that the 
explanation of the discordance between values obtained for the Sun’s general 
magnetic field at different epochs is that the Sun is a magnetic variable, with 
period 23 years, say. If this is so, the magnetic variation can be due to an 
oscillation only if the effect of the mechanical forces is diluted far more even than 
in Babcock’s star. On the other hand, the motions accompanying the oscillation 
must be sufficiently violent first to bury and then to uncover the lines of force. 
One may ask if the very dilute mechanical forces are strong enough to produce 
such motions. One possibility which appears to be open is that the oscillation is 
a relaxation oscillation, due to some thermal mechanism. Direct evidence of the 
motions which should accompany such an oscillation would appear highly 
desirable. 

In conclusion I would like to express my appreciation for inspiring dis- 
cussions with Dr H. W. Babcock and Dr M. Schwarzschild last year, when 
enjoying the hospitality of the astronomical departments of the California 
Institute of ‘Technology and of Princeton University. 


University of Leeds: 
1952 July 19. 
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REMARKS ON HALE’S DETERMINATION OF THE GENERAL 
MAGNETIC FIELD OF THE SUN 


H. von Kliiber 
(Received 1952 July 4) 


Summary 


A short discussion is given of the Zeeman patterns of the Fraunhofer lines 
used by Hale in his determination of the general magnetic field of the Sun, 
the object being to detect a possible systematic influence on his measurements 
arising from the various forms of pattern. Although considerably more 
information is available to us than to Hale, no trace of such an effect could be 
found. 





As is well known, until now the only extensive published measurements of 
the general magnetic field of the Sun which appear to establish the existence of 
a dipole field of 25 to 50 gauss are those due to Hale (1, 2). Later measurements 
by various investigators (3, 4, 5, 6, 7) have not succeeded in confirming this 
field, in spite of the use of modern interferometric and electronic devices. 
Thiessen’s 1946 measurements (8) are an exception, but they are contradicted 
by nearly simultaneous Potsdam observations (4). It thus appears that if such 
a magnetic field exists, the magnetic axis of the Sun being assumed to be near 
the rotational axis, and if it is measurable by a study of Zeeman effects, its strength 
must lie at present close to the limits of observational accuracy, i.e. it cannot 
exceed about 2 gauss. 

The origin of the discordance between Hale’s measurements and those of 
the recent observers has not yet been explained. In addition we must note that 
Hale himself drew attention to a group of Fraunhofer lines which, in contradiction 
to the other lines he measured, did not show any measurable magnetic field. 

At the time when Hale completed his two fundamental papers (1, 2) in 1913 
and 1918, the theory of atomic spectra was in its infancy and very little had been 
discovered of the various complicated types of Zeeman pattern which are now 
known to exist. ‘These developments came after 1919 (9), so that in selecting 
his lines Hale had to rely entirely on empirical work on the Zeeman effect, done 
with sunspots or in the laboratory. At various times therefore it has been suggested 
that Hale’s results, and in particular the discordant behaviour of a number of his 
lines, may be due to his insufficient knowledge of the actual individual Zeeman 
patterns. In investigations nowadays the spectral lines are most carefully 
selected on the basis of the well-established theory of the Zeeman effect (11, 12). 
In what follows we have attempted a critical examination of Hale’s work with 
respect to his choice of lines, from this point of view. 

Today we know the electron configurations, and therefore the quantum 
numbers, for nearly all of Hale’s lines (10), and can deduce their Zeeman types 
and intensities immediately (11). All the Zeeman patterns occurring in his 
measurements (2, Tables III, IV, and X) are illustrated in Nos. 2 to 22 of our 
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Fig. 1, which follows a well-known convention and gives for each case the correct 
relative spacing and the relative intensities of the components. ‘The intensity 
scale is adjusted so that the sum of the 7 components always equals the sum of 
the o components. The Zeeman patterns for individual lines given in Figs. 2 
and 3 differ from their respective types in Fig. 1 only in their individual separation 
and intensities. Investigations during the past few years have differed from 
Hale’s in concentrating almost exclusively on lines having very large Landé 
factors and a simple Zeeman pattern. An example is 45250-218Fe with the 


oO , , ' 
pattern cs illustrated in Fig. I (13). 


Fic. 1. 


Fig. 1 shows that Hale, not entirely aware of the complexity of the Zeeman 
phenomenon, for the most part used lines with a complicated Zeeman pattern, 
namely those which seemed to show rather large and suitable Zeeman separations 
in sunspots and in laboratory spectra. In investigations of spot fields and also 
in those made up till now of the general magnetic field of the Sun, only the 
longitudinal effect, i.e. the o components, of the Zeeman pattern are really useful 
for measurement. But the fine structure of the components cannot be resolved 
even in a strong spot field (about 4000 gauss), still less in the weak general field 
(<50 gauss). Hale could therefore at best measure only the centres of gravity 
(intensity) of the o components. Because of the symmetry of the pattern, these 
centres of gravity can easily be determined from the distances and the intensities 
of the components. The intensities of the o components of a Zeeman pattern 
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for AJ =o are simply symmetrical around the pattern itself. For AJ = —1 the 
intensity ratios of the o components are given by 


(J +J)(J+J—1): (J+J —1)J+J—2):.... 


2J(2J —1) : (24/ —1)(2J —2):....12:6: 2. 
‘The distance of each component from its neighbour is always equal. Applying 
the principle of moments and taking the distance between the weakest and the 
strongest component as unity, one gets for the position of the centre of intensity, 
measured from the weakest component, 


I (ee) will 
=; 





2J —2 4(2d —1)2J(2J + 1) 


independent of the value of J. From the known g-values of the pattern the 
distance of this centre of intensity from the normal position of the line can easily 
be calculated. 
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In Fig. 1 these centres of gravity are marked by small vertical strokes above 
each typical pattern. The distance of the centre of gravity from the normal 
position of the line may be called A, by analogy with the usual Landé factor. 
This quantity is then connected with the field strength H and with the actually 
measured Zeeman displacement AA through the well-known equation 


AA=C.H.A*.A,. 


With these centres of gravity and with this law we can immediately recalculate 
the polar field strengths from Hale’s published measurements (i.e. from his A,. 
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in 2), on the basis of modern well-established Zeeman effect theory. Hale, on 
the other hand, was compelled to perform this reduction with the help of special 
and rather complicated laboratory experiments. If there are really systematic 
errors in Hale’s results due to the imperfect knowledge of Zeeman patterns 
existing at that time, this recalculation should disclose them. ‘To avoid the 
publication of extensive tabular material the results are shown graphically in 
Figs. 2 and 3, together with other data concerning the individual lines. 

The abscissae of Figs. 2 and 3 are the wave-lengths and identifications of 
Hale’s lines according to (10). J is the intensity of the line from (15); E.P. is 
the excitation potential of the lower level from (10); A, denotes the Zeeman 
splitting for the above-mentioned centres of gravity, in the usual Lorentz units. 
‘The number attached to each point indicates the Zeeman type according to 
Fig. 1. Lines having two points in the diagram are blends. As to the value of 
H,, the filled circles give Hale’s published polar field strengths for a general 
dipole field(2) ; the open circles show the results of the re-calculation just described. 
(For a few lines Hale gave no values of the polar field strength, and of course 
none for the separately shown lines of Fig. 3, which gave negative results.) 
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We notice at once how extremely well, as a rule, Hale’s empirically reduced 
H, values agree with those calculated from the theoretical Zeeman splitting. 
In view of the very complicated, tedious, and correspondingly less certain 
procedure which he was compelled to follow, this agreement is another proof 
of the extreme care with which his work was done. The only serious and not 
immediately explicable differences are exhibited by the lines AA 5929-700 Fe, 
607902 Fe and 6128-990 Ni, but they do not influence the final result to any 
appreciable extent. ‘The discordance in the last-mentioned line may arise because 
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it is of type No. 17, making it particularly unsuitable for measurements of this 
kind. ‘There is no indication whatever of any systematic difference depending 
on Zeeman type between Hale’s results and these new reductions. The mean 
H, from the 30 lines assembled in Fig. 2 comes out from Hale’s measurements 
to be 24:3 gauss and from the present re-calculation 24-1 gauss, practically the 
same value. It may be noticed that several of these lines giving positive results 
belong to the same multiplet and their agreement with each other is actually 
quite good, as one might expect. There seems also to be no obvious correlation 
with the line intensities, or excitation potentials, apart from a certain independence 
of the line intensity, which had already been noticed by Hale. His explanation 
postulated a decrease of the field with height, but there seems to be no support 
for this view at present. 

On the other hand, discussion of the A, of those lines which had been marked 
by Hale as showing no effect does not provide any real clue as to the reasons for 
their behaviour. Several of these lines even belong to multiplets of which other 
members show a positive effect, and their physical behaviour is in every other 
respect very similar to these others. These lines are 


no effect positive effect 
A5848-09 Fe®D, —7D, A5827°89 Fe*®D, —‘D, 
5997°808 Fe °F, —>*F, 5927°798 Fe *F, —°F, 
6081-421 V *Daj —*P ae 6039°690 V 4Ds9—*P se 
6119°305 V *Ds_q—'Pape 


For the lines A6012-251 Ni and A 6042-092 Fe the quantum numbers were not 


available and the Zeeman effect therefore could not be calculated. But from 
sunspot observations with quarter-wave plate and polaroid it can be seen clearly 
that this line has a very wide and clear Zeeman separation, its A, being probably 
about 2-0, or even more. It is also remarkable that the line A6302-507 Fe, which 
splits into an extremely large, simple, and clear pattern in spot spectra, is found 
among those lines giving no effect in Hale’s measurements (as pointed out already 
by Hale himself). Also careful examination of the lines from Figs. 2 and 3 with 
the aid of the Utrecht Photometric Atlas, the Mt. Wilson Sunspot Atlas, and also 
our Own sunspot observations, gave no explanation why the one group shows a 
positive effect and the other not. The mean of the values of A, from the lines 
of Fig. 3 (showing no effect in Hale’s measurements) is I-36; on the other hand, 
this mean for lines with positive effect (Fig. 2) is only slightly greater, being 1°53. 

Summarizing, we can say that the lack of a comprehensive knowledge of the 
exact pattern and intensity ratios of the Zeeman components at the time of 
Hale’s work, which led him sometimes to choose lines not the most suitable for 
his purpose, seems to have introduced no systematic error whatever into his 
results. Unfortunately, however, our re-discussion has been unable to disclose 
why a number of apparently entirely suitable lines used by Hale failed to show 
any effect in his measurements. 


The Observatories, 
Madingley Road, 
Cambridge : 
1952 June 30. 
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IN'TERFEROMETRIC MEASUREMENTS OF WAVE-LENGTHS 


1, DEVELOPMENT OF THE METHOD OF CIRCULAR CHANNELS AND ITS 
APPLICATION AT 5080A AND 6020A 


VU. G. Adam 
(Communicated by the Director of the University Observatory, Oxford) 


(Received 1952 July 22) 


Summary 


Collisional shifts are now known to contribute to the ditferences between 
solar and laboratory wave-lengths. In view of this, it is in practice essential 
to determine the absolute wave-lengths in each of the sources independently. 
Differences reaching as much as o-o10 A between the accepted earlier 
measures show that the absolute values of the solar wave-lengths, at least, 
have not yet reached the accuracy required for comparisons with theory. 

A new method of wave-length measurement has been developed, as 
suggested by Treanor in 1949. ‘This method of “‘ circular channels ’’ enables 
solar and laboratory wave-lengths, in any spectral region, to be determined 
by direct comparison with the primary standard. ‘The methods of observa- 
tion, measurement and reduction are described in detail. In common with 
the usual method of “‘ parabolic channels” particular attention has to be paid 
to the displacement of the centre of the interference pattern, the determination 
of the focal length of the camera lens, and the correction for phase change. 
all as functions of wave-length, as well as to departures from the “‘ square 
law”’. Measures have so far been made at 5080 A and 6020 A, for solar lines 
and corresponding lines in the vacuum arc. 

The “ red shift” of the solar lines is 0-0093 A at 5080 A, but only 0-0065 A 
at 6020 A. ‘These results cannot be explained as being due to a combination 
of the relativity shift (0-0107 A at so80 A and 0:0126 A at 6020 A), with any 
reasonable collisional shifts. 





Introduction.—Vhe absolute wave-lengths of the Fraunhofer lines are 
fundamental observational quantities for solar physics. The highest possible 
accuracy is required in the values if they are to be of use in investigating the 
present-day problems of the solar atmosphere. It is now more than twenty years 
since standard values were recommended for wave-lengths in the visible region 
ot the solar spectrum, and during that time the recommended values have been 
tacitly regarded as final. ‘There are obvious reasons why the standard solar 
wave-lengths should now be reconsidered. In the first place the recommended 
wave-lengths can only be regarded as tertiary standards, since in none of the 
three series of observations from which they were derived was there any direct 
comparison with the primary standard. Furthermore, the reliability of both 
the light sources used as secondary standards has since been questioned. Finally, 
if we examine the agreement between the three series of values, we find that while 
in the red there is good agreement, in other regions the extreme differences 
between the three observations can reach as much as 0-010 A. The existence of 
these discordances can leave no doubt that there must be some re-observation 
of absolute values before any satisfactory interpretation of solar wave-lengths 
can be reached. 
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In 1949 ‘Treanor (19) described a new interferometric set-up for solar 
spectroscopy in which the interferometer is placed behind the focal plane of the 
spectroscope. He pointed out the advantages of the arrangement for wave-length 
determinations and demonstrated the high precision which could be obtained for 
absorption lines. ‘The overwhelming advantage of the method of circular 
channels is that it makes possible a direct comparison of the primary wave-length 
standard with any region of the spectrum, combined with simultaneity of 
exposure for the Sun and standard source through the same optical train. In 
principle the solar wave-lengths so determined may themselves rank as secondary 
standards. ‘T'reanor’s suggestions for wave-length measurement have been 
developed in the manner now to be described, and the method has been applied 
to the determination of absolute wave-lengths in the Sun and in laboratory 
sources. ‘his departure from the standard experimental arrangement has 
necessitated a reconsideration of wave-length determination from interferometric 
patterns, and in the course of reduction several factors important in all 
interference methods have revealed themselves. Measures have so far been 
made in two regions of the spectrum, at 5080 A and 6020 A. Parts I and II of 
the paper give respectively the observational procedure and the methods adopted 
in measurement and reduction of the plates. The wave-lengths obtained are 
given in Part III together with a general discussion of methods for the 
measurement of solar wave-lengths. ‘The investigation is primarily concerned 
with the problem of the determination of absolute wave-lengths for solar lines, 
but a brief discussion is given in Part IV of the theoretical interpretation of the 
solar red-shift obtained in these first two regions. 

Part I, OBSERVATIONAL PROCEDURE 

1.1. Experimental arrangements.—The general arrangement of the apparatus 
is shown in Fig. 1 and Plate 5. ‘The spectroscope is used primarily as a 
monochromator, and the interferometer is placed in a parallel beam behind its 
focal plane. Each spectrum line is an off-axis slit source for the interferometer, 
and gives in the camera focal plane a chordal section (appropriate to its position) 
of the full pattern corresponding to its wave-length. Light from a standard source 


is simultaneously passed through the interferometer by means of the cube with 


a partially reflecting diagonal face. ‘The interference pattern of the cadmium 
lamp and that given by the light through the spectroscope can be photographed 
simultaneously without superposition, because of complementary diaphragms, 
in the shape of a cross, D, effectively on the field lens F, of the spectroscope, and 
D, on the field lens IF, of the standard lamp arm. (See Plate 6.) 

In Fig. 1 the camera lens was an 18 in. Cooke Series 1X f to apochromatic 
process lens, and the collimating lens C was a 23cm doublet. ‘The spectrum 
was thus imaged with a magnification of two in the focal plane of the camera, 
where an area approximately 2 in. square could be photographed. fF, and F, 
were identical 50 cm scratch-free spectacle lenses. F, carried the diaphragm D, 
on its surface, and the lens and diaphragm could each be adjusted so that the 
two diaphragms D, and D, were simultaneously in focus on the plate, and the 
images coincident. ‘lhe beam-splitting cube was formed by cementing together 
two right-angled prisms with a layer of zinc sulphide evaporated on to the diagonal 
face, giving about 10 per cent intensity in the reflected beam.* ‘The cube was 

* | am very grateful to Dr K. M. Greenland, Chief Physicist of the British Scientific Instrument 
Research Association, who made the cube. 
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carried on an adjustable mount inside the box, which also carried the lens F,. 
The reflecting-cube box and the interferometer box carrying the collimating 
lens C were connected by flexible metal bellows which could easily be removed 
and permitted focusing adjustments. The same type of removable bellows was 
used in the cadmium arm between F, and the filters. 
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Plate 


Fic. 1.—Arrangement of apparatus in the method of circular channels. 


In setting up the optical train the camera (ring lens) was first focused for 
infinity for red light by Hartmann tests using the parallelism of interfering 
beams crossing the interferometer. The interferometer was then removed 
from its box until the lining up of all the components was completed. The 
solar spectroscope was focused for each region to bring its focal plane into the 
plane of F,. The ring lens was then lined up with the combined camera and 
collimating lens of the spectroscope. ‘The interferometer box carrying C was 
next put at approximately the right distance from F, and aligned with the ring 
lens and spectroscope lens. The reflecting cube was then inserted, and turned 
and levelled so that it did not alter the alignment. Finally the exact position 
for the collimating lens C was found, from Hartmann tests for the focusing of 
the solar spectrum re-imaged in the camera focal plane. 

The interferometer plates were crystalline quartz and were silvered to 
approximately go per cent reflection in the green, and diaphragmed to a circular 
area 14 mm in diameter. The plates had angles of 15’ between the silvered and 
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Photographs of apparatus for determination of wave-lengths by the method of circular channels. 


M. G. Adam, Interferometric measurements of wave-lengths. I. 
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(a) Measuring machine for interference spectra. 

(b) Solar spectrum. 5080A region. 24 mm spacer. 

(c) Solar spectrum. 6020 A region. 24 mm spacer. 

(d) Manganese vacuum-arc. 6020 A region. 5 mm spacer, 
(e) Iron vacuum-arc. 6020A region. 24 mm spacer. 


(f) Over-exposure of cadmium pattern to show ghost patterns. 


[Note: (6), (c), (d), (e) and (f) are negative prints.| 
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M. G. Adam, Interferometric measurements of wave-lengths. |. 
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unsilvered surfaces. ‘The first plate gives a ghost pattern which is a part of the 
main interferometric pattern, and effectively increases the area of the source; 
the line of slope for this plate was kept perpendicular to the spectroscope 
dispersion. ‘The ghost pattern of the second plate is not a part of the main 
pattern, and the line of slope was turned to different orientations for a set of 
exposures, so that errors introduced by this ghost pattern had a random 
distribution. Plate 6f is a heavy over-exposure of the cadmium pattern to bring 
out the different ghost patterns. ‘The interferometer was mounted inside the 
light-tight box (Fig. 1); it could be turned and levelled to bring the centre of the 
ring pattern to the centre of the image of D,, and the plates could be made parallel 
in situ. No pressure or temperature control was attempted for the interferometer, 
though temperatures were read for each exposure, and barometric pressures 
were recorded during the day. All the measured wave-lengths depended upon 
simultaneous exposure of sources. ‘l'ables for the refractive index of air were 
drawn up from the results of Barrell and Sears for the dispersion of air (§), but 
over the range of pressures and temperatures used, the corrections at 5100 A 
and 6000 A for simultaneous exposures were negligible. 

The standard source was an “‘ Osira’’ cadmium lamp run at 1-21-45 amp a.c. 
Barrell and Puttock (4) have shown that the wave-length of the cadmium red line 
in the Osira lamp with currents between I and 1-5 amp is the same (to well within 
I part in 10’) as the wave-length in the accepted standard Michelson source. 
The red line was isolated by Wratten filter 27 in the position shown. No other 
line with an intensity as great as I per cent of the cadmium red line is transmitted 
by this filter, and there is no systematic disturbance of the effective wave-length 
of the standard line on this account. 

1.2. Lines selected for measurement.—Wave-lengths at the centre of the 
solar disk have so far been measured in two regions of the spectrum, 10 lines 
round 5080 A and 7 lines round 6020 A. ‘The shorter wave-length region was 
chosen because the standard solar wave-lengths of 1928 show outstanding 
discrepancies there ; the 6020 A region because it is close to the neon secondary 
standard lines and the cadmium primary standard, and the inter-agreement between 
the three observers is good. ‘The 6020 A region also gives a check on the 
wave-lengths measured with the Oxford equipment in 1948 (2). ‘The lines 
selected in both regions are shown in ‘Table I, together with the relevant portions 
of the 1.A.U. 1928 ‘lables (25), and the wave-lengths measured at Oxford in 1948. 
The Allegheny Observatory-Bureau of Standards wave-lengths measured by 
Burns and Meggers in column 7 were interferometric determinations using neon 
standards; Babcock’s measures in column 6 were made interferometrically using 
iron-are standards; and St. John’s values in column 5 were obtained with a plane 
grating using iron-arc standards. Columns 1, 2 and 3 in ‘lable I are from the 
Rowland Revision (24); all the solar lines are there marked “‘ m’’, indicating three 
or more independent measures. 

1.3. Details for solar plates.—The solar exposures were made at the centre of 
the 18 cm image in focus at the slit. ‘The out-of-focus image on a guiding disk 
about one metre in front of the slit was used for guiding during the exposures. 
In continuous-spectrum interferometry where the interferometer must be 
‘“‘crossed”’ with the auxiliary dispersion, ‘Treanor (19) has shown that there is 
an upper limit to the length of spacers which may be used without undue 
overlapping of orders. ‘This limit is set by the spectroscope slit pattern, and is 
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rarely greater than 10 mm for powerful solar spectroscopes. We shall see later 
that instrumentally it is the short spacers which are the most serious factor in 
limiting the accuracy obtainable for solar wave-lengths. ‘Treanor estimates 9 mm 
as the maximum spacer permissible with the Oxford spectroscope with a 
narrow slit. In the present work spacers of 2-5 mm and 5mm have been used; 
the best slit widths for non-overlapping of the heterochromatic channels and 
reasonable exposure times were found to be 0-10 to 0-12mm for the 2-5 mm 
spacer and 0-055 to 0-060 mm for the 5mm spacer. With these slit widths the 
exposures in the green, where the spectroscopic dispersion is 1-1 mm A~', were 
10-15 minutes for the 2-5 mm spacer, while an exposure of 30-35 minutes was 
barely sufficient for the 5mm spacer. Because of the long exposure times the 
wave-lengths at 5080A depend chiefly on the 2-5mm spacer. At 6020A 
(spectroscopic dispersion 0-62mmA~'!), the exposures were approximately 
4 minutes for the 2-5 mm spacer, and 15—18 minutes for the 5mm spacer. The 
wave-lengths at 6020 A depend equally on both spacers. 


TABLe I 


Measured lines 





Rev. Rowland Table LA.U. 1928 
Oxford 


1948 
Recommended | s+ John | Babcock |A.O.-B.S. 
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5068°773 Fe 
5072-082 Fe 
5072°679 Fe 
5°74°755 Fe 
5078-983 Fe 
5079°232 Fe 
5079°746 Fe 
5080°541 Ni 
5081-121 Ni 
5083 °347 Fe 
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6007°973 Fe 
6008 °572 Fe 
6013°503 Mn 
6016°653 Mn 
6021-808 Mn 
6024'073 Fe 
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The minimum-deviation link-work of the spectroscope was moved occasionally 
between exposures to bring the selected lines into different positions on the 
plate, and so differently placed with respect to the pattern centre. In the 
‘“‘normal”’ position all the lines appeared on the plate, in the “ shifted” positions 
either the long or the short wave-length end of the range moved off the plate. 

Astra VII or Ilford Rapid Process Panchromatic plates were used for the 
solar exposures, and all the development was in Kodak Press Contrast developer. 
The plates in the 5080 A region were obtained in 1950 June, July and August, and 
those at 6020A in 1951 July and August. 

1.4. Vacuum-arc plates.—For iron and manganese the lines selected in the 
Sun were also measured in a vacuum-are source. Unfortunately none of the 
selected lines is strong in the vacuum arc, and measurement proved to be a matter 
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of some difficulty. ‘lhe manganese lines can be excited with reasonable intensity, 
but they are very broad in both arc and Sun, and very different in character from 
the iron lines (except those marked “‘ nebulous”, such as 5074-747 in Table II). 
Burns and Walters have observed vacuum-arc wave-lengths (7, 8); their values 
tor the selected lines and their listed intensities are given in ‘Table II. 

The vacuum arc used has adjustable horizontal pole-pieces which during 
the exposures were kept not less than 10mm apart. Mild steel pole-pieces 
12mm in diameter were used for the iron lines. For the manganese lines various ‘ 
methods of excitation were tried. Manganese chloride previously fused and 
desiccated was placed inside a cored graphite pole-piece; the manganese lines 
were then bright, but their appearance was somewhat erratic and short-lived, 
and the discharge could only be maintained with the graphite pole-pieces close 
together. High manganese-content steel pole-pieces sputtered a great deal in 
the arc, and the manganese lines would have required long exposure times. ‘The 
method finally adopted was that recommended by Kiess and Meggers (13); 
the iron pole-pieces were replaced by copper ones, with metallic manganese 
fused into a cavity in the positive pole-piece. Occasionally graphite instead of 
copper was used as the negative pole-piece. ‘he manganese lines also appeared 
in the long-exposure iron-arc plates, and the line 6021-794 A has been measured 
in both sources. 

TasLe Il 


Vacuum-arc lines (Allegheny Observatory wave-lengths) 





Wave-length in A Int. Wave-length in A Int. 





Tron lines Tron lines 
5068-766 6007°959 
5072°077 6008 °558 
5074°747 6024'060 
5079°225 6027°048 
5°79°74! 
5083°338 Manganese lines 
6013°490 

6016°638 

6021°794 




















Direct current at 100 volts was supplied to the arc, and controlled at about 
6 amp. ‘The pressure was maintained at approximately 6 cm of mercury by means 
of a rotary vacuum oil pump. No doubt lower currents and pressures would 
theoretically have been desirable. ‘The currents and pole separations are, 
however, within the limits of the conditions given by Burns and Walters (8) for 
which no variations were found in the iron-arc wave-lengths, and pole effect is 
reported to be negligible up to 10cm pressure. Even with the conditions used, 
the exposures were 2}-4 hours for the iron arc with a 25mm spacer. ‘The 
exposures for the manganese arc were very much shorter than for the iron arc; 
with the conditions above, 10-20 minutes was sufficient. A manganese exposure 
was obtained at a pressure of the order of I mm, but there was no marked change 
in the appearance of the very broad manganese lines. 

The experimental set-up for the arc plates is shown in Fig. 2. The 10 mm arc 
was imaged by the lens L, at unit magnification over a hole 2-5 mm wide in the 
screen S, so that not more than the central quarter of the discharge was used. 
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The lens L, could be moved horizontally and vertically, and throughout the 
exposures, as the arc wandered, the image was guided to fall over the hole. The 
hole in S was imaged on to the spectroscope collimator-camera lens by the 
15cm lens L,; the extra-focal patch at the slit was just longer than the Ig mm slit 
height which was used for both Sun and arc. 

During the long iron-arc exposures it was necessary to stop the arc and 
destroy the vacuum several times to put in fresh pole-pieces. ‘The arc would 
normally run for } hour to 1 hour continuously; the time required for the 
iron-arc plates was thus considerably longer than the stated exposure times. 
‘The cadmium source was screened off during any interval for which the arc was 
temporarily stopped, in order to make the exposures strictly simultaneous. As 
the exposure estimated for the 5 mm spacer for the iron arc was 8 hours, only the 
2°5mm spacer was used, and the plates were all in the “ normal”’ position of the 
minimum deviation. For the manganese arc, with the much shorter exposures, 
equal numbers of 25mm and 5mm plates were obtained, and the lines were 
photographed in both the ‘normal’ and the “shifted” position. With the 
short times it was perfectly possible to obtain solar and are plates on the same day 
with the optical train entirely untouched. 





Fic. 2.—Experimental arrangement for arc spectra. 


‘There is no limitation on the spectroscope slit-width for emission spectra, 
such as there is for the continuous spectrum (Section 1.3). Short of overlapping 
of lines it is advantageous to use wide slits; this does not shorten the exposure 
times, but it gives a greater extent of the arc of the chordal section in the pattern, 
which is correspondingly easy to measure. Slit widths of 0-12 mm were used for 
the arc spectra in the 5080 A region. At 6020 A the width was increased to 0-2 and 
0-3 mm to facilitate the measurement of the very broad manganese lines. 

Ilford HP3 plates were used for all the arc exposures. ‘These plates are much 
faster than Astra VII or Ilford Rapid Process Panchromatic in both regions, 
but their low contrast makes them unsuitable for continuous-spectrum work. 


Part I]. D&TERMINATION OF WAVE-LENGTHS 

Treanor has already described his method of measuring the arcs using a 
travelling microscope and a revolving stage. ‘The plate is fixed on to the 
revolving stage with the pattern centre O as nearly as possible coincident with O’, 
the centre of rotation of the stage; let the distance OO’ =p (Fig. 3). The line of 
travel MM’ of the cross-wire intersection of the microscope should pass 
through O’. If it does so within a distance o, then (as ‘l'reanor shows) o may be 
neglected if it is less than about 0:03 mm. ‘The machine now in use (Plate 6a) 
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for the measurement of the plates has been designed to give the maximum 
stability over the rather long period required for the measurement of each plate. 
The most important adjustment is that OO’ should remain invariable, and to 
secure this as far as possible, the plate is clamped horizontally on a table mounted 
on the outer axis from a theodolite. The rotating table is traversed by a 
Cambridge travelling microscope, and the two are mounted together on a special 
casting. The microscope of the Cambridge instrument can be moved 
perpendicular to its line of travel, which makes it comparatively easy to secure 
that MM’ passes through O’ within the tolerable limit; estimates of o for typical. 
plate adjustments were all under about 0-005 mm. 





Fic. 3.—Plate measurement. 


As the table is turned, any arc 5 on the plate cuts MM’ in two positions, 
whose separation 27 is given by the microscope readings. We also read the 
angles at which these intersections take place, and so obtain ¢, the angle O'S 
makes with some fiduciary line on the plate through O’. For this fiduciary line 
we simply take the line through O’ coincident with MM’ when the table 
reading iso”. ‘The position of the pattern centre O must be determined. Suppose 
it is given by p and a, where « is the angular reading of OO’. ‘Then for each arc 
the triangle OO’S must be solved from r, ¢, p and « to give the diameter ) of 
the interference ring. 

The linear dispersion in the interference patterns is, at 5000 A, 1omm A ! 
for D=10mm (the lower limit approximately to the measured diameters), and 
22mm A’! for the largest diameters. ‘The microscope settings were recorded 
to 0-001 mm, each the mean of three settings, and the angles 4 were read to o 1. 

2.1. Cadmium patterns.—From the cadmium patterns we can obtain p and ~ 
independently of the theory of formation of the rings, except that they are 
assumed to be concentric circles. We measure straight across the arcs in arms 
a and ¢ at angles d, and 180° + ¢,, so obtaining 2r, and 2r, for each ring (Fig. 4). 
Similarly, we measure across arms b and d at ¢, and 180° + ¢g, to get 27, and 2ry. 
The differences 2r,—2r7, for each ring should be constant across a and c, and 

39 
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21%» — 2r%4 across b and d. Let 2Ar be the mean value of 27,—2r,, and 2Ar’ the 
mean value of 27,,—27r4, then 
_ cos ¢, —(Ar/Ar’) cos ¢, 

(Ar/Ar’) sind, — sind, 





tan @ 


ts 2Ar a 2Ar’ 
4cos(¢,;—«) 4c08(p.—«) 





p 


‘The true diameters of the cadmium rings are then given by 


D? =2r, . 27, + 4p” =2r, . 27g + 4p"; 


a mean value of the products 2r, . 27, and 27, . 2rg was used. For the cadmium 
patterns five to seven rings were measured on each of the four arms, as far as the 
plate allowed, over a range of ring diameters comparable with those used for the 
are and solar lines. 























Fic. 4.—Determination of pattern centre. 


The theory of the interferometric ring patterns is well known, but as we shall 
later need to reconsider the approximation involved in the usual treatment, 
it will be set down briefly here. Suppose that the optical thickness of the etalon 
is t, that the rings in the pattern are numbered outwards from the centre, and 
that the sth bright ring corresponds to rays which cross the etalon at angle 6, 
to the normal to the silvered surfaces (Fig. 5). For the innermost ring, where 
s=1, let the order of interference for light of wave-length A be p; then if ¢ is the 
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fractional order of interference at the centre of the pattern, we have the following 
relations : 
Centre: pte=2t/a. 
First ring: p=2t/A. cos 0. 
sth ring: p—s+1=2t/A. cos 4. 
For any two successive rings we get 
I =(cos 0,—cos 6,,,)2t/A 
=2t/A. A cos 6. 
i.e. A cos @ is a constant =A/2t, and for the (s+1)th ring we have, subtracting 
trom the “‘centre’’ equation, 
s+e=(1—cos @,,,)/A cos 0. (1) 
In reducing patterns D, is effectively measured on the plate, f, the focal length 
of the ring lens, is constant for the whole pattern, and tan 6,=D,/2f. For small 
angles we write tan 0,~0,, and cos 0#,=1-—6,?/2, and so obtain the well-known 
linear relation between squares of ring diameters: 

D? =D? +(s—1)AD?, (2) 
where AD? = 4f?A/t is the constant difference in the squares of any two successive 
diameters. ‘The fractional order « can be obtained from any ring from 

D?2/AD* =s—1 +e. (3) 
‘The agreement of the diameters determined in the cadmium patterns with the 
theoretical relation (2), and the constancy of the « values from the different rings 
using equation (3) give a valuable check on the performance of the camera lens. 








Fic. 5.—Theory of ring formation. 


When this test was made it appeared that the values of AD? determined from 
the outer rings were about 0-3 per cent larger than those from the inner rings and 
that the « values showed a small but systematic change with ring diameter. ‘The 
values of « from the inner and outer rings agreed, but there was a difference 
amounting to 0-0015A in wave-length for the intermediate rings. An exactly 
similar change in « with ring diameter had been found with a different camera 
lens which was in use in 1949, so that it was unlikely that the anomalies were due 
to the camera lens. ‘These departures from the “ square law”’ (equations (2) and 
(3)) arise from the approximation which has been made for cosine @ in equation 
(1). ‘The approximation cos 6=1— } tan® @ is only sufficiently accurate for @ less 
than about 1°°5 (cf. C. V. Jackson (12)). ‘The extreme arcs measured in this work 
correspond to = 2°-6, and the usual equations must be modified to take account 
of these larger angles. 


39* 
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We can obtain'a higher approximation for cos # most easily as follows : 
cos 6=cos (tan ' D/2f)=[1+(D/2f)?} 2 =1- }(Di2f)? 
+ #(D/2f)*—B(D/2f)>+ .... (4) 


In the usual procedure we take only the first two terms, and with the present 
values of 8 we can obtain ample accuracy by including the third term in the 
expansion. Remembering that for the chordal sections we get D? from the 
triangle OO’S (Fig. 3), it is simplest to adopt the following procedure. All the 
cosines in equation (1) are close to unity, and as (1) depends on differences of 
cosines we can work most easily in versine 6. Let 
versin @=1—cos 0= D*?/8f?. 

On the usual two-term approximation then D*=D, the measured quantity, 
but when we include the third term in (4) D*? = D?(1-D*/f?). It is a simple 
matter to calculate the slowly varying quantity (1 —%D?/f*) once and for all, 
for a given camera lens, asa function of D? (five figures are sufficient). We then 
multiply each observed D? by the appropriate factor to give D®**, the 
equations (2) and (3) then hold with D** used instead of D*, and the reduction 
proceeds exactly as in the well-known way. 

Since the peculiarities which had been found in AD?® and in « across the 
pattern were, within the errors of observation, due to neglect of the higher terms 
in equation (4), it was concluded that the camera lens was satisfactory from the 
point of view of field effects. Values of \D* and of « for the cadmium patterns 
were deduced from a least squares solution for the straight line (cf. equation (2)) 
D,* = D,*? +(s—1)AD®. 


2.2. Displacement of the pattern centre.-or the vacuum-arc and solar plates 
we measure 2r and ¢ for the ring sections of each line. ‘Three or four sections 
above and below the pattern centre were measured on the solar plates, and, 
in general, four or more on the arc lines. The principle of the reduction is that 
when p and « are known from the cadmium patterns we can convert values of 27 
into diameters from 


D? = 4p" + (2r)? — 4p . 2r cos (f—«). 


We next proceed to D** and calculate the fractional order «. From the 2¢ value 
derived from the cadmium « and the wave-length 6438-4696 A, we then calculate 
the arc and solar wave-lengths. When this was done the wave-lengths obtained 
showed differences from the Rowland Revision values, as was to be expected, 
but these differences varied quite rapidly across the spectrum, and also in 
comparing the ‘‘normal’’ and “shifted” positions, the wave-length of any one 
line depended upon the position on the plate at which it was measured. When a 
second set of plates, taken with a different orientation of the second plate of the 
interferometer, was reduced, it was found that the effects varied across the 
photographic plate in the opposite direction. ‘The anomalies must then be 
associated with the orientation of the second interferometer plate, and they were 
finally recognized as being due to the wedge angle between the silvered and 
unsilvered surfaces of the plate; the result is that the cadmium pattern centre 
we have deduced in the red will not be exactly the pattern centre for other 
wave-lengths. 

The inclination of the unsilvered surface of the second interferometer plate 
gives a turn to the cone of rays, forming the main pattern, in the direction of the 
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line of greatest slope. The differential turn for the red cadmium light and, say, 
the green light, will be (u,—j,)A, where A is the wedge angle, and p, and p, 
are the refractive indices of the plate material for the two colours; the centre of 
the green pattern is displaced from the centre of the red pattern by this angle. 
Using the values of the refractive index of quartz (ordinary ray), we find a 
differential deviation of 0°-00164 between the red and green normals (to the 
silvered surtaces), for a wedge angle of 15’._ With a camera focal length of 470 mm 
this corresponds to a shift of 0-0134mm on the plate. The direction of the 
displacement will be shown by the position of the second-plate ghost-pattern 
centre, which indicates the direction of the line of greatest slope for the second 
interferometer plate. For each orientation of the second plate an over-exposure of 
the cadmium pattern must be made to show this direction (Plate 6 f), and the true 
centre of the green pattern can then be deduced, say at angle £, distant p’ from 
()' (Fig. 3). It is p’ and £8, not p and «, which must be used in the equation for D?. 

A check on the centre shift as calculated above was made by taking cadmium 
patterns simultaneously in the light of 6438A and 5086A, using a special 
double-cross diaphragm as shown, with pieces of Wratten gelatine filter 25A 
and 61 to isolate the red and green lines. ‘The red and green centres could then be 


Fic. 6.—Double-cross diaphragm. 


Filled cross : Wratten filter 61. 
Hatched cross: Wratten filter 25A. 


located quite independently, and the mean displacement obtained for the green 
trom the red centre was 0-013mm, in directions within a few degrees of the 
roughly estimated ones. ‘his is the displacement applicable to the 5080A 
region; at 6020A the calculated displacement is only 0°0031mm. ‘The 
correction was applied at 6020 A, but it produces a maximum effect of only 0-001 A. 

2.3. Determination of camera focal length.—-\n deriving « from D** we need 
the quantity AD*(=4f*A/t), For an accuracy of o-oorA at 5000A, 
i.e. AA/A=2 x 10-7, we can see from the expression 


_, fr A(AD*) AD 
AA A= | AD? _ >| 


that for @=2 we need an accuracy of 1 in 10‘ on each of the terms in the bracket, 
i.e. 2 in 10* on AD*. Now f is not known with sufficient accuracy for this 
purpose, and \/)*? must be found from the interference patterns. For the 
cadmium patterns \))** can be obtained with the required accuracy, but the 
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uncertainty in the AD**® from the solar lines and the vacuum-arc lines is 
considerably greater than this limit. The mean value of AD**, using 
}(D,** — D,**) for the ten solar lines on a typical plate at 5080 A, has a standard 
deviation of the order 7 in 10‘, and for the arc plates (5080 A) the uncertainty is 
greater than this. It was therefore decided to deduce AD* on any plate from 
the well-determined cadmium pattern value. In doing this we cannot assume 
that the camera focal length is the same at all wave-lengths. A photographic 
lens is not normally corrected for the wide range in wave-lengths now in use, 
and we must consider possible chromatic variations in magnification. If we are 
working in the green at wave-length A, we need to know the ratio of the camera 
focal length f, at A, to the focal length f, in the red (A,), since AD*,,,.,: 
AD* 04 =Aghe : Af’. Now for the region 5080 A we can use the cadmium green 
line at 5086 A to determine f,?: f2 from AD,**: AD,**. This was first attempted 
using the 2-5 mm and the 5 mm spacer; the ratio was apparently well determined 
in each case, but the 2-5 mm spacer gave a value 0-9994 and the 5 mm spacer 0-g999 
for f,* : f,2. 

This variation in the measured ratio must be attributed to the unresolved 
multiplicity of the cadmium green line causing a slight change in the effective 
wave-length across the patterns. The determination was accordingly repeated 
using a 10 mm spacer for which the range is greater than the hyperfine structure 
pattern width (16), and the two satellites are sufficiently separated from the 
main line for the first seven or eight rings. Simultaneous exposures for the red 
and green patterns were taken with the double-cross diaphragm and a mean value 
for (AD*?/X) green: (AD**/A) oq Of 0°9986 was obtained from eight of these plates. 
The exposures were made for various orientations of the second interferometer 
plate, and of the double-cross diaphragm ; and using two pairs of interferometer 
plates, the second pair more heavily silvered than the usual go per cent reflection 
in the green. At 6020A no sufficiently monochromatic line was available from 
which to deduce fgo29: f,, so we have assumed that the ratio varies linearly with 
wave-length separation from the standard line, and taken a ratio f%go.9:f?? of 
©9997. Again, as the region is so much closer to the standard, this correction 
has a much smaller effect on the wave-lengths than in the green. 

2.4. Phase change.—The fractional orders we deduce must be corrected for 
the chromatic variation of phase change on reflection at the silvered surfaces. We 
can take this into account by adding a certain small correction k, to 2t, the optical 
path difference. The equation for light of wave-length A which crosses the inter- 
ferometer at angle @ to the normal then becomes nA = (2t + k,) cos 4, where n is the 
order of interference. For long spacers k, has but little effect upon the value of A 
deduced, and theoretically it can be eliminated by employing two spacers of different 
lengths, but this is only possible with sources of the highest spectroscopic quality. 
If hyperfine structure is present in a line, then the two different spacers may give 
different effective wave-lengths under the different resolving powers, and we cannot 
distinguish between the apparent changes in wave-length due to this cause, and 
those due to phase change. 

Since the phase change correction at 5080 A amounts to a wave-length 
correction of about 0°006 A for the 2°5 mm spacer, and 0:003A for the 5 mm 
spacer, it is essential to determine its precise value. Several sets of phase change 
plates were obtained and reduced using the cadmium red and green lines with the 
2°5and 5mm spacers. But these determinations gave anomalous results for the 
phase change, presumably because the effective wave-length of the green line 
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was different under the different resolutions of the two spacers. As no better 
laboratory source was then available, a value was deduced from the solar lines at 
5080 A using the 2-5 and 5mm spacers. ‘I'he solar lines are broad and not 
really suitable for this determination, but there is some hope that, being so 
broad, their effective wave-lengths will remain unchanged for the two spacers. 
The result obtained from nine lines for k,/A (relative to the cadmium red region), 
was 0-012 + 0:0026 to be subtracted from the observed fractional orders. Williams 
(22) gives a value 0-010 in this region. Barrell and Teasdale-Buckell’s results 
(6)* give 0-008, and a value deduced from Schultz’s recent experimental results 
(17) is 0-018. Lacking any more certain value we have used a phase change 
correction of 0-012 at 5080 A. 

At 6020 A the phase change correction amounts to a wave-length correction 
of only 0-001 A and less; a value for k,/A (relative to the cadmium red line 
region) of 0-002 has been adopted for the 6020 A region, from the results of 
earlier workers (22, 6, 17). 


Part III. OBSERVED WAVE-LENGTHS AND COMPARISONS 


3.1. Results for wave-lengths.—The wave-lengths obtained for the solar lines 
and the arc lines in the two regions are given in Table II] and Table IV. ‘The solar 
lines are corrected for diurnal rotation and the Earth's orbital motion. At 5080 A the 
wave-lengths are based onf,? : {2 = 0-9986 and a phase change correction of — 0-012. 
The “normal” (N) position of the lines is with the minimum deviation adjusted so 
that 5076 A is approximately central on the plate (as in Plate 6 6), for shifted position 
S I the plate centre corresponds to 5084 A, and for shifted position S II to 5070 A. 
At 6020 A the wave-lengths are based on f*g929 :/ =0°9997 and a phase change 
correction of —0-002. ‘The normal line position (N) for the 6020A region is 
with the plate centred at 6018 A (as in Plate 6c), and in the shifted position (S) 
it is centred at 6028 A. ~ 

The results obtained for all the lines from the individual plates are shown, to 
demonstrate possible systematic changes of the wave-length deduced with spacer, 
position on the plate at which a line is measured, or orientation of the second 
interferometer plate. ‘T’his last is shown by the angle given in ‘Tables III and 1V at 
the head of the columns. ‘The standard deviation in column 4 is the standard 
deviation of the mean in column 3 deduced from the scatter of the individual plate 
values. On the whole, the agreement under the different experimental conditions 
is satisfactory. At 6020A the solar wave-lengths with the 5mm spacer are 
slightly shorter than with the 2:5 mm spacer, but the effect is negligible in the arc 
spectra. If the difference really exists it may mean either that the phase change 
correction should be smaller than the 0-002 we have used, or that the lines give 
different effective wave-lengths under different resolutions. ‘The iron-arc lines at 
5080 A were very difficult to measure, owing to their intrinsic faintness and the low 
photographic density, even with the four hour-exposure. The value 5078-973 A 
on plate A1g3 is obtained from one ring only (measured above and below the 
centre). 

If we multiply the standard deviations in column 4 by +/n, where n is the number 
of values from which the mean is obtained, we obtain the standard deviation of a 
single measurement of wave-length o,. We can also arrive at an estimate of this 
quantity from the scatter of the « values from which a wave-length is derived, 


* I am very grateful to these authors who permitted me to see their paper before publication. 
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combined with the uncertainty in the « from the cadmium pattern, and it is of some 
interest to compare these two estimates. ‘I’o make this comparison typical solar and 
arc plates have been examined, and the results are given in Table V. o,in column 5 
is the standard deviation of a single determination of wave-length expressed in 
angstroms. It is a mean value for the plate deduced from the scatter of the « 
values for each line combined with the scatter among the ¢ values for the cadmium 
pattern. In calculating the tabulated values, the uncertainty introduced by the 
cadmium pattern proved to be negligible; expressed in angstroms the standard 


Taswe III 
Measured solar wave-lengths (individual plate values) 
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deviation of the mean ¢ for the cadmium pattern ranged from 0-0003 A to 0-0001 A. 
The quantity o, in column 6 is a mean value, 0-0019 for all the solar plates in 
Table IT] at 5080 A, 0-0017 A for all the solar plates at 6020 A, and so on. 

The differences between o, and o, are in the direction one would expect, i.e. 
the internal consistency of any one line is better than the inter-agreement for that 
line from plate to plate. This difference probably arises from slight random 
disturbances in the apparatus from plate to plate, but it might also indicate a real 
variability in the wave-lengths of the lines. There has been some suggestion that 
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the solar wave-lengths are variable (10), but this could only be assumed to be so for 
these observations, in so far as (o,—@,) for the solar lines is greater than the value 
for the arc lines. ‘l'able V shows some indication of such an effect, but it is of the 
order of 0-001 A and is probably not significant. 

The experience gained in the micrometer measurement and reduction, in 
this first application of the method of circular channels, suggests that certain 
precautions should be especially emphasized. The usual interferometric 
precautions are well known, and are all applicable. In particular the illumination 
over the field of the pattern should be as uniform as possible; certain small 


TasBle IV 


Measured vacuum-are wave-lengths (individual plate values) 
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systematic differences in the cadmium pattern arms a-c and b~—d appear to be 
dependent on the illumination of the field in cadmium light. For all spectra 
it is obviously desirable to keep as close as possible to the pattern centre, where the 
dispersion is high and the effect of most of the corrections is small. In emission 
spectra the position of the line on the plate does not affect the ease of measurement, 
but for solar spectra measurement is very much easier where the line is on the same 
side of the pattern centre as the centre of the heterochromatic fringes. ‘The 
arcs of the absorption lines are then formed on much broader heterochromatic 
channels, and cut these channels nearly at right angles. On the other side of the 
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pattern centre the arcs are formed on narrow channels of continuous background, 
and cut these channels obliquely. ‘This makes micrometer setting difficult, and 
may lead to systematic differences in the settings on the arcs above and below the 
pattern centre. The effect is comparable with the systematic differences between 
direct and reversed measures in the usual cross-wire measurement of spectra, and 
should disappear for the combined measures. ‘The fact that the solar wave-lengths 
are not dependent upon line position shows that it does disappear, but obviously 
better results can be obtained on the other side of the pattern. 

3-2. Comparison of results with earlier values.—T here are two factors which may 
cause systematic differences between the wave-lengths obtained by different 
observers. The first is the “scale error” for a region, which means that the wave- 
lengths as a whole in that region are not properly related to the primary standard. 
Scale error could be introduced in the present work, for example, if the phase 
corrections adopted are in error, and in tertiary standard wave-lengths it could 
arise from similar errors in the secondary standards. ‘The second factor is the 
spectroscopic resolution which is used in wave-length measurements. Williams 
and Middleton (23) stress the fact that it is not justifiable to try to determine the 


TABLE V 
Standard deviations of a single wave-length measurement 





Source and Region Plate Spacer No. of 7 
‘ » . Z Space 
No. P lines in A 
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A271 mm ‘ “0008 





























wave-length, except by detailed contour measurements, to an accuracy greater than 
about 4 of the half-width of a line. “‘ With practice, one can set on the same part 
of a fringe with remarkable consistency, but there is no evidence to show that it is 
the centre of gravity of the line ’”’ (23). 

It is well known that the solar lines are intrinsically broad; from the Utrecht 
atlas we can estimate, for an iron line of Rowland intensity 5 in the 5000 A region, a 
half-width of o-150A. If then the resolution in wave-length determination is 
comparable with the atlas spectroscopic resolution, we should have, on Williams 
and Middleton’s criterion, a wave-length uncertainty of 0-004A. In this 
connection we must note that in the I.A.U. 1928 Tables centre and integrated 
disk wave-lengths are combined together. The integrated disk lines will be 
considerably broadened on account of solar rotation, and their wave-length 
uncertainty must be correspondingly greater than that of the centre lines; in 
addition there are possible systematic differences between centre and integrated- 
disk wave-lengths due to limb effect (26). 

Table VI shows the differences between the mean values in Tables II! and 1V 
and the earlier measures given in Tables I and II, in the sense present value minus 
earlier value. 
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At 6020 A, where the earlier values are accordant, the solar wave-lengths are 
in fair agreement with those values, and the vacuum-are wave-lengths now 
measured agree, in the mean, with those of Burns and Walters. There is little to 
be gained from attempting to analyse the cause of the small difference in the solar 
wave-lengths from the earlier values ; it probably arises from the different resolving 
powers with which the measures have been made. 

At 5080 A we note that the solar values now obtained all fall within the range of 
the earlier discordant ones, closer to the wave-lengths of St. John and Babcock 
than to the Allegheny Observatory—Bureau of Standards values. The differences 
of the Allegheny Observatory (Burns and Walters) vacuum-are values from the 
present vacuum-are wave-lengths are in the same sense and of approximately the 
same magnitude as the solar differences. If the present absolute wave-lengths are 

Taste VI 


Comparison of wave-lengths 
Differences (unit 0-001 A) from values now obtained 
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correct, we must conclude that there is a scale error of about 0-005A in the 
A.O.-B.S. wave-lengths in this region, and of 0-002 A in the opposite direction in 
Babcock’s values. 

However, scale errors might well arise in the present absolute wave-lengths 
from the value adopted for the phase change correction (k,/A) and the camera 
focal length (f,). A source of the highest spectroscopic quality is needed for an 
accurate determination of these quantities, and it is proposed to use a mercury 198 
isotope lamp to make a careful determination of k, and f, over the visible spectrum. 
Once k, and f, are determined as functions of A, the values can be applied to any 
laboratory and solar wave-lengths measured with the same camera lens (f,) and 
the same interferometer plate coatings (k,). ‘he determination will give a check 
on the values now used, and is an essential subsidiary investigation in extending the 
wave-lengths to other regions of the solar spectrum. 
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3.3. Comparison of the method of circular channels with other methods.—The 
wave-lengths now obtained have all been determined by direct comparison with 
the primary standard, and in that sense the solar wave-lengths may for the first 
time be regarded as secondary standards. We have encountered several new 
factors which have had to be taken into account in the method of circular channels ; 
the question naturally arises to what extent these factors are also present in the 
older methods, and whether wave-lengths of comparable accuracy can be obtained 
by these methods. 

The experimental arrangement adopted in the earlier interferometric deter- 
minations of solar wave-lengths has been to “cross”’ the interferometer with a 
prism or grating spectrograph by placing the interferometer in front, using a ring 
lens or mirror to focus the interference pattern on to the slit. With a linear 
spectroscopic dispersion this gives the well-known parabolic channels in the 
continuous spectrum. ‘This experimental arrangement we will refer to as the 
method ot parabolic channels, in contrast to the present method of circular channels. 
For emission spectra the names are simply indicative of an experimental arrange- 
ment, and do not describe any appearance of the spectrum. 

We must first of all note that, using the same interferometer, the resolution 
obtainable is necessarily smaller in the method of parabolic channels than in the 
method of circular channels. In the method of parabolic channels when the ring 
pattern is first formed on the spectroscope slit and then re-imaged in the 
spectroscope focal plane, the resolution in the measured pattern is dependent 
upon the performance of the spectroscope as well as upon the interferometric 
resolving power. On Williams’ criterion above, if we assume the spectroscopic 
resolution up the slit is comparable with that parallel to the dispersion, then for 
large solar spectroscopes the instrumental pattern alone (1g) introduces a wave- 
length uncertainty of 0-oor A. The ring pattern measured in the spectroscope 
focal plane may be distorted by line curvature, and by variations in the spectroscope 
magnification up the height of the slit; the wave-lengths deduced will be 
affected by any chromatic variations in these effects. 

In the method of parabolic channels, as in the method of circular channels, we 
are limited by the spectroscope slit pattern to short spacers for work in the 
continuous spectrum. ‘The use of short spacers sets a limit to the resolution ; 
it means that the phase change correction is of importance, and also involves the 
use of wide-angle interference patterns. A short spacer gives a few widely 
separated fringes, and to find the three or four that we need, we are forced to use 
large values of 6, which means decreasing dispersion in the pattern, and departures 
from the square law. 

The highest available auxiliary dispersion must be used in solar spectroscopy, 
which means that witha prism spectroscope, using the method of parabolic channels, 
we are forced, in general, to measure the solar wave-lengths against secondary 
standards. With a grating spectroscope the primary standard can be photographed 
simultaneously with other spectral regions by the use of subsidiary cameras, 
though this means doubling many of the sources of error. ‘The recommended 
solar wave-lengths in Table I were all measured against secondary standards, and 
the accuracy of the values obtained depends upon the reliability of these standards. 
Burns and Meggers used neon standards, and St. John and Babcock the iron arc in 
air. Subsequently C. V. Jackson (11, 12, 27) found that the neon standards, on 
account of the presence of the isotope atom, were unreliable for small spacers, 
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since it was difficult to apply the corrections necessitated by the hyperfine structure 
to observations made with etalons of 5mm to 10mm separations. Williams and 
Middleton (23) have abandoned an attempt to determine accurate wave-lengths tor 
iron in the standard Pfund air arc because of the breadth of the lines, which is of 
the order of a hundred times greater than the accuracy of setting. 

It is in the setting up of the secondary standards that we encounter, in the 
method of parabolic channels, the problem of f, and k,, and the shift of the pattern 
centre with wave-length. The scale of the measured pattern depends upon the 
magnification of the spectroscope m,, as well as upon the focal length of the ring 
lens, since now AD? = 4Af,?m?/t. As in the method of circular channels we avoid a 
determination of f, and m, by determining AD?, but again this, and the deter- 
mination of k,, can only be done with sufficient accuracy with highly mono- 
chromatic sources. The displacement of the pattern centre by the wedge angle of 
the second interferometer plate also enters into the method of parabolic channels. 
The spectroscope slit can only take a diametral section of the ring pattern for one 
wave-length, the other patterns must be slightly displaced on the slit, and the section 
will necessarily be a chordal one. With wedge angles of only a few minutes the 
effect is probably negligible, but it must be considered if larger wedge angles are 
used, or the observed region is far removed from the wave-length of the primary 
standard. 

We may then conclude that the method of circular channels provides a method 
of wave-length measurement against the primary standard which is more direct 
and, because of the higher resolution obtainable, capable of higher accuracy than 
could be achieved in the method of parabolic channels. ‘The difficulties we have 
encountered are not peculiar to the method, but are also inherent in the method of 


parabolic channels. The origin of the discordances between observers in the 
earlier measurements may well be found among the sources of error we have 
discussed. 


Part IV. ‘THE SOLAR RED SHIFT 

The red shift of the solar lines with respect to the vacuum-arc lines is given in 
Table VII. ‘The shifts in column 4 are based on four-figure values of the means, 
which in Tables III and IV are given to only three figures. It is important to 
remember that the red shift is much more certainly determined than the absolute 
wave-lengths for the Sun and arc. Uncertainty in phase correction and in 
chromatic variation of magnification, for example, affects wave-lengths from both 
sources in precisely the same way, and changes in the values adopted would not 
alter the red shift. 

The properties of the emitting atoms likely to be relevant to the interpretation 
of the red shift are also tabulated. The multiplet data are from Miss Moore’s 
1945 Table(15). ‘The equivalent widths are those of Allen (3) and the pressure 
classes are from the Rowland Revision (24). ‘These classes, first assigned before 
the days of modern atomic theory, are based on the appearance and behaviour 
of lines under pressure and under different conditions of excitation. Whatever the 
interpretation of this behaviour, the classes give a useful indication of the character 
of the lines. Lines of class @ are low energy-level lines, always sharp and 
symmetrical, and showing small displacements under pressure. Class bare medium 
energy-level lines, not quite so unvarying as a, but both a and 4 are regarded as 
suitable for secondary wave-length standards. Classes d and e are higher energy- 
level lines, regarded as unsuitable for standards of any order; they are described 
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as asymmetrical, d towards the red, and e towards the violet. Under pressure and 
at the pole of the arc they show displacements in the same direction as the 
asymmetries. 
Taste VII 
The solar red shift 


(Solar wave-length minus vacuum-arc wave-length) 
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* Marked as probably blended in the Sun, in Mult. Table 1933 (14). 
+ A faint iron line is predicted in the Sun at 6021-82 A. 


Within each of the two regions the red shift is remarkably constant, especially 
when we note that 50790 A is doubtful because of blending in the Sun, and that 
line 6021-8 A may be affected by different amounts of blending with the faint iron 
line in the vacuum arc and in the Sun. Also, within each region, the red shift 
does not appear to be dependent upon excitation level, parity of the lower state 
(an effect considered by Carter (g)), or line strength, though within the groups the 
range in strength is small. Nor, in the 6020 A region, is the shift dependent upon 
the element, or even upon the unresolved multiplicity of the manganese lines. 
‘The manganese lines as now observed appeared single in both Sun and arc, but 
in emission, under different conditions of excitation from those now used, 6013 A 
and 6016 A have been found by other observers to be complex (20, 21), and Abt (x) 
has deduced the presence of hyperfine structure in other manganese lines in the 
Sun. 

In fact the only regularity one can note is that in both regions the largest shift 
is given by a line of classe. The larger shift might be due to a slight violet shift in 
the laboratory wave-length, which is in accord with the description of the class, 
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but we might equally expect a red shift in the laboratory for the class d lines, 
whereas the solar red-shift now observed for these lines is not peculiar. ‘The 
most striking result however is the difference between the two regions, since the 
red shift is systematically larger in the green than in the orange region. ‘The mean 
shift for the six lines at 5080 A is 0-0093A and for the six lines at 6020A only 
0°0065 A. Onthe basis of Table VII this difference in shift can be considered to be 
correlated either with wave-length or with the equivalent widths of the measured 
lines. 

The results in Table VII are obviously insufficient for a decision on the 
question of the interpretation of solar wave-lengths, but we can nevertheless 
consider to what extent they can be reconciled with existing theories. In the light 
of present-day knowledge of atomic and solar physics, we must consider three 
mechanisms which might cause shifts in the solar lines relative to the laboratory 
lines : they are the gravitational potential, perturbations in atomic energy levels 
by collisions, and Doppler effects. In a recent discussion of the interpretation of 
solar wave-lengths Spitzer (18) has surveyed our knowledge of collisional shifts, 
and we will summarize his conclusions. ‘The dominant collisional shifts in solar 
lines must arise from hydrogen-atom collisions. We have no laboratory evidence 
for the shifts thereby produced, but it seems theoretically reasonable to expect a 
red shift for resonance lines, and probably a violet shift for an absorption line with 
an upper state near to the ionization limit. For non-resonance lines between 
intermediate energy levels the sign of the shift is not known, but the order of magni- 
tude of the shift is 20-80 per cent of the collisional half-width of the line. From 
indirect evidence Spitzer assumes that the shift for intermediate lines will in fact 
be to the violet and then finds that the “ solar wave-lengths seem consistent with a 
complete absence of Doppler shifts on the solar surface’. He tentatively concludes 
that “‘the measured shifts are simply the sum of the relativity red-shift and a 
collisional violet shift”. 

On Spitzer’s interpretation we can see from ‘Table VII that at 5080 A a mean 
violet shift of 0-o014 A due to collisions would be required, but at 6020 A a mean 
violet shift of 0-o061 A (or 0-0056 A for the iron lines only at 6020 A). ‘The iron 
lines in both regions however arise from transitions between intermediate 
atomic energy levels. Moreover the two lines 5079°7 A and 5083-3 A, arising from 
a lower state very near to the ground state, require shifts of very much the same 
amount as those arising from states two volts higher. We may try to avoid these 
anomalies by re-introducing Doppler shifts, though it has been shown (2) that only 
very small shifts, corresponding to velocities say not greater than 0-2 km sec, can 
be reconciled with existing theories of the solar atmosphere. We must take the 
same streaming velocity for both spectral regions. If for example we suppose there 
is an upward streaming velocity of 0-2 kmsec™', then assuming the relativity shift is 
present, the collisional shifts required are a red shift of 0-0020 A at 5080 A, but a 
violet shift of 0-o021 A at 6020A._ If, on the other hand, we assume there is a 
downward streaming velocity of 0-2 km sec™, the collisional shifts required would 
be a violet shift of 0-0048A at 5080A and a violet shift of o-o101 A at 6020 A, 
These calculations are based on the mean observed red shifts in each region, 
00093 A at 5080 A and 0-0065 A at 6020 A. 

Neither of these modifications lessens the difficulties in Spitzer’s interpretation. 
Indeed, on the basis of these observations, with so many uncertainties in the theory 
of collisional shifts, we cannot exclude the simplest interpretation of solar red 
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shifts, namely that they are the result of collisional shifts alone. On this hypothesis 
the effect of hydrogen collisions is to produce a red shift which is a function either 
of wave-length or of intensity. Such an interpretation leaves out of account the 
existence of “limb effect”’, but as there is some indication (2) that limb effect is 
much more confined to the extreme limb than was formerly supposed, it may 
legitimately be regarded as a separate problem. 

When the displacements of the solar lines were attributed simply to a combin- 
ation of the relativity shift and radial currents, it was justifiable to consider only 
differences between laboratory and solar wave-lengths, because neither of these 
effects is operative in the laboratory. Now that collisional shifts must be included 
in the interpretation of the observed displacements, this procedure no longer 
suffices. The laboratory lines as well as the solar lines will suffer collisional 
shifts, and the factors producing the shifts are different in the two sources. As the 
theory for both sources is incompletely understood, the most direct procedure 
is to determine absolute wave-lengths in each of them independently. From 
laboratory absolute wave-lengths we can test the theory of perturbations of atomic 
energy levels under laboratory conditions, and deduce the ideally undisturbed 
energy levels. From the solar wave-lengths we shall then be able to determine 
the displacement of a Fraunhofer line relative to its undisturbed position, and this 
is the displacement that has to be interpreted in terms of the physical conditions 
of the solar atmosphere. In this approach to the problem, the first step is the 
determination of absolute wave-lengths for laboratory and solar lines, with an 
accuracy greater than has apparently so far been achieved in most spectral regions. 
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MEAN AREAS AND HELIOGRAPHIC LATITUDES OF 
SUNSPOTS IN THE YEAR 1947 
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The following results are in continuation of those given in M.N., 111, 422, 
1951, and are derived from the measurement at Herstmonceux of photographs 
taken at the Royal Observatories of Greenwich and the Cape and at the 
Kodaikanal Observatory, India. 

Table I gives the mean daily areas of umbrae, whole spots and faculae for 
each synodic rotation of the Sun included in the year 1947; the means for each 
year as a whole are included in Table II, which summarizes the yearly values 
since the last sunspot minimum of 1944. 


TABLE I 
Mean Daily Areas 


A 





ry _ 


Projected* Corrected for Foreshorteningt 
No. Rotation Days peg Se 


of Commenced Photo- ‘Cideee Whole pi a Whole 
Rotation U.T. graphed Spots Spots 
1946-47 


1248 Dec. 25°18 28 530 3270 2522 391 2438 3000 
1249 Jan. 21°52 27 602 3959 2113 415 2768 2605 
1250 ~=6©Feb. 17°86 27 879 5470 2089 621 4048 2574 
1251 Mar. 17°19 28 783 5392 1841 568 4080 2317 
1252 Apr. 13°48 27 507 2917 2172 378 2212 2759 
1253 May 10°72 27 756 4699 2743 579 3551 3301 
1254 June 6°93 27 385 2430 2440 292 1881 3048 
1255 July 4:13 27 540 3428 2780 397 2595 3482 
1256 July 31°34 828 643 4468 = 2743 497 3238 «= 3485 
1257 Aug. 27°57 27 429 2756 2575 313. -2030«= 3290 
1258 Sept. 23°83 27 529 3175 2233 375 2300 2818 
1259 Oct. 21°12 27 386 2444 1822 285 1805 2274 
1260 Nov. 17°42 28 387 2446 2303 273 1763 2815 

* Expressed in millionths of the Sun’s disk. 

t Expressed in millionths of the Sun’s hemisphere. 





= 


Faculae 


Table III gives for each rotation in the year 1947 the mean daily area of 
the whole spots (corrected for foreshortening) and the mean _ heliographic 
latitude of the spotted areas for both northern and southern hemispheres. 

The mean heliographic latitude of the entire spotted area and the mean 
distance from the equator of all spots are also tabulated. The mean values 
for the year 1947 are included in ‘Table IV. 

Tables II and IV are in continuation of similar tables in Monthly Notices: 
for the years 1874 to 1888: 49, 381, 1889; 1889 to 1902: 63, 465, 1903; I90I to 
1914: 76, 402, 1916; 1913 to 1924: 85, 1007, 1925; 1923 to 1933: 94, si 19343 
and 1933 to 1945: 110, 501, 1950. 
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The rotations in Tables I and III are numbered in continuation of 
Carrington’s series. 

The chief features of the record are as follows :— 

(1) The continued increase of sunspot activity culminating in an exceptionally 
high maximum unrivalled in the Greenwich records. As compared with the 
year 1946, the increase in mean daily sunspot area (corrected for foreshortening) 
was 45 per cent and that for faculae 32 per cent. 


TABLE II 
Mean Daily Areas 


ma —_ 


No. of Days Projected* Corrected for Foreshorteningt 


~ a ~ icerietanieamenpith 


Photo- Without Whole . Whole 
graphed Spots Umbrae Spots Faculae Umbrae Spots 
366 157 30 160 281 23 126 338 
365 14 102 560 774 78 429 940 
365 ° 389 2381 1793 291 1817 2188 
365 ° 558 3559 2326 405 2637 2894 
* Expressed in millionths of the Sun’s disk. 
+ Expressed in millionths of the Sun’s hemisphere. 





—, 


Faculae 


TaBLe III 


Spots North Spots South 
of the Equator of the Equator Mean 
rane ee Distance 
a Mean Latitude from 
No. Rotation Dail Dail Helio- of Entire Equator 
aily : aily ‘ : 
of Commenced tots graphic pice graphic Spotted of all 
Rotation U.T. Latitude Latitude Area Spots 
1946-47 - 3 
1248 Dec. 25°18 1195 19°85 1243 14°77 + 2'20 17°26 
1249 Jan. 21°52 354 15°06 24131767 — 13°48 17°33 
1250 Feb. 17°86 g12 16°62 3136 19°49 —11°36 18°85 
1251 Mar. 17°19 356 15°36 3724 21°62 —18+3g 21°07 
1252 Apr. 13°48 788 17°69 1424 15°89 — 3°93 16°53 
1253 May 10°72 1875 18°15 1676 19°44 + O41 18-76 
1254 June 6°93 1320 17°45 561 17°52 7°02 17°47 
1255 July 4°13 1684 15°66 ‘gII 20°38 + 3°O% 17°32 
1256 July 31°34 1385 14°72 1853 15°85 2°77 15°37 
1257 Aug. 27°57 1054 =: 15°14 977 13°54 me 14°37 
1258 Sept. 23°83 1367 18:21 933 14°37 + 4°99 16°65 
1259 Oct. 21°12 721 13°97 1084 14'27 2°99 14°15 
1260 Nov. 17°42 427 15‘10 1336 16°17 — 8-60 15°92 


For two rotations, viz. Nos, 1250 and 1251 (see ‘Table I) the mean daily 
area (corrected for foreshortening) exceeded 4000 millionths of the Sun’s 
hemisphere. ‘The largest spot group ever photographed (perhaps the largest 
in sunspot records) occurred in rotation No. 1251. ‘This group had a mean 
area during disk passage of 5520 millionths of the Sun’s hemisphere and a 
peak area of 6132 units. In the previous rotation, this group had a mean area 
of 3637 and a peak area of 4554 units and ranks as the fifth largest group in the 
Greenwich records. 





Mean areas and heliographic latitudes Vol. 112 
Taste IV 


Spots North Spots South 
of the Equator of the Equator Mean 
= A Rhee , Mean Distance 
Mean I Latitude from 
Days Mean Helio- age i of Entire Equator 
Photographed graphic i Spotted of all 
° Latitude Area Spots 





° 


1944 366 19°00 22°81 ~8-70 21°53 


Old Cycle 4°17 7°61 ~1°14 5°72 
1944 3664 YS, Cycle 35 22°18 77 24°02 —9°63 23°45 
1945 365 121 20°13 309 20°26 —8-92 20°22 
1946 365 1127 20°74 690 18°79 +5°73 20°00 
1947 = 365 992 16°58 1645 17°86 —4°91 17°38 


The mean daily spot area of 2637 millionths for the year 1947 (see Table I) 
may be compared with those of the maximum years of preceding cycles as 
follows :— 

Mean daily 
Year spot area 


1947 2037 
1937 2074 
1928 1390 
1917 1537 
1905 IIgI 
1893 14604 
1883 1154 


The ratio of faculae/whole spots for the year 1947 is 1-10, the corresponding 
ratio for the year 1946 being 1-20. ‘These values are both much below average, 
the mean value for the peak years of the six preceding cycles being 1-75. 

Faculae associated with individual spot groups were examined for a number 
of sample years near earlier maxima and it was found that the area-ratios, 
faculae/spots gave a statistical decrease from the smaller to the larger spot 
groups. Thus, the mean ratio for spots around 1000 millionths of the Sun’s 
hemisphere is about one-half the value for groups about 500 units, and double 
that for spots of 1500 units and upwards. 

It would appear that the occurrence of several unusually large sunspots 
(including the three largest on record) during 1946 and 1947 with their 
individually small values of the ratio faculae/spots, tended to diminish the 
mean ratio for these two years. 

The stability of the measurement of faculae was tested by remeasurement 
of 72 plates throughout the year and found to be satisfactory. The possibility 
of an increased proportion of poor plates on which the definition might smudge 
out fainter areas is also not borne out by a general examination of the plates. 
It must be concluded therefore that the major factor producing the decreased 
ratios was an intrinsic, smaller increase in the amount of faculae disproportionate 
to the abnormal increase of sunspots in 1946 and 1947. 

(2) On no day in 1947 were sunspots and faculae absent from the Sun’s disk. 
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(3) The mean daily area on the southern hemisphere was considerably in 
excess of that on the northern. But of the 653 units of the southern excess, 
336 units were alone due to the two disk passages of the giant groups referred 
to above. In actual numbers of groups, the southern was the more prolific, 
as shown next. 

(4) The number and distribution, northern and southern hemispheres, of 
spot groups of (a) two days’ duration or longer; (5) one day only, are as 


follows :— 
(a) (6) 
Northern spots 250 52 
Southern spots 304 59 
Total 554 III 


The range in latitude for the longer-lived spots was from 1° to 42°. The 
latitudes of the 1-day spots ranged from 2° to 49°. 

(5) The mean weighted latitude of all spots was 17°-4. The latitudes of the 
two giant spots were —23°:1 and —24°-4 respectively. 

A short descriptive survey of the sunspots of 1947 has been given in 
Monthly Notices, 108, 122, 1948. 

An appended table gives the mean daily areas of sunspots and of faculae 
(corrected for foreshortening and expressed in millionths of the Sun’s 
hemisphere) for each calendar month of 1947. 


Monthly Mean Daily Areas of Sunspots and Faculae 
1947 

Month Spots Faculae Month Spots Faculae 
Jan. 2308 2984 July 2451 3402 
Feb. 2721 2526 Aug. 3214 3594 
Mar. 3596 2394 Sept. 1992 3091 
Apr. 3950 2507 Oct. 2380 2890 
May 3741 3182 Nov. 1686 2209 
June 1934 3239 Dec. 1643 2663 


Royal Greenwich Observatory, 
Herstmonceux Castle, Sussex: 
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